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

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675
                                                                                                                                                       |                                                                                                                                                     
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                                      <  1:32474/151‑1 (MQ=255)
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                                      <  1:48863/151‑1 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGGATTTGGCGCCCgggggg                                                                                                                                                     >  1:717126/1‑150 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTATTGGGCGGCCgggggg                                                                                                                                                     >  2:866949/1‑150 (MQ=255)
    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGGGTGTTTGGGGCCCggggggggt                                                                                                                                                  >  2:728750/1‑148 (MQ=255)
     cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGgttt                                                                                                                                                 >  2:661455/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCAGGGGGGgtttt                                                                                                                                                >  2:785957/1‑151 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                                                                                                 <  2:269689/149‑1 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                                                                                                          <  2:563811/151‑1 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTAc                                                                                                                                        >  1:395108/1‑146 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                    <  1:670585/150‑1 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTTTTGGGGGCCCGGGGGGGtttttttttttcccgtg                                                                                                                                   >  2:451937/1‑141 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                <  2:197542/150‑1 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                <  1:115505/150‑1 (MQ=255)
                         gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                             <  2:984952/151‑1 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTCTTTTCACCCGGGAGAAgggg                                                                                                                           >  2:701964/1‑150 (MQ=255)
                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTCTCCCGTGAGAAggg                                                                                                                            >  1:13397/1‑149 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGGGTTTTTTTTTTACCAGGGAGAAggg                                                                                                                            >  1:671568/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                        <  2:262288/149‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                      <  1:237897/151‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGCGGCGGTTTGCGGATTTGGGGGCCGGGGGGTTTTTCTTTTTCCCAGCGAGGCGGGCCACAGCTggt                                                                                                                 >  2:905719/1‑148 (MQ=255)
                                       aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtt                                                                                                                >  1:651338/1‑150 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGCATTTGGAGCCCGGGGGGGTTTTTTTTTCACCCGTGAGAGGGGCAACAGCCGAtag                                                                                                               >  2:32108/1‑145 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                            >  1:358991/1‑150 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGGCTTTGGGGGCCGGGGGGGTTTTCTTTTTCCCCGGGGGAAGGGGAACAGCTGGTTGcccct                                                                                                          >  2:48024/1‑146 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGCTTGGGGGCCCGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGACTGGCCCTCAcc                                                                                                      >  2:120724/1‑150 (MQ=255)
                                                  cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTGTAGAGGGCCTGGGGGGGTTTTCTTTTTTCCCAAGAGATGGGCAACAGCttgggttccttccgc                                                                                                      >  1:756185/1‑135 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  1:116134/151‑1 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  1:786151/151‑1 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCCGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGccc                                                                                             >  2:978510/1‑151 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAACTGATTTTCCTTCACCCCCTGGccc                                                                                             >  1:410209/1‑148 (MQ=255)
                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                         <  1:98807/149‑1 (MQ=255)
                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCCGATTTCCCCTCACCCCCTGGccccg                                                                                           >  1:638666/1‑144 (MQ=255)
                                                                 ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                                     <  1:145049/151‑1 (MQ=255)
                                                                  cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                                    <  1:636222/151‑1 (MQ=255)
                                                                    gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTACGAATTGGGGGCCCGGGGGGTTTTTTTTTTCACCACGGAGACGGGCAACAGGCGATTGCCCCTCACCGCCTCGCCCCcagaggcgt                                                                                   >  1:259008/1‑146 (MQ=255)
                                                                     tCGGGAAACCTGTCGTGCCAGATGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTCTCGGATAGGGGGCCCGGGGGGTGTTTTTTTTCCCCAACGAGACGGGCCACACCTGCTTGCCCCTCTCCCCCCCGCCCttcgag                                                                                       >  2:697663/1‑140 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                <  1:176545/151‑1 (MQ=255)
                                                                        ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCACTGGGACGGGCAAAACCTGATTTCCCCTCACCCCCCGGGCCTgcgaggggttca                                                                                >  1:952228/1‑142 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  1:456897/151‑1 (MQ=255)
                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTCTTGGGGGCCCGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCCGAGAGAGTTgga                                                                                >  2:974106/1‑145 (MQ=255)
                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGTTTTTTTTTTCACCAGTGAGACGGGCAACAACTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                            >  1:788399/1‑150 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                         <  1:452055/151‑1 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGGCCGGGGGGGTTTTTTTTTTTCCAATGGGACGGGGAACAGCTGATTTCCCCTTACCCCCCGGCCCTGGGAGAGATGCGGCAAACgg                                                                        >  1:420720/1‑150 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTGGGCGGCCGGGGGGGTTTTTTTTTTCCCAGGGAGACGGGGCAAAGCTGAGTGGCCTCCCCCCCCCGGCCCCGGGAGGGTTGCCGGAAGCgg                                                                        >  2:370408/1‑150 (MQ=255)
                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                       <  2:424842/150‑1 (MQ=255)
                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTgg                                                              >  1:330890/1‑150 (MQ=255)
                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGGCGGGGCACAGCTGGTTGCCCTTCACCGCCTGGCCCTGGGAGAGGTGCAGCAAGCGGTCCACGCTgg                                                              >  1:256121/1‑150 (MQ=255)
                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCAccccgggtt                                                           >  1:766314/1‑143 (MQ=255)
                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAAAAGCTGATTGGCCCTCACCGCCTGGCCCCGAGAGGGTTGGAGGAAGCGGGCCCCGCTGGgtt                                                           >  1:721622/1‑149 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCCTCCCCCCCTGGCCCTGCGAGAGGTGCAGCAAGCGGTCCCCCCTGGTTTTcccc                                                      >  1:496577/1‑150 (MQ=255)
                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                 <  1:770709/151‑1 (MQ=255)
                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                                <  1:228792/151‑1 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                               <  1:348380/150‑1 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                              <  2:992893/151‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                             <  1:205273/151‑1 (MQ=255)
                                                                                                                  gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                                    <  2:824678/151‑1 (MQ=255)
                                                                                                                                  tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  1:774128/150‑1 (MQ=255)
                                                                                                                                   gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGCGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                   >  1:963860/1‑151 (MQ=255)
                                                                                                                                    cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                  <  2:293459/151‑1 (MQ=255)
                                                                                                                                     gTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCTCGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                 <  2:75928/151‑1 (MQ=255)
                                                                                                                                          tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg               >  1:604010/1‑149 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  1:564740/151‑1 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata            >  1:58048/1‑150 (MQ=255)
                                                                                                                                              cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata          >  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)
                                                                                                                                                       |                                                                                                                                                     
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  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: