Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,661 A→T 10.3% intergenic (‑59/+307) yhiI ← / ← yhiJ putative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,630,6610AT10.3% 294.3 / 20.7 124intergenic (‑59/+307)yhiI/yhiJputative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base A (47/62);  new base T (7/6);  total (56/68)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 125 total aligned reads displayed.

TGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACG  >  NC_000913/3630513‑3630806
                                                                                                                                                    |                                                                                                                                                 
tGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAacgac                                                                                                                                                  <  2:828778/150‑1 (MQ=255)
    gggCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTc                                                                                                                                              <  1:1295353/150‑1 (MQ=255)
       cgcAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGc                                                                                                                                           >  1:620163/1‑150 (MQ=255)
        gcAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCg                                                                                                                                          <  2:1250399/150‑1 (MQ=255)
         cAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACACCCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCgg                                                                                                                                         <  1:929928/150‑1 (MQ=255)
         cAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGCCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGCCGTTCAGCgg                                                                                                                                         >  1:1251576/1‑150 (MQ=255)
         cAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCgg                                                                                                                                         >  1:748545/1‑150 (MQ=255)
         cAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCAAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCgg                                                                                                                                         <  2:917281/150‑1 (MQ=255)
          aaCAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACCACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGCCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCTGAACGTCGTTcggcggc                                                                                                                                        <  1:1042786/150‑2 (MQ=255)
          aaCAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCTGAACGTCGTTcggcggc                                                                                                                                        <  1:628467/150‑2 (MQ=255)
           aCAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCGCTGTAAGCAAACGACGCGTGTTACAACTGCCGCTGAACGTCGTTcggcggca                                                                                                                                       <  2:949287/150‑3 (MQ=255)
             aGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACACCCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTCTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGAttc                                                                                                                                     >  1:1503444/1‑150 (MQ=255)
             aGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGAttc                                                                                                                                     <  1:1105570/150‑1 (MQ=255)
              cccacccAGCCACGATACCCGCCACCGCCTGTAACCCGACAACCCACCACCCCATCTGGCGCATACTCATATCCATAACCCACCCGCCATAATCCCAGTAAGCAAACCACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGAttct                                                                                                                                    <  1:853848/149‑1 (MQ=255)
               ccaccaCGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCTGAACGTCGTTcggcggcagttg                                                                                                                                   <  2:428944/150‑7 (MQ=255)
                caccacGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCTGAACGTCGTTcggcggcagttgt                                                                                                                                  <  2:787391/150‑8 (MQ=255)
                caccacGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGAttcttc                                                                                                                                  <  2:1716234/150‑1 (MQ=255)
                     cgccacgATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATcc                                                                                                                             <  1:367937/150‑1 (MQ=255)
                     cgccacgATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATcc                                                                                                                             <  1:563274/150‑1 (MQ=255)
                     cgccacgATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCGTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATcc                                                                                                                             >  1:1123937/1‑150 (MQ=255)
                          cgATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGAt                                                                                                                        >  2:257492/1‑150 (MQ=255)
                           gATAGCCGCCACCGCCAGTAACCCGACACCCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATa                                                                                                                       <  2:375004/150‑1 (MQ=255)
                           gATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGCCCCGCCATAATCCCTGTAAGCAAACGCCGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATa                                                                                                                       >  1:1606513/1‑150 (MQ=255)
                           gATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATa                                                                                                                       <  1:482224/150‑1 (MQ=255)
                                ccgccagcggcaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTTAGCAAACTACGCTTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaagaag                                                                                                                  <  2:234063/150‑1 (MQ=255)
                                ccgccaccgccaGTAACCCGACAACCCACCTCGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaagaag                                                                                                                  <  2:344060/150‑1 (MQ=255)
                                ccgccaccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaagaag                                                                                                                  <  1:1124628/150‑1 (MQ=255)
                                ccgccaccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaagaag                                                                                                                  <  2:1405539/150‑1 (MQ=255)
                                   ccaccgccaGTAACCCTACAACCCACCACGCCTGATGGCGCTTACTCTTATCCATAATCGTCCCGCCATAATCCCTGTAAGCAAACGACGTGTGTTACAACTGCCGCCGATCGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGaa                                                                                                               <  1:470159/150‑1 (MQ=255)
                                    caccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGaaa                                                                                                              >  2:1054841/1‑150 (MQ=255)
                                     accgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAAc                                                                                                             <  2:134086/150‑1 (MQ=255)
                                      ccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATTATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGTCGTTCAGCGGATTCTTCAATAAGGGATTAGAAGGAAACg                                                                                                            >  2:790555/1‑150 (MQ=255)
                                       cgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACgg                                                                                                           >  2:1021619/1‑150 (MQ=255)
                                         ccaGTAACCCGACCACCCACCACGCCAGATGGCGCTTCCTCTTATCCATACTCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCa                                                                                                         <  2:856808/150‑1 (MQ=255)
                                         ccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCa                                                                                                         <  2:1341548/150‑1 (MQ=255)
                                         ccaGTAACCCGACAACCAACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCa                                                                                                         <  2:857095/150‑1 (MQ=255)
                                           tgTAACCCGACTAACAACCACGCCAGATGGCGCTTAATCTTATCAATAATCGACCCGCCATAATCCCTGTAAGCAATCGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGAAAAGAAGGAAACGGCaaa                                                                                                       <  2:185893/149‑1 (MQ=255)
                                           aGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGCCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGACCGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCaaa                                                                                                       >  1:1624665/1‑150 (MQ=255)
                                                cccGACATCTCACCACTCCAGATGGCGCTTACTCTTATCCATAATCGACCCGACATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGa                                                                                                  <  1:1635843/150‑1 (MQ=255)
                                                cccGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTTTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGa                                                                                                  <  2:980476/150‑1 (MQ=255)
                                                cccGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGa                                                                                                  >  2:1402470/1‑150 (MQ=255)
                                                 ccGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGAc                                                                                                 <  1:616945/150‑1 (MQ=255)
                                                 ccGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAGAAGGAAACGGCAAATACGAcc                                                                                                >  1:1142452/1‑150 (MQ=255)
                                                       aCCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGAACCGCCATAATCCCTGTAAGCATAAGACGCGTGTTTCAAATGACGCCGAACGACGTACAGAGGATTCTTACATCCGGGATAATAAGGAATCGGCAAATACGACCAgagc                                                                                           >  2:310759/1‑150 (MQ=255)
                                                        cccACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAgcgcg                                                                                          >  2:846523/1‑150 (MQ=255)
                                                         ccaccaCGCCAGATGGCGCTTACTCTTATCCATAATCGACCCTCCATAATCCCTGTAAGCTAAAGACTCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAgcgcgc                                                                                         <  1:161361/150‑1 (MQ=255)
                                                         ccaccaCGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAgcgcgc                                                                                         >  1:954525/1‑150 (MQ=255)
                                                          caccacGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCa                                                                                        <  1:846523/150‑1 (MQ=255)
                                                            ccacGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCaga                                                                                      <  1:1564125/150‑1 (MQ=255)
                                                                  aaGATGGAGCATACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTg                                                                                <  2:952701/149‑1 (MQ=255)
                                                                           cTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGCCGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGCGCATTGAGCGCCAAc                                                                       >  1:979946/1‑150 (MQ=255)
                                                                              aCTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACggg                                                                    <  1:686883/150‑1 (MQ=255)
                                                                               ctctTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACgggg                                                                   >  1:711169/1‑150 (MQ=255)
                                                                                 ctTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGaa                                                                 >  2:1668893/1‑150 (MQ=255)
                                                                                 ctTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGaa                                                                 >  2:147688/1‑150 (MQ=255)
                                                                                 ctTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGaa                                                                 >  2:967691/1‑150 (MQ=255)
                                                                                     tCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGCCGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAACCGGCAACTACGACCAGCGCGCAGAGCAGTGAGCGACAACGGGGAAGGCt                                                             >  1:261688/1‑150 (MQ=255)
                                                                                       cATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTa                                                           >  1:654809/1‑150 (MQ=255)
                                                                                        aTAATCGACCCGCCATAATCCCTGTACGCAAACGACGCGGGTTACAACTGCCGCCGACCGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                          <  2:1700361/150‑1 (MQ=255)
                                                                                        aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                          >  1:1443135/1‑150 (MQ=255)
                                                                                        aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                          <  1:1494326/150‑1 (MQ=255)
                                                                                         tAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGt                                                         <  1:1054841/150‑1 (MQ=255)
                                                                                             cGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGATACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGaaa                                                     <  2:1517403/150‑1 (MQ=255)
                                                                                                cccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACg                                                  >  1:59872/1‑150 (MQ=255)
                                                                                                 ccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACgg                                                 >  1:571523/1‑150 (MQ=255)
                                                                                                  cGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCATATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACggg                                                <  1:5066/150‑1 (MQ=255)
                                                                                                  cGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACggg                                                <  1:1021619/150‑1 (MQ=255)
                                                                                                   gCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgt                                               <  2:620163/150‑1 (MQ=255)
                                                                                                       tAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgtagta                                           >  1:560011/1‑150 (MQ=255)
                                                                                                           cccTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtgt                                       >  1:445425/1‑150 (MQ=255)
                                                                                                               gTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtc                                   <  2:730766/150‑1 (MQ=255)
                                                                                                                tAAGCAAACGACGAGTGTTACAACTGCCGCCGAACGAAGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtca                                  >  1:1236876/1‑150 (MQ=255)
                                                                                                                  accAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtcatc                                >  1:770757/3‑150 (MQ=255)
                                                                                                                       acgacgCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGAGTCTTCCATCCGGGATAAGAAGGAAACGGAAAATACGACCAGCGCTCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTAATCATCAg                           <  2:514609/150‑1 (MQ=255)
                                                                                                                         gacgCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCa                         >  2:23557/1‑150 (MQ=255)
                                                                                                                          acgCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCAt                        >  2:1726924/1‑150 (MQ=255)
                                                                                                                           cgcgTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCAtt                       <  2:1597911/150‑1 (MQ=255)
                                                                                                                            gcgTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGATCGCCTACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTc                      <  2:1503444/150‑1 (MQ=255)
                                                                                                                            gcgTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTc                      >  1:1337167/1‑150 (MQ=255)
                                                                                                                             cgTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTcc                     <  2:790425/150‑1 (MQ=255)
                                                                                                                              gtgtTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCt                    <  1:892422/150‑1 (MQ=255)
                                                                                                                                gtTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGATGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTgg                  <  2:368403/150‑1 (MQ=255)
                                                                                                                                     aaCTGCCGCCGAGCGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTa             <  1:1394500/150‑1 (MQ=255)
                                                                                                                                     aaCTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTa             <  2:1487745/150‑1 (MQ=255)
                                                                                                                                      aCTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTaaa           <  2:1142452/150‑1 (MQ=255)
                                                                                                                                        atccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAAt          >  1:1371389/3‑150 (MQ=255)
                                                                                                                                        atccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGGAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAAt          >  1:97907/3‑150 (MQ=255)
                                                                                                                                        atccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATg         >  1:1340667/3‑150 (MQ=255)
                                                                                                                                         tccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATg         >  1:687602/2‑150 (MQ=255)
                                                                                                                                         gccgccGAACGACGTGCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATg         <  2:1084808/150‑1 (MQ=255)
                                                                                                                                          ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCCACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCCTCATCAGCATTCCTGGCCGTAAATGa        >  1:1314315/1‑150 (MQ=255)
                                                                                                                                           cgccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAAGGGGTAGTAAGGTTGTCATCATCAGCATTCCTGGCCGTAAATgta       >  2:272634/1‑148 (MQ=255)
                                                                                                                                           cgccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGaa       >  2:1641756/1‑150 (MQ=255)
                                                                                                                                           cgccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGaa       >  2:1605465/1‑150 (MQ=255)
                                                                                                                                           cgccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGaa       <  1:323057/150‑1 (MQ=255)
                                                                                                                                             ccgAACGACGTTCAGCGTATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAac     <  1:668625/150‑1 (MQ=255)
                                                                                                                                             ccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAac     <  1:359517/150‑1 (MQ=255)
                                                                                                                                             ccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAac     >  1:872071/1‑150 (MQ=255)
                                                                                                                                              cgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAaca    >  2:1174070/1‑150 (MQ=255)
                                                                                                                                                aaCGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACg  <  1:1009862/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAACTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACG  >  NC_000913/3630513‑3630806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: