Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 36.0% 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,9850TC36.0% 213.5 / 170.3 164intergenic (‑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 (49/56);  new base C (53/6);  total (102/62)
Fisher's exact test for biased strand distribution p-value = 1.50e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 166 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083
                                                                                                  |                                                                                                      
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                       >  1:771350/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                     >  1:1062191/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                    <  2:1022773/100‑1 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                    >  1:1707957/1‑100 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                  <  2:1707948/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                  >  1:1778813/1‑100 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                  <  2:2580812/100‑1 (MQ=37)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                 <  1:48697/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                <  2:923265/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                <  2:544539/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                <  2:369294/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                <  2:1216467/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                <  1:2325448/100‑2 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                               <  2:2298309/100‑3 (MQ=37)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                             <  1:2336721/99‑1 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                            <  1:1607919/100‑3 (MQ=25)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                           <  1:805886/100‑7 (MQ=25)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                        <  2:1007152/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                     >  2:876076/1‑100 (MQ=255)
                   cGCAGAGGCACAGGGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                    <  2:902983/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                    <  1:2418044/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                    >  2:757397/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                   <  1:1145859/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                  >  1:815253/1‑100 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                  <  1:2553411/100‑1 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                <  1:1614245/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                             <  1:496457/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            >  1:540830/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                     <  1:2378914/100‑1 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                     >  1:419434/1‑100 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACCAAATCGTGAAAATTccctacct                                                                      >  1:1516866/1‑91 (MQ=17)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt                                                             >  1:1127505/1‑100 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t                                                           >  1:2007420/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑ttt                                                         >  2:2564313/1‑100 (MQ=255)
                                            cttccgatcttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑ttt                                                         <  2:41584/91‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                     <  2:698255/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                     <  2:1062186/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                     >  2:1442099/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                     >  2:2053708/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑Gt                                                    >  2:1505116/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTa                                                   <  1:2456670/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTa                                                   <  1:1295404/100‑1 (MQ=255)
                                                   cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGc                                                >  2:1375860/6‑100 (MQ=18)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGa                                            >  2:458724/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGa                                            >  2:616009/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGa                                            >  1:1376531/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGAt                                           >  2:857977/1‑100 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATaa                                         <  2:1954618/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                                      <  1:887673/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                                     <  1:697698/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                  >  2:697753/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                 <  2:1215278/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  2:166188/1‑100 (MQ=255)
                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGCCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAagcgcagcgc                                >  2:1151671/1‑100 (MQ=255)
                                                                   cttttgtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAgcgc                                >  2:337654/10‑100 (MQ=11)
                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                              >  2:782768/8‑100 (MQ=17)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            <  2:469172/94‑1 (MQ=21)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCA                            <  1:1824567‑M1/100‑1 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:1443529/6‑100 (MQ=17)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAGTCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAG                           <  2:1352379‑M1/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAG                           <  1:2178491‑M1/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAG                           >  2:652013‑M1/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:2460344/5‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:2246027/5‑100 (MQ=11)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAg                           >  1:385782/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:1947520/4‑100 (MQ=25)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1432070/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCCATATATTGCA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:711191/3‑100 (MQ=11)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1474869/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:557520/3‑100 (MQ=17)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  1:1617054/1‑100 (MQ=255)
                                                                           ggacggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:296812/11‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:411202/2‑100 (MQ=11)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1317099/2‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:2300905/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  1:1305373/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1915451/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGAAGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGTA‑TGAG‑CACC‑‑GTAGGACTAATAAGCGTAGAGCATc                             >  2:732663/1‑94 (MQ=12)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                     >  2:2342837/9‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                     >  1:2139329/9‑100 (MQ=21)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1898047/1‑100 (MQ=255)
                                                                               ngataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:2597082/7‑98 (MQ=21)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  2:1764999/7‑98 (MQ=21)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2450120/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1001213/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagct                    >  1:220144/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                    >  1:1988629/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2474515/7‑100 (MQ=18)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                   >  1:871028/6‑95 (MQ=25)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1177594/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:1697993/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  2:190575/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  2:608934/100‑1 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt             >  1:332203/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1493987/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  2:1919742/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  1:262755/1‑100 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtn       <  1:2427260/99‑2 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:626750/100‑1 (MQ=255)
                                                                                                 ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  >  1:734498/1‑100 (MQ=255)
                                                                                                  |                                                                                                      
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083

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