Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG41.6% 230.9 / 21.4 114T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (11/55);  new base G (47/0);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 1.44e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366377‑366674
                                                                                                                                                      |                                                                                                                                                     
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtg                                                                                                                                                      <  4:1006853/151‑1 (MQ=255)
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggt                                                                                                                                                       <  2:1029825/150‑1 (MQ=255)
 aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtgg                                                                                                                                                     <  1:114870/151‑1 (MQ=255)
 aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtgg                                                                                                                                                     <  2:559932/151‑1 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTC‑TTTGGCGCCAggggggg                                                                                                                                                    >  1:117725/1‑150 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtggt                                                                                                                                                    <  2:242266/151‑1 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCAGGGGGGtttttt                                                                                                                                               >  2:556617/1‑150 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGtttt                                                                                                                                                 <  4:1351522/149‑1 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGGGTA‑TTTGGGGCCCGGGGGGgttttg                                                                                                                                               >  2:1412521/1‑150 (MQ=255)
           tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGttttctttt                                                                                                                                           >  4:110511/1‑151 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGCGCCCGGGGGGGttttcttttt                                                                                                                                          >  3:138050/1‑150 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGGttttcttttta                                                                                                                                         >  3:237807/1‑149 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGttttctttttt                                                                                                                                         >  3:269374/1‑149 (MQ=255)
              cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCCGGGGGGTTTTTCTTTTCAc                                                                                                                                        >  2:464190/1‑151 (MQ=255)
               tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                                                                                                       <  3:355630/151‑1 (MQ=255)
                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTCAcc                                                                                                                                       >  4:165893/1‑149 (MQ=255)
                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAg                                                                                                                                     <  1:319754/151‑1 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTC‑TTGGGCGGCCGGGGGGTTTTTCTTTTCTCCa                                                                                                                                      >  3:442170/1‑148 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                    <  4:811354/150‑1 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                   <  2:759330/151‑1 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGTTTTCTTTTTACCAggg                                                                                                                                   >  4:295903/1‑147 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                               <  2:693999/151‑1 (MQ=255)
                              cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGTTTTCTTTTCACCAGGGAGACGGGCaa                                                                                                                         >  1:227328/1‑150 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGGAACAGc                                                                                                                     >  2:1197425/1‑151 (MQ=255)
                                  caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGAAGGGCAACAg                                                                                                                      >  1:221663/1‑149 (MQ=255)
                                  caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGATTTGGGCG‑CAGGGGGGTTTTTCCTTTCACCAGGGAGACGGGGAACAGc                                                                                                                     >  3:545612/1‑150 (MQ=255)
                                     tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTGGGGGCCAGGGGGGGTTTTTTTTTCACCAGGGAGACGGGGAACAGCTGAt                                                                                                                 >  1:560198/1‑151 (MQ=255)
                                       aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGCAACAGCTGGTTg                                                                                                               >  4:510560/1‑151 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                           <  4:209889/151‑1 (MQ=255)
                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAAAAGCTGATTGCCCTTc                                                                                                         >  2:614839/1‑150 (MQ=255)
                                                cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                                       <  3:284684/150‑1 (MQ=255)
                                                   tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGCGGCAGGGGGGGTTTTCTTTTCCCCAGTGAGAGGGGGAAAAACCGATTGCCCCTCAcccgc                                                                                                   >  4:865699/1‑148 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑Gcc                                                                                                  <  4:946211/150‑1 (MQ=255)
                                                      cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACAACTGATTGCCCTTTACC‑CCCttg                                                                                               >  4:946571/1‑149 (MQ=255)
                                                       tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGc                                                                                              <  3:26409/151‑1 (MQ=255)
                                                        gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTgg                                                                                               <  4:253016/149‑1 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCCGGGGGGTTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGccc                                                                                            >  2:7607/1‑151 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTga                                                                                         <  4:582091/151‑1 (MQ=255)
                                                             cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTgag                                                                                        >  3:11734/1‑151 (MQ=255)
                                                              tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACC‑CCCTGGCCCTgggg                                                                                       >  2:1337460/1‑148 (MQ=255)
                                                              tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGCGCCCGGGCGGTTTTTCTTTTTACCAGTGGGAGGGGCCACAGCTGATTGCCCCTTACC‑TCCTGGCCCTgcgg                                                                                       >  2:123819/1‑148 (MQ=255)
                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTgagag                                                                                      <  3:93199/151‑1 (MQ=255)
                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTgagag                                                                                      <  2:680878/150‑1 (MQ=255)
                                                                 ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTgagaga                                                                                     <  4:1242724/150‑1 (MQ=255)
                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGC‑TTTGGCGACAGGGGGGGTTTTCTTTTCACCAGTGAGGCGGGGCACCGCTGATTGCCCTTCCCC‑CCCCGGCCCTgcgtaaatt                                                                                  >  4:461193/1‑143 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑CCCCGGCCCTGAGAGAGGTgca                                                                               >  3:88787/1‑151 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTgca                                                                               <  4:460311/151‑1 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCCAGGGGGGTTTTTTTTTTTACCAGTGAGAAGGGGAACAGCTGAATGGCCTTCACCGGCCTGG‑CCTGAGAGAGGTggagca                                                                            >  1:21856/1‑151 (MQ=255)
                                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGAGAAGGGCAAAAGCTGATTGCCCTTCCCC‑CCCTGGCCCCGAGAGAGGTGCAGCAAGc                                                                         >  4:386078/1‑151 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑ATGGGGGCCAGGGGGGTTTTTCTTTTTCCCCCTGAGACGGGCAACAGCTGGTTTCCCTTTACC‑CCCTGGCCCCTCGGGAGGTgcagcacacg                                                                        >  4:1083850/1‑147 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                                       <  3:358998/121‑1 (MQ=255)
                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                       <  4:380301/149‑1 (MQ=255)
                                                                                   gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGCGAGAAGGGCAACAGCTGATTGGCCTTCACC‑CCCCGGCCCTGTGGGAGTTGCAGCAAGCGGTCCAc                                                                  >  4:1382628/1‑151 (MQ=255)
                                                                                      ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTA‑TTTGGCGCCACGGGGGGTTTTCTTTTCTCCCGTGGGACGGGCCACAGCTCATTTCCCCTCCCC‑GCCCGGCCCCGGGAGAGGTGCAGCAAACGGGCCACCCt                                                               >  4:642236/1‑151 (MQ=255)
                                                                                        aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                             <  2:1334057/151‑1 (MQ=255)
                                                                                        aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                             <  4:28488/151‑1 (MQ=255)
                                                                                          cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCAGGGGGGTTTTTCTTTTCACCAGGGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCCGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCCgggt                                                           >  4:213067/1‑149 (MQ=255)
                                                                                          cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                             <  1:486747/149‑1 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                         >  1:391580/1‑151 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTg                                                         >  1:125050/1‑151 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                          <  4:539148/150‑1 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtttt                                                         >  1:465255/1‑150 (MQ=255)
                                                                                              aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcc                                                       <  2:510083/151‑1 (MQ=255)
                                                                                                   tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGTTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGGGAGAGTTGCAGCAAGCGGGCCACGCCGGTTTGCCCcagc                                                  >  2:178227/1‑151 (MQ=255)
                                                                                                   tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGCT‑TTGGGCGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACCGGCCACAGCTGATTGCCCTTCACC‑CCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTTCACGCTGGGTTGCCCcagc                                                  >  4:79721/1‑151 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                              <  2:1355578/150‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAACGGGCCACGCTGGGTTGCCCCAGCAGGCg                                             >  2:166859/1‑150 (MQ=255)
                                                                                                          gCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaa                                           <  2:503919/151‑1 (MQ=255)
                                                                                                           ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                          <  2:123540/151‑1 (MQ=255)
                                                                                                           ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                          <  2:317288/151‑1 (MQ=255)
                                                                                                               cgcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                                      <  4:723195/151‑1 (MQ=255)
                                                                                                                gcgcgGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGGGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTGCCCCCGCAGGCGAAAATCCt                                     >  4:866898/1‑151 (MQ=255)
                                                                                                                gcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                     <  4:1352133/151‑1 (MQ=255)
                                                                                                                         gagGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                            <  2:510556/151‑1 (MQ=255)
                                                                                                                               gTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  2:1262976/151‑1 (MQ=255)
                                                                                                                                 ttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                    >  2:31008/1‑151 (MQ=255)
                                                                                                                                  tGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                   >  2:291783/1‑151 (MQ=255)
                                                                                                                                  tGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                   >  4:148714/1‑151 (MQ=255)
                                                                                                                                  tGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                    <  2:1121832/150‑1 (MQ=255)
                                                                                                                                    cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                 <  4:171723/151‑1 (MQ=255)
                                                                                                                                       a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               <  4:1057497/150‑1 (MQ=255)
                                                                                                                                       a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  2:1454171/1‑151 (MQ=255)
                                                                                                                                          tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa             <  4:1401558/150‑1 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          <  4:1166521/151‑1 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          >  4:1245478/1‑151 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          <  2:1078172/151‑1 (MQ=255)
                                                                                                                                             gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata         <  2:110510/151‑1 (MQ=255)
                                                                                                                                              cgcCAGGGTGGTTTTCCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa        <  2:393987/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc       <  1:11734/150‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      <  4:123821/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      <  4:545610/151‑1 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     <  2:212220/151‑1 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     <  2:165891/151‑1 (MQ=255)
                                                                                                                                                  aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  2:587981/151‑1 (MQ=255)
                                                                                                                                                  aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  2:344241/151‑1 (MQ=255)
                                                                                                                                                   gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    >  2:1337049/1‑150 (MQ=255)
                                                                                                                                                    ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     <  3:391575/148‑1 (MQ=255)
                                                                                                                                                    ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     <  3:311321/148‑1 (MQ=255)
                                                                                                                                                    ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAAATg    >  4:110712/1‑149 (MQ=255)
                                                                                                                                                     gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAAATGAg  >  3:148120/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACC‑GCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366377‑366674

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