Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG46.4% 160.8 / 51.4 114T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (12/48);  new base G (52/0);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 1.15e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675
                                                                                                                                                           |                                                                                                                                                      
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                                        <  4:1375785/150‑1 (MQ=255)
 aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtgg                                                                                                                                                      <  2:32099/151‑1 (MQ=255)
     cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTACGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGGATTTGGCGGCCGGGGGGgttt                                                                                                                                                  >  4:1252196/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCCGGGGGGGtttt                                                                                                                                                  >  4:21661/1‑151 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                                                                                                  <  2:878720/149‑1 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                                                                                                  <  4:512243/149‑1 (MQ=255)
        gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGgtttttt                                                                                                                                               >  4:1543547/1‑149 (MQ=255)
        gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                               <  3:218361/151‑1 (MQ=255)
        gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                               <  2:804569/151‑1 (MQ=255)
          gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                             <  1:400455/151‑1 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCAGGGGGGGttttttttt                                                                                                                                            >  3:476946/1‑150 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGGGGCCGGGGGGGtttttttttt                                                                                                                                            >  4:1569122/1‑151 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTTGGGGGCCGGGGGGGTtttttttt                                                                                                                                             >  2:101646/1‑149 (MQ=255)
               tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTGTTGGGGGGCCGGGGGGGtttttttttttcc                                                                                                                                         >  2:640931/1‑148 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGAATTGGGGGCCGGGGGGGGTTTTTTTTTcccc                                                                                                                                        >  4:145673/1‑151 (MQ=255)
                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCAGGGGGGGTTTTCTTTTTAc                                                                                                                                         >  3:301595/1‑148 (MQ=255)
                  tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGTTTTTTTTTTTACCCGt                                                                                                                                     >  2:863743/1‑151 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGTTTTGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                        >  1:286208/1‑147 (MQ=255)
                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                     <  4:1081819/150‑1 (MQ=255)
                      tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTTTGGCGCCCGGGGGGGTTTTCTTTTTAcc                                                                                                                                        >  1:254931/1‑144 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGCATTTGGGGCCCGGGGGGGTTTTTTTTTTAc                                                                                                                                         >  2:1264873/1‑142 (MQ=255)
                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTTTTTGGGGGCCGGGGGGGTTTTTTTTTTTACCAgggg                                                                                                                                   >  1:130021/1‑145 (MQ=255)
                         gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                              <  4:1210003/151‑1 (MQ=255)
                         gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACgg                                                                                                                              >  4:391190/1‑151 (MQ=255)
                          cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGATTGGGGCCCGGGGGGGtttttttttttccctggaggaggg                                                                                                                              >  2:1071267/1‑137 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAggg                                                                                                                             >  4:861338/1‑150 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGACCAACGCGCGGGGAGAGGCGGCTTGCGGATTTGGGGCCAGGGAGGGTTTTCTTTTTCACAGTGAGACggg                                                                                                                             >  4:1334321/1‑150 (MQ=255)
                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTTTTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGTGAGACggg                                                                                                                             >  4:1172587/1‑149 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                          <  1:114118/151‑1 (MQ=255)
                             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACggg                                                                                                                             >  3:264526/1‑148 (MQ=255)
                               tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTATTTGGCGGCAGGGGGGGTTTTTTTTTTACCAGCGAGACggg                                                                                                                             >  2:1642984/1‑146 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                      <  4:1616019/151‑1 (MQ=255)
                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGTGAGACGGGCaaaa                                                                                                                        >  4:1512313/1‑145 (MQ=255)
                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                                                                                                    <  4:637404/151‑1 (MQ=255)
                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                                                                                                    <  2:438808/151‑1 (MQ=255)
                                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGTTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGGGACAGCTg                                                                                                                    >  2:367363/1‑150 (MQ=255)
                                             ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGAAGGGGAACAACTGATTGgcc                                                                                                             >  2:452574/1‑148 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                                        <  3:9857/151‑1 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAATGAGAAGGGCAACAGCTCATTGCCCtt                                                                                                           >  4:120472/1‑146 (MQ=255)
                                                  cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGGGTTTTTTTTTTACCAGGGAGACGGGCAACAGCCGATTGCCCCTCCCCGc                                                                                                     >  3:450585/1‑151 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                   <  2:1506288/151‑1 (MQ=255)
                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCCGTg                                                                                                                                    >  4:256696/1‑117 (MQ=255)
                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                               <  2:256284/122‑1 (MQ=255)
                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                  <  4:285482/151‑1 (MQ=255)
                                                       tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTGTATTTGGCGGCACGGGGGGTTTTTTTTTCACCAAGGGGAGGGGAAAAAAATGATTGTCCTTCCCCGcc                                                                                                    >  4:1554446/1‑147 (MQ=255)
                                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCTCCAGGGTGGTTTTTCTTTTCGCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                             <  4:498190/151‑1 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTga                                                                                           <  2:1669926/151‑1 (MQ=255)
                                                              tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGGCAGGGGGGGTTTTCTTTTTACCCGTGAGAAGGGCAAAAGCTGATTGCCCTTCACCGCCTGGccc                                                                                              >  1:441215/1‑146 (MQ=255)
                                                                 ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGGCCGGGGGGTTTTTTTTTTCACCAGGGAGACGGGGAACAGCTGATTTCCCTTCACCGCCTGGCCCTgaggg                                                                                        >  3:41143/1‑147 (MQ=255)
                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                    <  2:1044497/151‑1 (MQ=255)
                                                                    gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTg                                                                                   <  2:933430/151‑1 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGGTTGCCCTTCACCCCCTGGCCCTGGGAGGGTTg                                                                                   >  1:287486/1‑149 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc                                                                                  <  4:809783/150‑1 (MQ=255)
                                                                                 tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGGGAGAAGGGCAACAACTGATTGCCCCTCACCCCCCGGCCCCGCGAGAGGTGCAGCAAGCGGGcc                                                                     >  2:943766/1‑151 (MQ=255)
                                                                                  cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGGGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGGTGCAGCAAGCGGGcc                                                                     >  2:1055521/1‑150 (MQ=255)
                                                                                  cGTGCCAGCTGCATTAATGAATCGGCAAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                                     <  2:453901/150‑1 (MQ=255)
                                                                                   gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCCGAGAGAGGTGCAGCAAGCGGGCCAc                                                                   >  4:761914/1‑151 (MQ=255)
                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCCGAGAGAGGTGCAGCAAACGGGCCACGCCGGgtt                                                           >  4:242067/1‑151 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAAGGAGAAGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTGGgtt                                                           >  1:140047/1‑150 (MQ=255)
                                                                                            gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                            <  2:279009/149‑1 (MQ=255)
                                                                                               ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGTTTTGGGGCCCGGGGGGGTTTTCTTTTTACCCGGGAGACGGGGCACAGCTGGTTTCCCCTCACCCCCCGGCCCTGGGAGAGTTGCAGCCAGCGGGCCCCCCCGGTTTGccc                                                       >  2:7352/1‑151 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGccca                                                      <  2:187521/150‑2 (MQ=255)
                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGATTTGGTGCTGTGGCGCAGGGCGGGTTTTTTTTTCAACAAAGAGACGGGCAACAGCCGATTGCCCTTCACCGCCCGGGCCTGAGAGAGCTGTAGCAAACGGCCCACGCCGGTTTTccccc                                                     >  3:67172/1‑150 (MQ=255)
                                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTCTTTTTACCCGTGGGACGGGGAAAAGCCGATTGGCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCCCCCCGGGTTTcccc                                                      >  2:1344626/1‑149 (MQ=255)
                                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                    <  4:159569/151‑1 (MQ=255)
                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTTTTTGGCGCCAGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCCACAGCTGATTGCCCCTCACCCCCTGGCCCCCGGGGAGTTGCCGCCAACGGTCCCCCCCGGTTTTccccc                                                     >  2:795965/1‑146 (MQ=255)
                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTCACCAATGGGAAGGGCAAAAGCTGATTGCCCCTCCCCGCCCGGCCCTGGGAGGGGTGCAGCAAACGGTCCACGCTGGTTTGGCCccgc                                                   >  2:105321/1‑149 (MQ=255)
                                                                                                       tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCGCCTGGCCCCGAGAGGGGTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCCgg                                                >  4:29304/1‑150 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                              <  2:467054/151‑1 (MQ=255)
                                                                                                          gCCAACGCGCGGGGAGAGGCGGATTTCGGTTTTGGCGCCGGGGGGGGTTTTCCTTTTTCCAGTTAGAAAGGCAACCACTGATTGCCCTTCACCCCCTGGCCCTGAGAGGGTTGTAGCAAGCGGGCCACGCTGGGTTGCCCCAGCAGGc                                               >  4:613900/1‑148 (MQ=255)
                                                                                                           ccAACGCGCGGGGAGAGGCGGGTTTGGTATTGGGCGCCCGGGTGGGTTTTTTTTTTTCCAGTGTGACCGGGAACAGCTGGTTTTCCTTCACCCCCCGGCCCTTAAAGAGATGCAGCAAAGCGGCCACCCTGGTTTTCCCCAGCAgg                                                >  4:408202/1‑146 (MQ=255)
                                                                                                             aaCGCGCGGGGAGAGGCGGTTTTGGTGTTTGGGGGCCGGGGGGGTTTTTTTTTTACCCGGGGGAAGGGGAACAGCTGATTGTCCCTCACCGCCTGGCCCTGAGAGGGGTGGAGCAAGCGGGCCCCGCTGGGTTTCCCCAGCAGGGGgaaa                                          >  2:1149069/1‑150 (MQ=255)
                                                                                                               cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                                       <  2:1250620/151‑1 (MQ=255)
                                                                                                                gcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                      <  4:314842/151‑1 (MQ=255)
                                                                                                                  gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                                     <  1:316766/150‑1 (MQ=255)
                                                                                                                   cgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                   <  4:249745/151‑1 (MQ=255)
                                                                                                                     gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                                 <  3:193705/151‑1 (MQ=255)
                                                                                                                      ggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                                <  2:1171108/151‑1 (MQ=255)
                                                                                                                               gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                        >  1:149049/1‑150 (MQ=255)
                                                                                                                                tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                       >  2:761963/1‑150 (MQ=255)
                                                                                                                                 ttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                      <  2:21622/150‑1 (MQ=255)
                                                                                                                                  tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  3:447593/150‑1 (MQ=255)
                                                                                                                                     agTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                  >  2:653072/2‑151 (MQ=255)
                                                                                                                                        tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                >  4:59748/1‑151 (MQ=255)
                                                                                                                                          aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                <  1:427243/150‑1 (MQ=255)
                                                                                                                                          aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg               <  4:453213/151‑1 (MQ=255)
                                                                                                                                               gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat             >  4:902642/1‑150 (MQ=255)
                                                                                                                                                  gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata          <  4:810677/151‑1 (MQ=255)
                                                                                                                                                   cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         <  2:145423/151‑1 (MQ=255)
                                                                                                                                                    gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  1:406289/1‑150 (MQ=255)
                                                                                                                                                     ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  2:1379862/1‑149 (MQ=255)
                                                                                                                                                     ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc        <  4:1072632/150‑1 (MQ=255)
                                                                                                                                                       aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg     <  1:114311/151‑1 (MQ=255)
                                                                                                                                                         ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa       >  2:299677/1‑145 (MQ=255)
                                                                                                                                                         ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt      >  4:320372/1‑148 (MQ=255)
                                                                                                                                                         ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg   <  2:1541634/151‑1 (MQ=255)
                                                                                                                                                          gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg   >  4:10669/1‑150 (MQ=255)
                                                                                                                                                          gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  <  3:140061/151‑1 (MQ=255)
                                                                                                                                                          gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  <  2:64911/151‑1 (MQ=255)
                                                                                                                                                           tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg   >  2:530441/1‑149 (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: