Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,301,6360AC26.0% 108.9 / 43.3 73D481E (GAT→GAG) mdtOputative multidrug efflux pump subunit MdtO
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base C (6/13);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTCCAGCGCCACGCGTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGCCCCGCCAAAACGCAGC  >  NC_000913/4301393‑4301882
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cTCCAGCGCCACGCGTTGGCACATCTCTTCACAGGCCTTAAACGCCGCATGCAGCCCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCCGACAAAGGTGTAAATCACCGCGGACACCACGGAACCAATGATGATACCCAGGGCGCGATCGCGa                                                                                                                                                                                                                                                  >  2:400699/1‑250 (MQ=255)
 cccAGCGCCACGCGTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAAt                                                                                                                                                                                                                                                <  2:102319/250‑1 (MQ=255)
     gcgcCACGCGTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTcc                                                                                                                                                                                                                                            <  1:399737/251‑1 (MQ=255)
             cgTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGGCTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATTCCCAGCGTGCCAGCCAGTTTTTGCGACAGTATGCGCGCTTCTCTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTc                                                                                                                                                                                                                                    <  1:585985/251‑1 (MQ=255)
             cgTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTc                                                                                                                                                                                                                                    <  1:205744/251‑1 (MQ=255)
                             cacaGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAAGGTCCACCAGGTCGTACACTGGGCCaaaa                                                                                                                                                                                                                    <  2:287817/251‑1 (MQ=255)
                                aGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCTCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCTCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACg                                                                                                                                                                                                                 <  2:397508/251‑1 (MQ=255)
                                  gCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGtt                                                                                                                                                                                                               <  2:443421/251‑1 (MQ=255)
                                   cATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGACGTGACTTCCTGCTGGAGTGGGATCCGCATTACTTTCCTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAAACAAAGGAGTAAATCACCGCGGACACCACGGAACAACTGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGACGTACACTGGGCCAAACACGttt                                                                                                                                                                                                              >  1:335567/1‑251 (MQ=255)
                                          cgccgcATGCACACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAACCGTTTTCGAGCg                                                                                                                                                                                                       >  2:572335/1‑251 (MQ=255)
                                                 tGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATGCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGACGTACACTGGGCCAAAAACGTTTTCGCGCGGGGCg                                                                                                                                                                                                  >  2:114657/1‑249 (MQ=255)
                                                        cGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCACTGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGa                                                                                                                                                                                           >  2:317574/1‑249 (MQ=255)
                                                         gATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTT‑CGAGCGTGGCGAGCGCGAAggg                                                                                                                                                                                       >  2:302384/1‑250 (MQ=255)
                                                                aTTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGaccac                                                                                                                                                                                 <  2:177650/251‑1 (MQ=255)
                                                                 tttGCAGATAAGGGCGCAGCGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGACTTTCCACCAGGTCATACACTGGGCCAAAAACGTTTTCGAGCGTGTCGAGCGCGAAGGTGaccccc                                                                                                                                                                                >  2:34911/1‑248 (MQ=255)
                                                                  ttGCAGATAAGTGCGCATAGCCGTGACGTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGTTACCAATGATGATACCCAGGGCGCGCTCGCGAATTGCCACCAGGTCGTACACTGGGCCAAAAACGGTTTCGAGCGTGGCGGGCGCGAAGGTGaccgcc                                                                                                                                                                                <  2:281843/250‑4 (MQ=255)
                                                                       gATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTCACATACCCCGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGACTTTCCACCAGGTCGTACACTGGGCCACAAACGTTTTCGAGCGTGGCGAGCGCGAAAGTGACCACCATCtgt                                                                                                                                                                          >  1:430406/1‑251 (MQ=255)
                                                                         tAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGCGACCACCATCtgtgg                                                                                                                                                                        >  2:131426/1‑250 (MQ=255)
                                                                          aaGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCATCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCtgtgt                                                                                                                                                                        <  1:478670/250‑1 (MQ=255)
                                                                           aGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTccacca                                                                                                                                                                                                                                        <  2:458889/185‑1 (MQ=255)
                                                                           aGTGCGCAGAGCCGTGACTTCCTGCTGGCGAGGGATCCGCATTGCTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTccacca                                                                                                                                                                                                                                        >  1:458889/1‑185 (MQ=255)
                                                                                gcagtgccttgtcTTCCTGCTGGCGTGGGATCCGCAGTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCATGTCAGTCAAAGGTGTAAATCACCGCGGACTCCACGGTACCAATGATGATACCCAGTGCGCGATCGTTAAGTGCCGCCAGGTCGTGCACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCTCCATCTGTGTGCCGat                                                                                                                                                                  <  1:556653/238‑1 (MQ=255)
                                                                                gcAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTAACTTTCAGGCCAGACAACGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGc                                                                                                                                                                                                        >  1:576032/1‑212 (MQ=255)
                                                                                 cAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGatat                                                                                                                                                                <  2:142049/251‑1 (MQ=255)
                                                                                  agagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGaa                                                                                                                                                                                                                                                 >  1:6487/1‑169 (MQ=255)
                                                                                  agagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGaa                                                                                                                                                                                                                                                 <  2:6487/169‑1 (MQ=255)
                                                                                  agagACGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGatat                                                                                                                                                                >  1:90400/1‑250 (MQ=255)
                                                                                       cgcgACTTCCTGCTCGAGCGGGCTCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTccacccgg                                                                                                                                                                                                                                      >  2:68019/4‑172 (MQ=255)
                                                                                       cGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAgg                                                                                                                                                                                                                                      <  1:68019/175‑1 (MQ=255)
                                                                                         tGACTTCCTGCTGGCGTGGGCTCCGCCTTACTTTACTTAACATACCCCGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGACTTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATaagaag                                                                                                                                                          >  1:345406/1‑249 (MQ=255)
                                                                                          gACTTCCTGCTGGCGTGGGATCCGCATTACGTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGc                                                                                                                                                       <  2:125379/251‑1 (MQ=255)
                                                                                           aCTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCg                                                                                                                                                      <  1:500049/251‑1 (MQ=255)
                                                                                               ccTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTc                                                                                                                                                   <  1:130305/250‑1 (MQ=255)
                                                                                               ccTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCa                                                                                                                                                  <  2:141861/251‑1 (MQ=255)
                                                                                                ctgctgGCGTGGGATCCGCCTTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTc                                                                                                                                                   >  1:436893/1‑249 (MQ=255)
                                                                                                  gctgGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTTCCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCag                                                                                                                                                 <  1:572335/249‑1 (MQ=255)
                                                                                                    tgGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCt                                                                                                                                             >  1:247242/1‑251 (MQ=255)
                                                                                                        gTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGCTATAAGAAGAGCGTTCAGAGCtggtg                                                                                                                                         >  1:28799/1‑251 (MQ=255)
                                                                                                           ggATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACATTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGc                                                                                                                                       >  1:429345/1‑250 (MQ=255)
                                                                                                                         ttacttAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGGGTGCGCGCTTCACTTTCAGGCCAGACAACGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCACCCCGGTCGTACACTAGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATTAGAAGAGCGTTCCGAACTGGTGGCAATCCATGCGCCCa                                                                                                                         >  1:540309/1‑250 (MQ=255)
                                                                                                                           acttAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAAc                                                                                                                       <  1:234390/250‑1 (MQ=255)
                                                                                                                               aaCATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGaa                                                                                                                  <  2:80450/251‑1 (MQ=255)
                                                                                                                                 cATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGaaaa                                                                                                                <  1:285356/251‑1 (MQ=255)
                                                                                                                                  aTACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCAt                                                                                                                                >  2:459401/1‑234 (MQ=255)
                                                                                                                                  aTACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCAt                                                                                                                                <  1:459401/234‑1 (MQ=255)
                                                                                                                                    aCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTTGTGGCAATCCATGCGCCCAACAGGAAAATCg                                                                                                             <  1:539573/251‑1 (MQ=255)
                                                                                                                                    aCCCAGCGTGCCAGCCAGTTGTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGGAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGTAAAATc                                                                                                              <  2:367171/250‑1 (MQ=255)
                                                                                                                                             gccagccaGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGca                                                                                                    <  2:440765/251‑1 (MQ=255)
                                                                                                                                              ccagccagTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGcac                                                                                                   >  1:185190/1‑251 (MQ=255)
                                                                                                                                              ccagccagTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTg                                                                                                          >  2:121780/1‑244 (MQ=255)
                                                                                                                                              ccagccagTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTg                                                                                                          <  1:121780/244‑1 (MQ=255)
                                                                                                                                                        tttGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATCCCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGATTTCGAGCGTGGCGAGCGCGAAGGAGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCCCAAACAGCAAt                                                                                         >  2:488978/1‑251 (MQ=255)
                                                                                                                                                             ggCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGAc                                                                                    <  1:501702/251‑1 (MQ=255)
                                                                                                                                                              gccgtgtGCGCGCTTCAATATCAGGCCAGACCAAGGTGTAAATCAACGCGAACCCAACCGTACCACTGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGCGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCAAACAGGAAAATCgg                                                                                                            >  2:380665/4‑226 (MQ=255)
                                                                                                                                                              gCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATGTCCACCAGGTCGTACACTGTGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCgg                                                                                                            <  1:380665/226‑1 (MQ=255)
                                                                                                                                                                      cgcgcTTCCCTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACTCCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCTACAATATTGTc                                                                           <  1:387289/251‑1 (MQ=255)
                                                                                                                                                                      cgcgcTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGCCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCACTCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTc                                                                           >  1:572649/1‑251 (MQ=255)
                                                                                                                                                                                        cAGACAAATGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGGTCATAGCAGGTTGCAATCCATGCGCCCAACAGGAAAATCGGGGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAg                                                         <  1:370957/251‑1 (MQ=255)
                                                                                                                                                                                             aaaGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCg                                                    >  1:176836/1‑251 (MQ=255)
                                                                                                                                                                                                      aaTCACCGCGGACACCACGGGACCAATGATGATACCCAGGGCGCGATCGCGACTTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAAc                                            >  2:15689/1‑250 (MQ=255)
                                                                                                                                                                                                       aTCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAAGAACAg                                          <  1:96200/251‑1 (MQ=255)
                                                                                                                                                                                                                 acacCACGGTACCAAGGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGGTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATc                                <  1:24638/251‑1 (MQ=255)
                                                                                                                                                                                                                   accacGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATcgc                              <  1:187510/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccacGGTACCAATGATGATACCCAGGGCGCGCTCGCGAATGTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGGCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCTCCAGAATcgcg                             <  2:17193/251‑1 (MQ=255)
                                                                                                                                                                                                                      acGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGcgcc                           >  2:474791/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACCCTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCCTGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACACTATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCaaa                      >  1:96350/1‑249 (MQ=255)
                                                                                                                                                                                                                                  atgatACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGcccc               >  2:479682/1‑251 (MQ=255)
                                                                                                                                                                                                                                  atgatACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGcccc               >  1:98449/1‑251 (MQ=255)
                                                                                                                                                                                                                                  atgatACCCAGGGCGCGATCGCGAATTGCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGcccc               <  1:382062/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tgatACCCAGGGCGCGATCGCGAATTTCCACCAGGACGTACCCTGGGCCAAAAACGTTTTCGAGCGTGGCGCGCGAGAAGGTGACCACCATATGTGTGCCGCTATACGAAGAGCGTTCAGAGCTGGTGGCACTCCATGCGCCCACCAGGAAAATCGGTGCCAGCACCAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGACTCGCGCCGCAACAGGccccc              >  2:537946/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cccAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGCGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGCCCCGCCaa          >  2:404977/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cccAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGCCCCGCCaa          >  1:104749/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cgcgATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGCCCCGCCAAAACgcagc  <  1:251178/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CTCCAGCGCCACGCGTTGGCACATCTCTTCACAGGCATTAAACGCCGCATGCAGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGCGGACACCACGGTACCAATGATGATACCCAGGGCGCGATCGCGAATTTCCACCAGGTCGTACACTGGGCCAAAAACGTTTTCGAGCGTGGCGAGCGCGAAGGTGACCACCATCTGTGTGCCGATATAAGAAGAGCGTTCAGAGCTGGTGGCAATCCATGCGCCCAACAGGAAAATCGGTGCCAGCACAAACAGCAATTCGACAATATTGTCCAGCCAGGGCATGACCAGTAGCGTGAATAACAGCGCCAGAATCGCGCCGCAAAAGGCCCCGCCAAAACGCAGC  >  NC_000913/4301393‑4301882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑