Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 29.8% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC29.8% 317.2 / 156.7 188intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (69/63);  new base C (49/7);  total (118/70)
Fisher's exact test for biased strand distribution p-value = 2.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 190 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082
                                                                                                   |                                                                                                        
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:2739913/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:2684643/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                        >  1:2990138/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      <  1:2782064/100‑1 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      >  2:1736738/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                    <  1:2865785/100‑1 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                    <  2:382320/100‑1 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTATGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                   <  1:1018245/100‑1 (MQ=37)
        gacggCGCTGCACGCAGCGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  1:3324301/97‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                 <  2:602555/100‑3 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                 <  2:3297967/100‑3 (MQ=37)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                           <  2:2890647/100‑1 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑aa                                                                                             <  1:420566/99‑2 (MQ=25)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aa                                                                                          >  2:3311698/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                         >  2:1212625/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                           >  1:918624/1‑98 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                      >  1:3248448/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                     >  2:1939794/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                     <  2:276078/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                   <  2:3148290/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                  >  1:3065622/1‑100 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                  <  1:2910478/100‑1 (MQ=255)
                       agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaa                                                                                <  2:2808960/100‑1 (MQ=255)
                        gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaa                                                                                <  2:2148281/99‑1 (MQ=255)
                        gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaaa                                                                               >  2:512624/1‑100 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                             <  2:2183261/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                            <  1:1939784/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                            <  2:148912/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                            >  2:1081574/1‑99 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                    >  2:1987695/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                   >  2:2532193/1‑100 (MQ=255)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                  >  1:3114958/1‑100 (MQ=255)
                                           tataCAGATTTTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                    <  1:512622/92‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑t                                                          >  1:1918168/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                                    >  2:2307703/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGcc                                              >  2:1767517/1‑100 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCt                                             <  1:2013829/100‑1 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTg                                            >  1:75724/1‑100 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTg                                            <  1:2940/100‑1 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTg                                            <  2:1353975/100‑1 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGa                                           >  1:3534522/1‑100 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGAt                                          >  2:1639038/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCTTTCTCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                                         <  1:2492857/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                                         <  2:2271243/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                      <  2:3050567/100‑1 (MQ=255)
                                                             aGTGCGCTTTGGTTATGCCGTATGCGGCGTTAACGCCTTAACCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAg                                       <  1:2134437/99‑1 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                     >  2:3193447/1‑100 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                    <  2:2713378/100‑1 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                  >  1:599879/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                  <  2:1789689/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                 <  1:2422677/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAAACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                 >  2:2501397/1‑100 (MQ=255)
                                                                  gcTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGTCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCATA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                 <  1:39616/100‑1 (MQ=255)
                                                                   tgatttgtaggCCGGATAACGCATTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGTCCTGATAAGCGTAgcgc                               >  1:1612746/11‑100 (MQ=18)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                >  2:72301/1‑100 (MQ=255)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGAGAAGCGGAGCGCAt                             >  2:2066863/1‑100 (MQ=255)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCAGAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                            >  2:3109539/7‑100 (MQ=17)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                            <  1:1136920/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                           <  2:3458914/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                           <  2:1421305/100‑1 (MQ=255)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                          <  2:2811572/100‑1 (MQ=255)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                          <  1:1949751/100‑1 (MQ=255)
                                                                          taggCCGGATAAGGCGTTTACGTCGCATCCGGCCT‑‑AC‑AAAAGAGTGCAAATTCAATACATTGCATG‑GG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  1:1260320/4‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCAAG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                         >  1:3007797/4‑100 (MQ=17)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                         <  2:2216283/100‑1 (MQ=255)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  1:2724005/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCAAT‑CT‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:955032/3‑100 (MQ=18)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATCCAATACATCGCAGA‑AT‑TTTC‑‑GTAGGCCTGACAAGCGGAGCGCATCAGGCa                       >  1:694511/1‑100 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:2903438/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCATG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:623107/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:3343404/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:2520086/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:2990578/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCAGA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:1240649/2‑100 (MQ=17)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAACTTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                        >  2:3389846/1‑98 (MQ=255)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  1:1510454/1‑100 (MQ=255)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  1:1659209/1‑100 (MQ=255)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  1:781198/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  1:1808540/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  1:2365248/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  2:1424177/1‑100 (MQ=21)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                     <  2:2101170/92‑3 (MQ=18)
                                                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                    >  1:666677/8‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                   >  2:437911/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                   >  2:421041/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                   >  2:870255/7‑96 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATATATTGCAGG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  1:1997370/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                   >  1:2603507/7‑96 (MQ=18)
                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                   >  2:1332835/1‑100 (MQ=255)
                                                                                 gataaggcgTTTACGCCGCATCCGGCATAAAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                    >  1:3251720/6‑100 (MQ=14)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATCAATTGCAGG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  2:849191/6‑99 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  1:1101177/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct              >  2:1382781/1‑100 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct              <  2:2557941/100‑1 (MQ=255)
                                                                                         tctAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          <  1:1201389/98‑1 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa        >  1:2317624/1‑100 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt       >  1:1349859/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      >  2:940539/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGCCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATc      <  2:1003516/100‑1 (MQ=255)
                                                                                              cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa     >  2:1561073/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑