Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG42.1% 140.1 / 22.1 76T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (8/36);  new base G (32/0);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 3.41e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670
                                                                                                                                                      |                                                                                                                                                 
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtg                                                                                                                                                  <  2:1089872/151‑1 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtggt                                                                                                                                                <  2:1246047/151‑1 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTA‑TTTGGCGCCCGGGTGGgtttt                                                                                                                                            >  1:612471/1‑151 (MQ=255)
           tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGA‑TTTGGCGCCAGGGGGGGttttttttt                                                                                                                                       >  2:576094/1‑151 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGGCCGGGGGGGTtttttttt                                                                                                                                        >  2:1024442/1‑149 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                  <  2:73167/151‑1 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                  <  2:825128/151‑1 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTCTTTTTTCCAGTg                                                                                                                               >  2:630983/1‑149 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGGCAGGGGGGGTTTTTTTTTCACCAGTGAGACggg                                                                                                                        >  1:404603/1‑150 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGGGAGACggg                                                                                                                        >  1:95691/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCAGGGGGGGTTTTTTTTTTCCCCGTGAGACGGGGaa                                                                                                                     >  2:687769/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                    <  2:1142593/149‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                   <  1:1016093/150‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                  <  1:639739/151‑1 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                 <  2:689241/151‑1 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                 <  1:1225264/151‑1 (MQ=255)
                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACa                                                                                                                   >  1:159250/1‑147 (MQ=255)
                                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTG‑TTTGGCGCCCGGGGGGTTTTTTTTTTTCCCAGTGAGACGGGGaa                                                                                                                     >  2:961618/1‑144 (MQ=255)
                                     tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                              <  1:641255/150‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                             <  1:1133466/150‑1 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtttgc                                                                                                         >  1:1421227/1‑148 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGGTTg                                                                                                           >  2:3718/1‑149 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                          <  2:1334053/150‑1 (MQ=255)
                                             ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCGGGAGACGGGGAACAGCTGATTGcccct                                                                                                      >  2:327147/1‑148 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGCGGGACGGGGGAAAGCTGATTGcccct                                                                                                      >  2:833754/1‑146 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGCA‑TTTGGCGGCCGGGGGGTTTTTTTTTTTTCCCGAGAGAAGGGGAAAAGCTGGtt                                                                                                            >  2:1277446/1‑140 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGAAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTA‑TTTGGCGGCGGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                                  >  2:397599/1‑150 (MQ=255)
                                                   tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTTCCACAGAGACGGGCCAAAGGTGATTGCCCTTCACCGc                                                                                                >  2:156630/1‑150 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGccc                                                                                         <  1:714846/151‑1 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGTGGCGGGGAACAGCTGATTTGCCCTCTCCGCCCGGgcc                                                                                         >  2:1005401/1‑151 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCAGGGGGACGGGGAACAACCGGTTGGCCGTCCCCCCCCGGGCCTgg                                                                                      >  2:1103138/1‑150 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGGA‑ATGGGCGTCACGGTGGTTTTTTTTTTTTCCAGGGAGAGGGGGGACAACTGATTTGCCCTCCCCACCCGGccc                                                                                         >  2:234233/1‑148 (MQ=255)
                                                              tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCCGGCCCt                                                                                        >  1:1308145/1‑147 (MQ=255)
                                                                     tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCGCCTGGCCCTGagagaggt                                                                               >  1:261846/1‑147 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                            <  2:936579/151‑1 (MQ=255)
                                                                       gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTTCCAAGGAGACGGGGAACAGCTGATTTGCCTTCACCCCCCCGCCCTgggaggggt                                                                               >  1:331349/1‑143 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCCGGGGGGTTTTTTTTTTTACCCGTGAGACGGGCAAAAGCTGATTGCCCTTCACCGCCCGGCCCCgagag                                                                                   >  1:109463/1‑141 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGG‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAAAAACTGATTTGCCTTCACCCCCCGGCCCCgcgaggggt                                                                               >  2:1273346/1‑141 (MQ=255)
                                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                      <  1:1103161/151‑1 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGGCCTGGGGGGGGTGCAGCaa                                                                        >  1:396440/1‑148 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                     <  2:461214/151‑1 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                     <  2:93683/151‑1 (MQ=255)
                                                                              cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                        >  2:980109/1‑147 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGAGTGCGTA‑ATTGGGGCCGCGGGGGGTTTTTTTTTCTCCAGCGAGAAGGGCAACAGCTTATTGCCCTTTTTCCCCTGGCCCTGAGCGAGTTGCCGGAAgggg                                                                    >  1:1016564/1‑150 (MQ=255)
                                                                                   gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTGGGGGCCCGGGGGGGTTTTCTTTTTACCAGAGAGACGGGCAACAGCTGATTGCCCTTTACCGCCTGGCCCTGGGAGAGTTGCAGCAAGCGGTcccc                                                               >  1:115527/1‑149 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACg                                                              <  1:158094/150‑1 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACg                                                              <  1:282159/150‑1 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                             <  2:939021/151‑1 (MQ=255)
                                                                                      ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                            <  1:615143/151‑1 (MQ=255)
                                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGTTTTTTTTTTTTCCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTTCCGCAAGCGGGCCAc                                                               >  1:925388/1‑147 (MQ=255)
                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCAACAGCTGATTGGCCCTCACCCCCCGGGCCCGAGAGAGGTGCAGCAAACGGGCCACGCT‑gg                                                         >  2:1284469/1‑148 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGtt                                                       <  1:744949/149‑1 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTg                                                     <  1:1180344/151‑1 (MQ=255)
                                                                                               ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGGGACGGGGAAAAACTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTGCAGCAAGCGGGCCCCGCT‑GGTTTGccc                                                  >  1:929006/1‑151 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGGGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGGTGCAGCAAGCGGGCCACGCC‑GGGTTGgccc                                                 >  2:311511/1‑150 (MQ=255)
                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAgg                                           <  1:3718/151‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTGGGTTTTccccgccgg                                            >  1:1104995/1‑143 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCCCCCC‑GGGTTTCCCCACCAGGCg                                         >  1:848209/1‑150 (MQ=255)
                                                                                                             accgcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAAt                                    <  2:545776/149‑1 (MQ=255)
                                                                                                                       gagaGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                          <  1:916881/151‑1 (MQ=255)
                                                                                                                        agagGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                         <  2:645026/151‑1 (MQ=255)
                                                                                                                          agGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                       <  1:576107/151‑1 (MQ=255)
                                                                                                                           ggCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                       <  2:57750/150‑1 (MQ=255)
                                                                                                                            gCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                     <  1:238913/151‑1 (MQ=255)
                                                                                                                             cGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                     >  2:1305267/1‑150 (MQ=255)
                                                                                                                              ggTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                   <  1:665448/151‑1 (MQ=255)
                                                                                                                                   gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               <  1:1024463/150‑1 (MQ=255)
                                                                                                                                   gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               <  1:687788/150‑1 (MQ=255)
                                                                                                                                    cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc             <  2:181274/151‑1 (MQ=255)
                                                                                                                                      tA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg           >  2:555973/1‑151 (MQ=255)
                                                                                                                                          tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa         <  1:321778/150‑1 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata       >  2:1330151/1‑150 (MQ=255)
                                                                                                                                             gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGCCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata     <  2:402377/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  <  2:928985/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  <  2:205619/151‑1 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  2:563367/1‑150 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGACATAACa  >  2:829619/1‑150 (MQ=255)
                                                                                                                                                  aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  2:593143/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670

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