Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC31.9% 389.1 / 262.1 278intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (80/84);  new base C (76/1);  total (192/86)
Fisher's exact test for biased strand distribution p-value = 3.32e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-03

Only 100 of 280 total aligned reads displayed.

AGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705881‑3706077
                                                                                                 |                                                                                                       
 gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                      >  2:561208/1‑100 (MQ=255)
     tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    >  2:2904049/1‑98 (MQ=255)
      gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                 >  2:1411145/1‑100 (MQ=255)
      gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgca                                                                                                 <  2:1042163/100‑4 (MQ=38)
           cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                            >  2:2611342/1‑100 (MQ=255)
            ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                           >  2:870413/1‑100 (MQ=255)
            ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                           <  2:2753336/100‑1 (MQ=37)
             gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:3055624/100‑1 (MQ=255)
             gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:1350471/100‑1 (MQ=255)
               tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                        <  1:1215384/100‑1 (MQ=37)
                ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                       <  1:1010223/100‑2 (MQ=25)
                  cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                     >  1:1393163/1‑100 (MQ=255)
                   gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                    >  2:3200668/1‑100 (MQ=255)
                    cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                   >  2:2927379/1‑100 (MQ=255)
                        aCGGAGAGGCACAGCGCCACCGAGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                 <  1:2413519/98‑1 (MQ=255)
                        aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                               <  2:1846455/100‑1 (MQ=255)
                        aCGCAGAGGCACAGCGCCACCGGGTATACAGATTCTTGACAAAGTGCGATTTGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                               <  1:2680874/100‑1 (MQ=255)
                         cGCCGAGGAACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                              <  2:2075076/100‑1 (MQ=255)
                         cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                              >  1:36698/1‑100 (MQ=255)
                           cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                            >  2:1385114/1‑100 (MQ=255)
                           cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                            >  1:2811676/1‑100 (MQ=255)
                             gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                          <  1:1108216/100‑1 (MQ=255)
                                 cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                      >  2:2206902/1‑100 (MQ=255)
                                   caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                    >  2:657288/1‑100 (MQ=255)
                                   caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                    <  2:414607/100‑1 (MQ=255)
                                    aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                   <  1:2702218/100‑1 (MQ=255)
                                        ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                               <  1:2178248/100‑1 (MQ=255)
                                          aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                             <  2:780017/100‑1 (MQ=255)
                                            cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑g                                                          >  1:2561218/1‑100 (MQ=255)
                                                tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t                                                     <  1:3343268/100‑1 (MQ=255)
                                                tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t                                                     >  2:1016478/1‑100 (MQ=255)
                                                    cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTc                                                  >  1:1693706/1‑99 (MQ=255)
                                                    cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                               <  2:1099787/100‑1 (MQ=255)
                                                    cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                               <  2:2147790/100‑1 (MQ=255)
                                                      gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTa                                             >  1:1985488/1‑100 (MQ=255)
                                                       ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAg                                            >  2:1976259/8‑100 (MQ=18)
                                                       aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAg                                            >  2:1950292/1‑100 (MQ=255)
                                                             gACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTg                                       >  2:511868/1‑99 (MQ=255)
                                                              aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGAt                                     <  1:377126/100‑1 (MQ=255)
                                                               cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATa                                    >  2:1754236/1‑100 (MQ=255)
                                                                   gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                                <  2:2119197/100‑1 (MQ=255)
                                                                   gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                                <  1:1888560/100‑1 (MQ=255)
                                                                     gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGa                              >  1:173554/1‑100 (MQ=255)
                                                                     gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGa                              >  2:1683094/1‑100 (MQ=255)
                                                                        cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                               >  1:2922746/1‑96 (MQ=255)
                                                                        cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                           <  2:200376/100‑1 (MQ=255)
                                                                         tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                          <  1:511862/100‑1 (MQ=255)
                                                                         attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                          >  2:1583615/10‑100 (MQ=18)
                                                                           tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                       >  1:1584114/8‑100 (MQ=17)
                                                                           tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                        >  2:1884199/8‑100 (MQ=21)
                                                                             tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  2:3032432/1‑100 (MQ=255)
                                                                              ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     >  1:2256146/1‑100 (MQ=255)
                                                                              ttATGCCGGAGGCGGCGTTAACGCCTTGTCCGGCCTACAAAATCGTGCAAATTGAATACATTGCA‑GAAT‑GTTC‑‑GTAGGCCTGAGAAGCGGAGCGCATCAg                     >  1:2027701/1‑100 (MQ=18)
                                                                              gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:2674041/5‑100 (MQ=17)
                                                                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:2786412/5‑100 (MQ=21)
                                                                               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                    >  2:1173356/1‑100 (MQ=255)
                                                                               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                    >  2:2713442/1‑100 (MQ=255)
                                                                                aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:2539254/3‑100 (MQ=18)
                                                                                aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:3042451/3‑100 (MQ=17)
                                                                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:2226113/3‑100 (MQ=21)
                                                                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:2502633/3‑100 (MQ=18)
                                                                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1230055/3‑100 (MQ=21)
                                                                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:3443466/3‑100 (MQ=21)
                                                                                aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑AAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                   >  2:2437221/1‑100 (MQ=255)
                                                                                aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc                   <  2:3233774/100‑1 (MQ=21)
                                                                                 ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2732408/11‑100 (MQ=21)
                                                                                 ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                 >  2:394775/2‑100 (MQ=18)
                                                                                 ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:246317/2‑100 (MQ=18)
                                                                                 ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑TGGG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1145500/2‑100 (MQ=17)
                                                                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCA‑GGAT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:843489/2‑100 (MQ=18)
                                                                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGACTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGAGCATCAGGCa                  >  1:2289709/2‑100 (MQ=18)
                                                                                 ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  2:449092/2‑100 (MQ=18)
                                                                                  gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  1:218837/1‑99 (MQ=18)
                                                                                  gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑AAT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgagat                 >  2:3332629/1‑96 (MQ=18)
                                                                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                >  1:102586/1‑100 (MQ=21)
                                                                                  gCCGGATAAGGCGGTCACGCCGCATCCGGCCTACGAAACCGCAAGAATTCAATATATTGCA‑GGAG‑CGGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  2:2586407/1‑99 (MQ=18)
                                                                                   ctggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCA‑CAGC‑TAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                >  2:59344/9‑98 (MQ=25)
                                                                                   ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt               >  2:3093070/9‑100 (MQ=18)
                                                                                   ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:238962/9‑98 (MQ=18)
                                                                                   ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  1:1359262/9‑98 (MQ=18)
                                                                                   ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:2280278/9‑98 (MQ=18)
                                                                                   ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                >  1:2781434/9‑98 (MQ=18)
                                                                                   ccGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  2:958944/1‑100 (MQ=18)
                                                                                    cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTACA‑GAGT‑TCAC‑‑GTAGGACTGATAAGCGTAGCGCATCAGGCagtt               >  1:2069487/8‑97 (MQ=17)
                                                                                    cggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:747371/8‑97 (MQ=21)
                                                                                     ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCCATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:1543336/7‑100 (MQ=17)
                                                                                     ggataaggcgTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:923982/7‑100 (MQ=18)
                                                                                     ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct              >  2:2864781/7‑96 (MQ=17)
                                                                                     ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagct               >  1:2986906/7‑96 (MQ=18)
                                                                                     ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt              >  2:2951293/1‑100 (MQ=255)
                                                                                     ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtgt               >  2:3478379/1‑97 (MQ=255)
                                                                                      gataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt             >  2:1538955/6‑95 (MQ=21)
                                                                                      gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt             >  2:3387718/6‑95 (MQ=21)
                                                                                      gataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCA‑GGAA‑CGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactac             >  1:720422/6‑95 (MQ=25)
                                                                                       ataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg            >  1:642990/5‑100 (MQ=17)
                                                                                         gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          <  2:2561252/100‑1 (MQ=255)
                                                                                           ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGAGAATCAGGCATTtttgc           >  1:3231777/1‑97 (MQ=255)
                                                                                            gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       <  2:2314379/100‑1 (MQ=255)
                                                                                            gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATGCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       <  1:3290504/100‑1 (MQ=255)
                                                                                             cgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GAGA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAATtttgcgtttg      >  1:131119/1‑100 (MQ=17)
                                                                                                 |                                                                                                       
AGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705881‑3706077

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: ‑