Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC41.6% 273.5 / inf 281intergenic (‑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 (92/72);  new base C (113/4);  total (205/76)
Fisher's exact test for biased strand distribution p-value = 4.41e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-02

Only 100 of 283 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                               
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                                >  1:2683104/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTGTGCCGGATGCGGCGTTAACGCCTTAt                                                                                                               >  1:1380772/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                               >  2:2130855/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                             >  2:1568755/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                         <  2:2060411/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                        >  2:1633559/1‑100 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                       >  2:1334344/1‑100 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                      <  2:37798/100‑1 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                      <  2:34747/100‑1 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                                      <  1:2187116/100‑2 (MQ=25)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                    >  1:3175436/1‑100 (MQ=255)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                                   <  1:1932104/100‑1 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                                  >  1:834200/1‑100 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                               >  2:2965499/1‑100 (MQ=255)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                             <  1:1218372/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                            <  2:128274/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                            <  1:2906089/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                           <  1:1050419/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                          >  2:821864/1‑100 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                          >  2:1303478/1‑100 (MQ=255)
                       agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                                         >  1:999389/1‑100 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                                      >  1:3233446/1‑100 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                     <  1:2840855/100‑1 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                             >  1:2396575/1‑100 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                                     <  2:315224/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                                     >  2:738546/1‑100 (MQ=255)
                                          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑At                                                                 <  2:1533940/100‑1 (MQ=255)
                                          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑At                                                                 <  2:2802759/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑g                                                      >  2:1735310/1‑100 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                                     <  1:662644/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                                     >  2:3192613/1‑100 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                                     >  1:879260/1‑100 (MQ=255)
                                                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGc                                                 <  2:613185/100‑1 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGccn                                               <  2:1170729/100‑2 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           >  1:2090669/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  2:1065198/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                        <  2:3245923/100‑1 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                       >  2:4128/1‑100 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                       >  1:589677/1‑100 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  2:2390336/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     >  2:1012508/1‑100 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     >  1:381876/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:2996701/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  1:2632174/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  1:2258854/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:2770920/100‑1 (MQ=255)
                                                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  2:510925/8‑100 (MQ=18)
                                                                      tGTTTATGCCGGATGGGGCGTTAGCGCCTTATCCGGTCTAC‑AAGATCGTGCAAATTCGATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  2:254476/99‑1 (MQ=18)
                                                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  2:1121180/7‑100 (MQ=18)
                                                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:1181973/7‑100 (MQ=18)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑C‑‑A‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  2:2345972/7‑100 (MQ=17)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             <  1:1509871/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             <  1:2765394/100‑1 (MQ=255)
                                                                        ngtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTG‑C‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  2:696428/6‑100 (MQ=18)
                                                                        cgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTG‑T‑‑G‑GG‑CT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  2:2242991/6‑100 (MQ=21)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTG‑C‑‑ATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:3022136/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑C‑‑A‑GA‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:780444/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTG‑CGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:794244/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTG‑C‑‑G‑GA‑GA‑A‑‑‑AAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:1477032/5‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:2283585/4‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:1747510/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTG‑CGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:200442/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                           >  2:2439937/4‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1661117/3‑100 (MQ=21)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑C‑‑AGGG‑AC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  2:2708847/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTG‑C‑‑A‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                         >  2:2281968/3‑99 (MQ=18)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑GTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                         <  1:1144524/100‑1 (MQ=255)
                                                                            ggacggataaggcgTTCACGCCGCATCCGTCCTAC‑AAGTTCATGCAATTTTAATAAGTTGCC‑‑A‑GA‑TA‑T‑‑‑CCT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2423221/11‑100 (MQ=25)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑C‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2650298/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑C‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3195863/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2169263/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑C‑‑A‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                         >  2:1777595/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:625132/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑C‑‑A‑GG‑GA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                         >  1:2664858/2‑100 (MQ=18)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:2120521/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑‑‑‑A‑AT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  1:1821118/1‑98 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTG‑C‑‑‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                         >  1:3137849/1‑100 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1724615/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑C‑‑A‑GG‑GA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                        >  1:3089461/1‑99 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑C‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:1253716/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTG‑T‑‑ATGG‑AA‑A‑‑‑TTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGCtttt                      >  2:479706/9‑100 (MQ=25)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:2720228/1‑100 (MQ=255)
                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑C‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  2:1371519/1‑98 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:2919671/7‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑C‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1488497/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:2831307/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                     >  2:184925/7‑96 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑C‑‑A‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:1425824/7‑100 (MQ=17)
                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  1:497200/1‑100 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:1159230/100‑1 (MQ=255)
                                                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑C‑‑AGGG‑AC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctttac                  >  2:2160413/4‑93 (MQ=21)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                <  1:53518/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:367114/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:1448762/100‑1 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:1576267/100‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:725969/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  2:1697269/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  1:2512069/1‑100 (MQ=255)
                                                                                                 cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc    >  1:2461578/1‑100 (MQ=255)
                                                                                                   tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:3223255/1‑100 (MQ=255)
                                                                                                   |                                                                                                               
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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