Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,258,963 +C 100% coding (778/1971 nt) efeU → ferrous iron permease EfeU

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,258,9631.→C100.0% 38.9 174coding (778/1971 nt)efeUferrous iron permease EfeU
Reads supporting (aligned to +/- strand):  new base (82/78):  ref base (0/14):  total (82/92)
Fisher's exact test strand distribution p-value = 9.47e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.30e-01

Only 100 of 189 total aligned reads displayed.

CAATGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.A.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGCCTCCTT  >  NC_014034/3258863‑3259058
                                                                                                       |                                                                                               
  aTGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.A.cc                                                                                               >  2:3398697/1‑101 (MQ=255)
      cTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.A.ccctgc                                                                                           >  1:2906663/1‑97 (MQ=255)
       ttggTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCCTGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.A.ccctgcc                                                                                          <  1:4183681/99‑6 (MQ=255)
            cgcgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGctct                                                                                      <  2:4122425/101‑1 (MQ=255)
              cgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTc                                                                                    <  2:1782467/101‑1 (MQ=255)
               gggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTcc                                                                                   >  1:3435996/1‑101 (MQ=255)
               gggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTcc                                                                                   <  1:1739653/101‑1 (MQ=255)
                ggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCt                                                                                  <  1:1094615/101‑1 (MQ=255)
                ggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCt                                                                                  <  1:4140199/101‑1 (MQ=255)
                ggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCt                                                                                  <  2:2651388/101‑1 (MQ=255)
                ggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCt                                                                                  >  1:1704446/1‑101 (MQ=255)
                 gCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTg                                                                                 <  1:3936289/101‑1 (MQ=255)
                 gCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTg                                                                                 <  1:1688601/101‑1 (MQ=255)
                 gCACCGTGCTTTCGGCCCTCTTTGGCTATCCCGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTg                                                                                 <  2:524812/101‑1 (MQ=255)
                  cACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTgc                                                                                >  1:70089/1‑101 (MQ=255)
                     cGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGcc                                                                             <  1:4060196/101‑1 (MQ=255)
                      gTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGccc                                                                            >  1:2793866/1‑101 (MQ=255)
                      atcttttCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGccc                                                                            <  1:186290/97‑1 (MQ=255)
                        gCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcg                                                                          >  1:3255917/1‑101 (MQ=255)
                        gCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcg                                                                          <  1:3613685/101‑1 (MQ=255)
                         cTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcgc                                                                         >  2:1354594/1‑101 (MQ=255)
                         cTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcgc                                                                         >  1:2713727/1‑101 (MQ=255)
                         cTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcgc                                                                         <  1:2394226/101‑1 (MQ=255)
                         cTTTCGGCGCTCTTTGGCTATCACGACACCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCcgc                                                                         <  2:1981992/101‑1 (MQ=255)
                          tttCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGcccagat                                                                        >  2:186252/1‑97 (MQ=255)
                           ttCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCgcgc                                                                       <  1:1985818/101‑1 (MQ=255)
                             cGGCGCTTTTTGACAATCACGTCAGCCGGGCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgccc                                                                     <  2:1995276/101‑1 (MQ=255)
                             cGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgccc                                                                     >  2:1899060/1‑101 (MQ=255)
                             cGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgccc                                                                     <  2:4033640/101‑1 (MQ=255)
                              ggCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgcccg                                                                    >  2:2873700/1‑101 (MQ=255)
                               gcgcTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgcccgc                                                                   <  2:2363612/101‑1 (MQ=255)
                               gcgcTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGcgcccgc                                                                   <  1:1696562/101‑1 (MQ=255)
                                 gcTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCgcccgccc                                                                 >  1:3205540/1‑101 (MQ=255)
                                   tctTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGccccc                                                               >  2:1964954/1‑101 (MQ=255)
                                    ctTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGcccccc                                                              <  2:3842670/101‑1 (MQ=255)
                                     tttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCt                                                             <  1:939639/101‑1 (MQ=255)
                                     tttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCt                                                             >  2:1660828/1‑101 (MQ=255)
                                      ttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCtt                                                            <  1:3821496/101‑1 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           >  2:3032578/1‑101 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           >  1:2043332/1‑101 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           >  1:2334480/1‑101 (MQ=255)
                                        ggCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTcc                                                          <  1:904490/101‑1 (MQ=255)
                                          cTATCACGACAGCCCGCCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGc                                                        <  2:76074/101‑1 (MQ=255)
                                          cTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGc                                                        >  2:3085994/1‑101 (MQ=255)
                                          cTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTCCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGc                                                        >  1:1906898/1‑101 (MQ=255)
                                           tATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCa                                                       <  1:3351259/101‑1 (MQ=255)
                                             tCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCaaa                                                     <  1:48508/101‑1 (MQ=255)
                                             tCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCaaa                                                     <  2:437978/101‑1 (MQ=255)
                                              cACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCaaaa                                                    >  2:1719471/1‑101 (MQ=255)
                                              cACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCaaaa                                                    >  1:2206742/1‑101 (MQ=255)
                                              cACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCCCGTTCCTGGCC.GTG.ACCCCGCCGCTCTTACTGCGCCCCGCCCCCGCCCCCCTTCCGCaaaa                                                    >  1:3856688/1‑101 (MQ=255)
                                               ccGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAg                                                   <  2:2325913/100‑1 (MQ=255)
                                                cGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAgg                                                  <  2:2575320/101‑1 (MQ=255)
                                                cGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAgg                                                  <  2:2210541/101‑1 (MQ=255)
                                                 gACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGc                                                 >  1:3690492/1‑101 (MQ=255)
                                                   cggCCCGACCCAGGCCGGGGTCGCCCTCTGGGCCGCGTTCCTGGCC.TTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCt                                               <  2:3516910/99‑1 (MQ=255)
                                                   cAGCCCGCCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC..TG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTg                                              <  2:2813453/101‑1 (MQ=255)
                                                   cAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCt                                               <  1:3534091/101‑1 (MQ=255)
                                                     gCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGa                                             <  2:939300/101‑1 (MQ=255)
                                                      cccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGAt                                            <  1:2382152/101‑1 (MQ=255)
                                                       ccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATc                                           <  2:2136890/101‑1 (MQ=255)
                                                         gACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCg                                         <  1:1365886/101‑1 (MQ=255)
                                                          ccccAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGa                                        <  1:1492796/100‑1 (MQ=255)
                                                          aCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGa                                        >  1:773229/1‑101 (MQ=255)
                                                          aCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGa                                        >  2:1520622/1‑101 (MQ=255)
                                                           cccAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGAt                                       <  2:1544555/101‑1 (MQ=255)
                                                           cccAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGAt                                       <  1:2885178/101‑1 (MQ=255)
                                                             cggggCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGa                                     <  1:2558079/99‑1 (MQ=255)
                                                             cAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGa                                     >  1:2847490/1‑101 (MQ=255)
                                                             cAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGa                                     >  1:572816/1‑101 (MQ=255)
                                                              aGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAc                                    >  2:2935621/1‑101 (MQ=255)
                                                                ggCGAGGTCGCCCTCGGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAccc                                  <  2:70074/101‑1 (MQ=255)
                                                                 gCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAcccc                                 <  1:1696434/101‑1 (MQ=255)
                                                                  aagaggcGCCCGTTGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCt                                <  2:811303/95‑1 (MQ=255)
                                                                   gAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCtt                               <  1:1775772/101‑1 (MQ=255)
                                                                   gAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCtt                               >  1:1970605/1‑101 (MQ=255)
                                                                    aGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTc                              >  1:4236260/1‑101 (MQ=255)
                                                                     ggTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTcc                             <  1:1660943/101‑1 (MQ=255)
                                                                     ggTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTcc                             <  1:724823/101‑1 (MQ=255)
                                                                      gccccccTCTGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCg                            <  1:1886472/97‑1 (MQ=255)
                                                                      gTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCg                            <  2:2474981/101‑1 (MQ=255)
                                                                      gTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCg                            >  1:2681925/1‑101 (MQ=255)
                                                                       tCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGc                           >  2:1127675/1‑101 (MQ=255)
                                                                       tCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGc                           >  1:3885120/1‑101 (MQ=255)
                                                                         gCCCTCTGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGcac                         <  1:1048931/101‑1 (MQ=255)
                                                                         gCCCTCTGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGcac                         >  1:2474958/1‑101 (MQ=255)
                                                                            ctctGGGCCGCGTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAg                      <  2:3356057/101‑1 (MQ=255)
                                                                              ctGGGCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAgcg                    <  1:3032803/101‑1 (MQ=255)
                                                                               tGGGCGGCGTTCCTGGCC.TTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGAGGACCCCTTCCGCACAAgcgc                   <  2:125165/101‑1 (MQ=255)
                                                                                gggCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGcc                  <  1:1312689/101‑1 (MQ=255)
                                                                                 ggCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCg                 >  1:1011642/1‑101 (MQ=255)
                                                                                  gCCGCGTTCCTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGt                <  2:1053460/101‑1 (MQ=255)
                                                                                    cgcgTTCCTGGCC.GTGCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTcc              <  1:2122038/101‑1 (MQ=255)
                                                                                       gTTCCTGGCC.GTGCC.CCGGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCggg           <  2:3804083/101‑1 (MQ=255)
                                                                                         tCCTGGCCGGGG.C.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGgcg         <  2:2813753/101‑1 (MQ=255)
                                                                                          ccTGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGgcgc        >  2:1366380/1‑101 (MQ=255)
                                                                                            tGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGCTCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcct      >  1:1935976/1‑101 (MQ=255)
                                                                                            tGGCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcct      >  2:1493079/1‑101 (MQ=255)
                                                                                             ggCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcctc     >  1:3747822/1‑101 (MQ=255)
                                                                                              gCC.GTG.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcctcc    >  1:1167971/1‑101 (MQ=255)
                                                                                                       |                                                                                               
CAATGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCC.GTG.A.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGCCTCCTT  >  NC_014034/3258863‑3259058

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG