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

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670
                                                                                                                                                      |                                                                                                                                                 
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                                  <  2:1089872/151‑1 (MQ=255)
  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtggt                                                                                                                                                <  2:1246047/151‑1 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTATTTGGCGCCCGGGTGGgtttt                                                                                                                                            >  1:612471/1‑151 (MQ=255)
           tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTTGGCGCCAGGGGGGGttttttttt                                                                                                                                       >  2:576094/1‑151 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGGCCGGGGGGGTtttttttt                                                                                                                                        >  2:1024442/1‑149 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                  <  2:73167/151‑1 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                  <  2:825128/151‑1 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCCGGGGGGGTTTTCTTTTTTCCAGTg                                                                                                                               >  2:630983/1‑149 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGGCAGGGGGGGTTTTTTTTTCACCAGTGAGACggg                                                                                                                        >  1:404603/1‑150 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGGGAGACggg                                                                                                                        >  1:95691/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCAGGGGGGGTTTTTTTTTTCCCCGTGAGACGGGGaa                                                                                                                     >  2:687769/1‑148 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                    <  2:1142593/149‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                   <  1:1016093/150‑1 (MQ=255)
                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                  <  1:639739/151‑1 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                 <  2:689241/151‑1 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                 <  1:1225264/151‑1 (MQ=255)
                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACa                                                                                                                   >  1:159250/1‑147 (MQ=255)
                                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGTTTGGCGCCCGGGGGGTTTTTTTTTTTCCCAGTGAGACGGGGaa                                                                                                                     >  2:961618/1‑144 (MQ=255)
                                     tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                              <  1:641255/150‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                             <  1:1133466/150‑1 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAtttgc                                                                                                         >  1:1421227/1‑148 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGGTTg                                                                                                           >  2:3718/1‑149 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                          <  2:1334053/150‑1 (MQ=255)
                                             ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTTGGGGGCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGGAACAGCTGATTGcccct                                                                                                      >  2:327147/1‑148 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGCGGGACGGGGGAAAGCTGATTGcccct                                                                                                      >  2:833754/1‑146 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGCATTTGGCGGCCGGGGGGTTTTTTTTTTTTCCCGAGAGAAGGGGAAAAGCTGGtt                                                                                                            >  2:1277446/1‑140 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGAAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTTGGCGGCGGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                                  >  2:397599/1‑150 (MQ=255)
                                                   tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTCTTTTTTCCACAGAGACGGGCCAAAGGTGATTGCCCTTCACCGc                                                                                                >  2:156630/1‑150 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGccc                                                                                         <  1:714846/151‑1 (MQ=255)
                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGGTTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGTGGCGGGGAACAGCTGATTTGCCCTCTCCGCCCGGgcc                                                                                         >  2:1005401/1‑151 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCCGGGGGGGTTTTTTTTTTACCAGGGGGACGGGGAACAACCGGTTGGCCGTCCCCCCCCGGGCCTgg                                                                                      >  2:1103138/1‑150 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGGAATGGGCGTCACGGTGGTTTTTTTTTTTTCCAGGGAGAGGGGGGACAACTGATTTGCCCTCCCCACCCGGccc                                                                                         >  2:234233/1‑148 (MQ=255)
                                                              tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCCGGCCCt                                                                                        >  1:1308145/1‑147 (MQ=255)
                                                                     tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCGCCTGGCCCTGagagaggt                                                                               >  1:261846/1‑147 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                            <  2:936579/151‑1 (MQ=255)
                                                                       gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTTCCAAGGAGACGGGGAACAGCTGATTTGCCTTCACCCCCCCGCCCTgggaggggt                                                                               >  1:331349/1‑143 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCCGGGGGGTTTTTTTTTTTACCCGTGAGACGGGCAAAAGCTGATTGCCCTTCACCGCCCGGCCCCgagag                                                                                   >  1:109463/1‑141 (MQ=255)
                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGGTTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAAAAACTGATTTGCCTTCACCCCCCGGCCCCgcgaggggt                                                                               >  2:1273346/1‑141 (MQ=255)
                                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                      <  1:1103161/151‑1 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGGCCTGGGGGGGGTGCAGCaa                                                                        >  1:396440/1‑148 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                     <  2:461214/151‑1 (MQ=255)
                                                                             ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                     <  2:93683/151‑1 (MQ=255)
                                                                              cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                        >  2:980109/1‑147 (MQ=255)
                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGAGTGCGTAATTGGGGCCGCGGGGGGTTTTTTTTTCTCCAGCGAGAAGGGCAACAGCTTATTGCCCTTTTTCCCCTGGCCCTGAGCGAGTTGCCGGAAgggg                                                                    >  1:1016564/1‑150 (MQ=255)
                                                                                   gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTGGGGGCCCGGGGGGGTTTTCTTTTTACCAGAGAGACGGGCAACAGCTGATTGCCCTTTACCGCCTGGCCCTGGGAGAGTTGCAGCAAGCGGTcccc                                                               >  1:115527/1‑149 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACg                                                              <  1:158094/150‑1 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACg                                                              <  1:282159/150‑1 (MQ=255)
                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                             <  2:939021/151‑1 (MQ=255)
                                                                                      ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                            <  1:615143/151‑1 (MQ=255)
                                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGTTTTTTTTTTTTCCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTTCCGCAAGCGGGCCAc                                                               >  1:925388/1‑147 (MQ=255)
                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCAACAGCTGATTGGCCCTCACCCCCCGGGCCCGAGAGAGGTGCAGCAAACGGGCCACGCTgg                                                         >  2:1284469/1‑148 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                       <  1:744949/149‑1 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                     <  1:1180344/151‑1 (MQ=255)
                                                                                               ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGGGACGGGGAAAAACTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTGCAGCAAGCGGGCCCCGCTGGTTTGccc                                                  >  1:929006/1‑151 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGGTGCAGCAAGCGGGCCACGCCGGGTTGgccc                                                 >  2:311511/1‑150 (MQ=255)
                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                           <  1:3718/151‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTGGGTTTTccccgccgg                                            >  1:1104995/1‑143 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCCCCCCGGGTTTCCCCACCAGGCg                                         >  1:848209/1‑150 (MQ=255)
                                                                                                             accgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                    <  2:545776/149‑1 (MQ=255)
                                                                                                                       gagaGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                          <  1:916881/151‑1 (MQ=255)
                                                                                                                        agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                         <  2:645026/151‑1 (MQ=255)
                                                                                                                          agGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                       <  1:576107/151‑1 (MQ=255)
                                                                                                                           ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                       <  2:57750/150‑1 (MQ=255)
                                                                                                                            gCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                     <  1:238913/151‑1 (MQ=255)
                                                                                                                             cGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                     >  2:1305267/1‑150 (MQ=255)
                                                                                                                              ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                   <  1:665448/151‑1 (MQ=255)
                                                                                                                                   gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               <  1:1024463/150‑1 (MQ=255)
                                                                                                                                   gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               <  1:687788/150‑1 (MQ=255)
                                                                                                                                    cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc             <  2:181274/151‑1 (MQ=255)
                                                                                                                                      tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg           >  2:555973/1‑151 (MQ=255)
                                                                                                                                          tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa         <  1:321778/150‑1 (MQ=255)
                                                                                                                                            ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata       >  2:1330151/1‑150 (MQ=255)
                                                                                                                                             gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGCCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata     <  2:402377/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  <  2:928985/151‑1 (MQ=255)
                                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  <  2:205619/151‑1 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  2:563367/1‑150 (MQ=255)
                                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGACATAACa  >  2:829619/1‑150 (MQ=255)
                                                                                                                                                  aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  2:593143/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  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: