Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG43.8% 104.2 / 24.9 64T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (7/29);  new base G (28/0);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 7.97e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675
                                                                                                                                                       |                                                                                                                                                     
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGggtg                                                                                                                                                      <  1:32474/151‑1 (MQ=255)
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGggtg                                                                                                                                                      <  1:48863/151‑1 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGGA‑TTTGGC‑GCCCgggggg                                                                                                                                                     >  1:717126/1‑150 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTA‑TTGGGC‑GGCCgggggg                                                                                                                                                     >  2:866949/1‑150 (MQ=255)
    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGGGTG‑TTTGGG‑GCCCggggggggt                                                                                                                                                  >  2:728750/1‑148 (MQ=255)
     cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGC‑GCCAGGGTGGgttt                                                                                                                                                 >  2:661455/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GGCAGGGGGGgtttt                                                                                                                                                >  2:785957/1‑151 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGtttt                                                                                                                                                 <  2:269689/149‑1 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTttttcttttc                                                                                                                                          <  2:563811/151‑1 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTTAc                                                                                                                                        >  1:395108/1‑146 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                    <  1:670585/150‑1 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTT‑TTGGGG‑GCCCGGGGGGGtttttttttttcccgtg                                                                                                                                   >  2:451937/1‑141 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                <  2:197542/150‑1 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                <  1:115505/150‑1 (MQ=255)
                         gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                             <  2:984952/151‑1 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCAGGGGGGGTTTTCTTTTCACCCGGGAGAAgggg                                                                                                                           >  2:701964/1‑150 (MQ=255)
                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCCGGGGGGGTTTTTTTTTCTCCCGTGAGAAggg                                                                                                                            >  1:13397/1‑149 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGG‑GCCAGGGGGGGTTTTTTTTTTACCAGGGAGAAggg                                                                                                                            >  1:671568/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                        <  2:262288/149‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                      <  1:237897/151‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGCGGCGGTTTGCGGA‑TTTGGG‑GGCCGGGGGGTTTTTCTTTTTCCCAGCGAGGCGGGCCACAGCTggt                                                                                                                 >  2:905719/1‑148 (MQ=255)
                                       aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                >  1:651338/1‑150 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGCA‑TTTGGA‑GCCCGGGGGGGTTTTTTTTTCACCCGTGAGAGGGGCAACAGCCGAtag                                                                                                               >  2:32108/1‑145 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                            >  1:358991/1‑150 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGGC‑TTTGGG‑GGCCGGGGGGGTTTTCTTTTTCCCCGGGGGAAGGGGAACAGCTGGTTGcccct                                                                                                          >  2:48024/1‑146 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGC‑TTGGGG‑GCCCGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGACTGGCCCTCAcc                                                                                                      >  2:120724/1‑150 (MQ=255)
                                                  cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTG‑TAGAGG‑GCCTGGGGGGGTTTTCTTTTTTCCCAAGAGATGGGCAACAGCttgggttccttccgc                                                                                                      >  1:756185/1‑135 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  1:116134/151‑1 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  1:786151/151‑1 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTTCCCGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGccc                                                                                             >  2:978510/1‑151 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGA‑TTTGGG‑GCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAACTGATTTTCCTTCACCCCCTGGccc                                                                                             >  1:410209/1‑148 (MQ=255)
                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                         <  1:98807/149‑1 (MQ=255)
                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCCGATTTCCCCTCACCCCCTGGccccg                                                                                           >  1:638666/1‑144 (MQ=255)
                                                                 ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                                     <  1:145049/151‑1 (MQ=255)
                                                                  cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                                    <  1:636222/151‑1 (MQ=255)
                                                                    gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTACGAA‑TTGGGG‑GCCCGGGGGGTTTTTTTTTTCACCACGGAGACGGGCAACAGGCGATTGCCCCTCACCGCCTCGCCCCcagaggcgt                                                                                   >  1:259008/1‑146 (MQ=255)
                                                                     tCGGGAAACCTGTCGTGCCAGATGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTCTCGGA‑TAGGGG‑GCCCGGGGGGTGTTTTTTTTCCCCAACGAGACGGGCCACACCTGCTTGCCCCTCTCCCCCCCGCCCttcgag                                                                                       >  2:697663/1‑140 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                <  1:176545/151‑1 (MQ=255)
                                                                        ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTTCCACTGGGACGGGCAAAACCTGATTTCCCCTCACCCCCCGGGCCTgcgaggggttca                                                                                >  1:952228/1‑142 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  1:456897/151‑1 (MQ=255)
                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTC‑TTGGGG‑GCCCGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCCGAGAGAGTTgga                                                                                >  2:974106/1‑145 (MQ=255)
                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCCGGGGGGTTTTTTTTTTCACCAGTGAGACGGGCAACAACTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                            >  1:788399/1‑150 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                         <  1:452055/151‑1 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGG‑GCCGGGGGGGTTTTT‑TTTTTTCCAATGGGACGGGGAACAGCTGATTTCCCCTTACCCCCCGGCCCTGGGAGAGATGCGGCAAACgg                                                                        >  1:420720/1‑150 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTGGGCGGCCGGGGGGGTTTTTTTTTTC‑CCAGGGAGACGGGGCAAAGCTGAGTGGCCTCCCCCCCCCGGCCCCGGGAGGGTTGCCGGAAGCgg                                                                        >  2:370408/1‑150 (MQ=255)
                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                       <  2:424842/150‑1 (MQ=255)
                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTgg                                                              >  1:330890/1‑150 (MQ=255)
                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGC‑GCCAGGGGGGGTTTTCTTTTCACCAGTGAGGCGGGGCACAGCTGGTTGCCCTTCACCGCCTGGCCCTGGGAGAGGTGCAGCAAGCGGTCCACGCTgg                                                              >  1:256121/1‑150 (MQ=255)
                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCAccccgggtt                                                           >  1:766314/1‑143 (MQ=255)
                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGA‑TTGGGC‑GCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAAAAGCTGATTGGCCCTCACCGCCTGGCCCCGAGAGGGTTGGAGGAAGCGGGCCCCGCTGGgtt                                                           >  1:721622/1‑149 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCCTCCCCCCCTGGCCCTGCGAGAGGTGCAGCAAGCGGTCCCCCCTGGTTTTcccc                                                      >  1:496577/1‑150 (MQ=255)
                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                 <  1:770709/151‑1 (MQ=255)
                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                                <  1:228792/151‑1 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                               <  1:348380/150‑1 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                              <  2:992893/151‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                             <  1:205273/151‑1 (MQ=255)
                                                                                                                  gcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                                    <  2:824678/151‑1 (MQ=255)
                                                                                                                                  tGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  1:774128/150‑1 (MQ=255)
                                                                                                                                   gCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGCGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                   >  1:963860/1‑151 (MQ=255)
                                                                                                                                    cGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                  <  2:293459/151‑1 (MQ=255)
                                                                                                                                     gTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCTCGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                 <  2:75928/151‑1 (MQ=255)
                                                                                                                                          tGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg               >  1:604010/1‑149 (MQ=255)
                                                                                                                                            ggC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  1:564740/151‑1 (MQ=255)
                                                                                                                                            ggC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata            >  1:58048/1‑150 (MQ=255)
                                                                                                                                              c‑gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata          >  2:103752/1‑150 (MQ=255)
                                                                                                                                                gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  1:44606/1‑150 (MQ=255)
                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg     >  1:663810/1‑150 (MQ=255)
                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa    <  2:30441/151‑1 (MQ=255)
                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa    <  1:99686/151‑1 (MQ=255)
                                                                                                                                                      gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  >  1:598805/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675

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