Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,640,9810CT52.9% 14.3 / 102.0 85F40F (TTC→TTT) dtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (4/36);  new base T (45/0);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 7.34e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-06

GGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTG  >  NC_000913/3640841‑3641131
                                                                                                                                             |                                                                                                                                                         
ggTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTtaaa                                                                                                                                                 >  1:1204990/1‑150 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGGTTaaca                                                                                                                                                 >  2:1307563/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTTTTTAACa                                                                                                                                                >  1:582995/1‑151 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTtaaaa                                                                                                                                                >  2:1126606/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTtaaaa                                                                                                                                                >  1:1405197/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACa                                                                                                                                                >  2:708896/1‑151 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAAc                                                                                                                                                 <  1:323419/150‑1 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGttttttttttttaaa                                                                                                                                                >  2:343093/1‑142 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTttcttttcttgaca                                                                                                                                                >  2:832906/1‑142 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTtaaaa                                                                                                                                                >  1:485497/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTtaaaa                                                                                                                                                >  1:1442412/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGGTaaaaa                                                                                                                                                >  1:1479498/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTCTTTGTTtaaaa                                                                                                                                                >  2:1616074/1‑149 (MQ=255)
   aTTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTTTTTAAACa                                                                                                                                                >  2:179721/1‑150 (MQ=255)
        aaGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGa                                                                                                                                         <  2:841512/151‑1 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTTTTTAAAAACCTTGGttt                                                                                                                                       >  1:1236205/1‑151 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAGCTTGGGttt                                                                                                                                      >  1:278190/1‑150 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTTTtaaaacccttggttt                                                                                                                                      >  1:675858/1‑139 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTtaaaacactgggttt                                                                                                                                      >  1:1397689/1‑139 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAACTTGGCttt                                                                                                                                      >  2:486691/1‑150 (MQ=255)
           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTTAACAGATTgggatt                                                                                                                                      >  2:855136/1‑147 (MQ=255)
              aataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTTAAAAACtttggttttc                                                                                                                                    >  1:1538244/1‑142 (MQ=255)
                  aataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAag                                                                                                                               <  1:1258594/151‑1 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTTTTTAACAGCTTGGGTTTTCGCCagag                                                                                                                             >  1:1521349/1‑151 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAGCTTTGGTTTTCGCAagag                                                                                                                             >  2:791391/1‑151 (MQ=255)
                     atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGc                                                                                                                            <  2:1306058/151‑1 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAACACTTGG‑TTTTCGCCAGAGCAg                                                                                                                          >  2:614664/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGGTTAAAAACTTGGGTTTTCGccaaaggag                                                                                                                          >  1:1395364/1‑143 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTTTTAACAACTTGGGTTTTCCCCAAAGCAg                                                                                                                          >  1:839430/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAACCTGGATTTTCGCCagagacg                                                                                                                          >  2:1499907/1‑148 (MQ=255)
                          tACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCt                                                                                                                       <  1:99498/151‑1 (MQ=255)
                          tACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCt                                                                                                                       <  1:315529/151‑1 (MQ=255)
                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAGCCTTGGTTTTCGccaaaggcgggttt                                                                                                                     >  2:669481/1‑138 (MQ=255)
                              acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTg                                                                                                                   <  1:1180911/151‑1 (MQ=255)
                              acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGCTTTTTTTGTTTAACAGCTTGGGTTCTCGCCAGAGCAGGCgtttt                                                                                                                   >  2:921621/1‑150 (MQ=255)
                                aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTc                                                                                                                 <  2:1131153/151‑1 (MQ=255)
                                aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTc                                                                                                                 <  2:111230/151‑1 (MQ=255)
                                  cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTTTTTAACAGCTTTGGTTCTTGCCAGAGGAGGGtttttttc                                                                                                                >  1:958908/1‑146 (MQ=255)
                                  cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAAAAACCTGGGTTTTCGCAAAAGCAGGGTTTTTTcc                                                                                                                >  1:1192356/1‑149 (MQ=255)
                                  cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTTGTTAAACAGCTTGGGTTTTCTCCCAAGCAGGGTTTTTTcc                                                                                                                >  2:28096/1‑149 (MQ=255)
                                   acCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTTAAAAACTTGGGTTTTCGCCAAAGCAGGCCTTTTTct                                                                                                                >  1:1421603/1‑148 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTCCTGGC‑GGGTTTTTTGGTTATAAATCTTTGGTTCTCGCCAAAGCCGGCCTTTTTCGCtt                                                                                                             >  2:530896/1‑151 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAACTTGGGTTTTCGCCAAAGCAGGGTTTTTTCACtt                                                                                                             >  1:145409/1‑151 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAGCTTTGGTTTTCGCAAAAGGCGGGtttttttcctt                                                                                                             >  1:953889/1‑144 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTCTTAAAAAGCTTGGATTCTCGCAAAAGCAGGGTTTTTTCACtt                                                                                                             >  1:704062/1‑151 (MQ=255)
                                        tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTTAACAGCTTGGGTTCTCGCCAGAGCAGGCTTTTTTTACcttttg                                                                                                         >  2:1457987/1‑149 (MQ=255)
                                          gggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAGATTGGGTTTTTGCAAAAGGAGGGTTTTTTTACtttttggg                                                                                                       >  1:736257/1‑147 (MQ=255)
                                              tgctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAACTTGGGTTTTCGCAAAAGGAGGGGTTTTTCCCTTTTTGTGGttt                                                                                                   >  1:664380/1‑151 (MQ=255)
                                                 tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgc                                                                                                <  1:119914/151‑1 (MQ=255)
                                                  gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAGCTTGGGTTTTCGCAAAAAGAGGCCTTTTTCACTTTTTGTGCttttt                                                                                                 >  1:245374/1‑148 (MQ=255)
                                                   cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAACTTGGGTTTTCGCCAGAGGAGGGTTTTTTCTCCTTTTGGGCttttt                                                                                                 >  1:1203418/1‑147 (MQ=255)
                                                       aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAACTTGGGTTTTCGCCAGAGGAGGCTTTTTTTCCTTTTTGGGCgtttttttcg                                                                                            >  2:1454828/1‑143 (MQ=255)
                                                         ccTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTCTTTTAACAGCCTTGGTTCTCGCCAGAGCCGGGTTTTTTTTCTTTTTTTGCTTTTGCTGCGCTgg                                                                                        >  2:268227/1‑151 (MQ=255)
                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                      <  2:1529109/151‑1 (MQ=255)
                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                      <  1:476793/151‑1 (MQ=255)
                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                      <  1:298963/151‑1 (MQ=255)
                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                      <  1:827695/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg                                                                                 <  1:570102/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg                                                                                 <  1:1518357/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg                                                                                 <  2:842081/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg                                                                                 <  1:1103062/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg                                                                                 <  1:866064/151‑1 (MQ=255)
                                                                 atgtttCATGATGTCTGTTGTCCAGTAAGGGGAGCGATTCGGCTACTACGTCGTGCAGGGCGGACAGGC‑GGTTTTCTTCGTAAAACAGCTTGGATTCTCGCAAGGGCAGGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑c                                                                               <  2:833263/148‑1 (MQ=255)
                                                                 aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑c                                                                               <  1:1392364/151‑1 (MQ=255)
                                                                  tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑cc                                                                              <  1:1212519/151‑1 (MQ=255)
                                                                        aTGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CAttt                                                                       <  1:448891/151‑1 (MQ=255)
                                                                        aTGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CAttt                                                                       <  1:374722/151‑1 (MQ=255)
                                                                            tCTTTTTTGTCGAGTTATGGGTGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCAt                                                                   <  2:880957/151‑1 (MQ=255)
                                                                            tCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGGTTAACAACTTGGATTCTCTCCAGAGCCGGGCTTTTTCCCTTTTTGTGCTTTTTCCTCGCCGGTCTATTG‑CCT‑CAttttctt                                                                   >  2:1383804/1‑147 (MQ=255)
                                                                                       gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTCGTTTAAAAGCTTGGATTTTCGCCAAAGCAGGCTTTTTTTACTTTTTGTGCTTTTTTTTCGCTGGGCTATGG‑CCTCCTTTTCCTTTGGGGGCttt                                                         >  1:82242/1‑149 (MQ=255)
                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAGCTTGGATTTTCGCAAGAGCAGGGTTTTTTCACTTTTTGTGGTTTTTCTGCGCTGGGCTATGGCCCC‑CTTTTCCTTTgggggcattt                                                        >  1:88361/1‑147 (MQ=255)
                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTGTTTAAAAGCTTGGGTTTTCGCAAGAACAGGGTTTTTTTCCTTTTTGTGCTTTTTCTTCGCTGGGCTTTGG‑CCT‑CCTTTTCATTGGGGGGTATGt                                                       >  2:872271/1‑151 (MQ=255)
                                                                                                gAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAAAAGCCTGGGTTTTCGCCAAAGAAGGCTTTTTTCACTTTTTGTGGTTTTTCTGCGCTGGGCTATGG‑GCC‑CCTTTCCCTTtgcggccctttttgcggc                                                >  1:1063052/1‑138 (MQ=255)
                                                                                                gAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGac                                                <  2:740018/150‑1 (MQ=255)
                                                                                                   cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGgccac                                             <  2:170000/150‑1 (MQ=255)
                                                                                                      ttCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACctgc                                         <  2:333017/151‑1 (MQ=255)
                                                                                                        cGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACctgc‑tg                                      <  1:1095356/151‑1 (MQ=255)
                                                                                                             actacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGAc                                 <  2:1488927/151‑1 (MQ=255)
                                                                                                                     gTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAGCTTGGGTTTCCGCAAGAGCAGGGTTTTTTTACTTTTTGTGCTTTTTCTGCGCTGGGCTATGG‑CCT‑CATTTTCATTTGCGGGTCTTTTGGCGCCCCCCCCCTTGGGGCCCAAAcc                           >  1:1283524/1‑149 (MQ=255)
                                                                                                                     gTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTTTTACTTTTTGTGGTTTTGCTGTGCTGGGCTTTTG‑CCT‑CCTTTTCATTTGCGGCTTTTTCTGCGGCCCCCCGCTTGGGGGCCaaccc                           >  2:550593/1‑147 (MQ=255)
                                                                                                                      tGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCAc                        <  1:1418720/151‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTc                <  2:958425/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTc                <  2:308558/150‑1 (MQ=255)
                                                                                                                                     c‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAg          <  1:57357/150‑1 (MQ=255)
                                                                                                                                     c‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAg          <  2:664063/150‑1 (MQ=255)
                                                                                                                                       ggTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCa        >  1:880472/1‑151 (MQ=255)
                                                                                                                                        gTTTTCTTCGTTAACCAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCa        <  1:359941/150‑1 (MQ=255)
                                                                                                                                           ttcttcGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCACTtg    <  2:1495555/151‑1 (MQ=255)
                                                                                                                                             cttcGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCACTtgtg  <  1:207540/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                         
GGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGG‑CCT‑CATTTCCATTGGCGGCTATGTCGGCGACCACCTGC‑TGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTG  >  NC_000913/3640841‑3641131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑