Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG45.3% 122.7 / 29.8 75T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (4/37);  new base G (34/0);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 3.25e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676
                                                                                                                                                     |                                                                                                                                                       
aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTTGGGGGCAgggggg                                                                                                                                                       >  1:712659/1‑151 (MQ=255)
  gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                                                                                                     <  2:718439/151‑1 (MQ=255)
  gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAACTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTTTGGCGCCCgggggg                                                                                                                                                       >  1:1118343/1‑149 (MQ=255)
        ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGATTGGGGGCCTGGGGGGgtttttt                                                                                                                                                >  1:541876/1‑148 (MQ=255)
          tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCAGGGGGGGttttttttt                                                                                                                                             >  1:459575/1‑151 (MQ=255)
           gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGttttctttt                                                                                                                                             >  2:1239685/1‑150 (MQ=255)
           gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                                                                                                            <  2:606327/151‑1 (MQ=255)
              tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGAATTGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                         >  1:1155956/1‑151 (MQ=255)
               aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                                          <  2:226635/149‑1 (MQ=255)
                aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTGGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                         >  2:880202/1‑149 (MQ=255)
                   agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                    <  1:1265044/151‑1 (MQ=255)
                   agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                    <  2:1110740/151‑1 (MQ=255)
                     tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGGTTTGGGGCCCGGGGGGGTTTTTTTTTTccccgtggg                                                                                                                                    >  1:41322/1‑145 (MQ=255)
                          tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTCCCCGTGGGACggg                                                                                                                              >  2:454477/1‑150 (MQ=255)
                              tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGGCggg                                                                                                                              >  2:216689/1‑146 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAAAAGGGGaa                                                                                                                           >  1:820936/1‑148 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                          <  2:238339/149‑1 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                         <  1:634149/150‑1 (MQ=255)
                                 caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                       <  2:1000488/150‑1 (MQ=255)
                                 caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                                      <  2:1136908/151‑1 (MQ=255)
                                   ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGAAGGGGaa                                                                                                                           >  1:218485/1‑144 (MQ=255)
                                    tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                                    <  1:338058/150‑1 (MQ=255)
                                        tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                               <  1:889595/150‑1 (MQ=255)
                                             tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTTGGCGGCCGGGGGGGTTTTCTTTTTTCCAGTGGGGCGGGGAAAAGCGGAATGcccgt                                                                                                           >  2:535521/1‑147 (MQ=255)
                                             tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTc                                                                                                          <  1:430527/150‑1 (MQ=255)
                                             tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTc                                                                                                          <  2:719129/150‑1 (MQ=255)
                                               cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGACTCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGGTTGGGCGCCCGGGGGGGTTTTCTTTTTACCCGTGAGGAGGGCAACAGCTGATTTCCCTTcccc                                                                                                       >  1:132723/1‑151 (MQ=255)
                                                gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCGACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGGCCGGGGGGTTTTTTTTTTTACCCGAGAGAAGGGCAACAGCAGGTTTCCCtt                                                                                                           >  2:905463/1‑146 (MQ=255)
                                                gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGGTTTGGGGCCGGGGGGGGTTTTTTTTTTAACAGAGGGACGGGGAACAGATAAtttgcc                                                                                                              >  1:1260707/1‑140 (MQ=255)
                                                  tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGTGGCTTGGGGGGCCGGGGGGGTTTTTTTTTTACCCGTGGGGAGGGGAACAGAGGATTGCACCTCCCCGcc                                                                                                    >  1:1030108/1‑151 (MQ=255)
                                                   cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                   <  2:235842/151‑1 (MQ=255)
                                                      tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTg                                                                                                  >  2:808810/1‑149 (MQ=255)
                                                          cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTg                                                                                            <  2:715889/151‑1 (MQ=255)
                                                          cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGAAGGGCAACAGCTGATTGCCCTTCACCCCCCGGccc                                                                                              >  1:924461/1‑149 (MQ=255)
                                                           gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTga                                                                                           <  1:178040/151‑1 (MQ=255)
                                                             tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGGAATGGGCGCCCGGGGGGGGTTTTTTTTTCCCAGGGAGAAGGGCGACAGCCGTGTGCCCTTCACCGCCCTGccc                                                                                              >  1:1035719/1‑146 (MQ=255)
                                                               tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTCTTTTTGCCAGGGAGACTGGCCAAAGCTGATTGCCCTTCACCCCCCGGCCCTgtgag                                                                                        >  1:1037760/1‑150 (MQ=255)
                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCCCTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCCgagag                                                                                        >  2:301374/1‑149 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCTgagag                                                                                        >  1:977759/1‑144 (MQ=255)
                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                 <  1:67928/151‑1 (MQ=255)
                                                                      gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcag                                                                                <  1:130687/151‑1 (MQ=255)
                                                                       ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTCCCAGGGAGACGGGGCACAGCTGATTGCCCCTCCCCCCCCGGCCCTGGGAGAGTTgca                                                                                 >  2:797370/1‑149 (MQ=255)
                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                              <  2:796888/151‑1 (MQ=255)
                                                                         aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                             <  1:152624/151‑1 (MQ=255)
                                                                           aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATTTGGGGGCCGGGTGGTTTTTTTTTTTACCCATGAGACGGGGAACAGCTGGTTGACCTTTACCGCCTGGCCCTGGGAGAGTTgtaagaa                                                                             >  1:917719/1‑143 (MQ=255)
                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGGCCTGAGAGGGTTGCAGCAAGCgg                                                                         >  1:280700/1‑151 (MQ=255)
                                                                               gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCCCCGCCCGGCCCTGCGAGAGGTGCAGCAAGCgg                                                                         >  1:771082/1‑149 (MQ=255)
                                                                               gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                         <  1:481623/149‑1 (MQ=255)
                                                                               gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                        <  2:748089/150‑1 (MQ=255)
                                                                                   tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTGGGGGGCAGGGGGGGTTTTTTTTTCACCAATGAGACGGGGAACAGCTGATTGTCCCTTCCCGCCTGGCCCCGGGAGAGGTGCAGCAAACGGGcc                                                                      >  2:270427/1‑148 (MQ=255)
                                                                                   tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGATTGGGCGCCAGGGGGGTTTTTTTTTTTACCCGTGAGACGGGCAACAGCTGATTTCCCCTCACCCCCCGGCCCCGCGAGAGGTGCCGCCAGCGGGcc                                                                      >  2:1197398/1‑148 (MQ=255)
                                                                                    gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                                  >  2:31388/1‑151 (MQ=255)
                                                                                      cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTg                                                                <  2:144429/151‑1 (MQ=255)
                                                                                        gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGGGAAGGGCCAAAGCTGATTGGCCTTTACCCCCTCGCCCCGGGAGAGGTGGAGCAAGCGGTCCACGCCgg                                                               >  1:1192084/1‑150 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTTCAGCAAGCGGGCCACGCTGGgtt                                                            >  2:243080/1‑150 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                           <  2:494788/151‑1 (MQ=255)
                                                                                            cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                          <  2:1238611/151‑1 (MQ=255)
                                                                                             aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGCAACAGCTGATTGGCCTTCACCCCCTGGGCCCGGGAGAGATGCAGCAAAAGGGCCACCCTGGTTTGcc                                                         >  2:1001666/1‑151 (MQ=255)
                                                                                                  tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGTTGGAGCAAACGGGCCCCGCTGGGGTGcccc                                                       >  2:993701/1‑148 (MQ=255)
                                                                                                  tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAACTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGGTGCAGCAAACGGGCCACGCCGGGTTGCCCccgc                                                    >  2:1172782/1‑151 (MQ=255)
                                                                                                   gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                   <  2:811101/151‑1 (MQ=255)
                                                                                                   gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagc                                                    <  2:134230/150‑1 (MQ=255)
                                                                                                        ggCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGCCCGGGGGGGTTTTTTTTTCACCAATGAGAAGGGCAACAGCTGATTGCCCCTCACCGCCTGGCCCTGGGGGAGTTGCAGCAAGCGGGCCACGCTGGTTTGCCCCAGCAgg                                                 >  2:499710/1‑148 (MQ=255)
                                                                                                           cAACGCGCGGGGAGAGGCGGTTTTGGGATTTGGGGGCCGGGGGGGTTTTTTTTTTACCCCGGGGAAGGGCAACAGGTGATTTGCCCTCACCCCCCCGCCCCGAGAGAGGTGTAACAAACGGGGCCCCCCGGGTTGCCCCCGCCgggggaa                                            >  1:488901/1‑145 (MQ=255)
                                                                                                           cAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                           <  1:643944/151‑1 (MQ=255)
                                                                                                              cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                                        <  2:29038/151‑1 (MQ=255)
                                                                                                                   ggggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                                   <  2:257670/151‑1 (MQ=255)
                                                                                                                      gagaGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                                <  1:135942/151‑1 (MQ=255)
                                                                                                                       agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                               <  2:361639/151‑1 (MQ=255)
                                                                                                                             ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                         >  1:1032637/1‑151 (MQ=255)
                                                                                                                                 ttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                       <  1:113774/150‑1 (MQ=255)
                                                                                                                                   gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                    <  2:186584/151‑1 (MQ=255)
                                                                                                                                    cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                   <  1:808798/151‑1 (MQ=255)
                                                                                                                                      tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                 >  1:919309/1‑151 (MQ=255)
                                                                                                                                          gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat              <  2:851527/150‑1 (MQ=255)
                                                                                                                                                 aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa        >  1:245148/1‑147 (MQ=255)
                                                                                                                                                  gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa     <  2:4959/151‑1 (MQ=255)
                                                                                                                                                   ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg    <  1:1095389/151‑1 (MQ=255)
                                                                                                                                                     tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt  >  2:1166915/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676

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: