Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,388,041 T→G T50P (ACG→CCG)  aaeB ← aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,0410TG91.7% 115.0 / ‑0.4 49T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (17/27);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTA  >  NC_000913/3387795‑3388284
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
caATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcg                                                                                                                                                                                                                                                   <  1:187624/251‑1 (MQ=255)
   tGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATACATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGGCCGGCCcggcg                                                                                                                                                                                                                                                >  2:211646/1‑251 (MQ=255)
     aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGCTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                              >  1:191368/1‑251 (MQ=255)
     aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                              >  1:321356/1‑251 (MQ=255)
             gggataACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAACGGCCGGACCGGCGGcaacaaacc                                                                                                                                                                                                                                      >  2:136094/4‑248 (MQ=255)
                        caacctcccacgcactctacttttCTAGTCGTCCTAGCAAGGAATTCCATGTAAAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGGGCGCTGATCATCGCAATGATGATCCCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGCCCGGCGGCATCAATCGCCGCTGTCAg                                                                                                                                                                                                                            <  2:278427/230‑1 (MQ=255)
                                   gCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCg                                                                                                                                                                                                                >  1:292560/1‑251 (MQ=255)
                                     aTACGAGTTTTCTATGCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCATCAATAATGGTGCGCGGATCATCGCAAGGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTg                                                                                                                                                                                                              <  2:292560/251‑1 (MQ=255)
                                                      gTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAg                                                                                                                                                                                             <  2:238210/251‑1 (MQ=255)
                                                             cGAGTATATCCAGGTACAAAATCCGGCCCTGATCCAGCACACCAGAATCATCATGACTGGTGCGCGGATCATCGCAATTATGAGCACCAGTCCGGCAATACAGCCAATAAATGGGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCGGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGCAAGTGAAAGc                                                                                                                                                                                      <  2:219000/251‑1 (MQ=255)
                                                                gagtatccatctACAAAAACCGTCCCAGAGACAGCAAACCAGAACCATCAAAAAGGGGTCGCGGATTGTCGAAATGATGACCACCAGTCCGTCAATACAGCCAATAAATGTGCCGATTATGCGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGTCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaa                                                                                                                                                                                   <  2:279575/240‑1 (MQ=255)
                                                                       cAGGTAAAAAAACCGGCCCAGAGAAAGCACACCAGAATCATCAAAAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAg                                                                                                                                                                            <  2:354841/251‑1 (MQ=255)
                                                                           tACAAAAACCGGCCCAGATACAGCACACCAGAATCATCGATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCTGGAAGTGAAAGCCAACAAACAACGcc                                                                                                                                                                        >  2:418460/1‑251 (MQ=255)
                                                                               aaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAACCAGCGCCAGt                                                                                                                                                                     >  1:212114/1‑250 (MQ=255)
                                                                                   ccGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                                <  1:496393/251‑1 (MQ=255)
                                                                                    cGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGGCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGAACAATg                                                                                                                                                               >  2:441665/1‑251 (MQ=255)
                                                                                     ggCCCAGATACAGCACACCAAAAACATAAATAATGGTGCGTGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGACGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                              <  2:131396/251‑1 (MQ=255)
                                                                                       cccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCACCAAACAGCGCCAGTACAATGGCg                                                                                                                                                            >  1:213682/1‑251 (MQ=255)
                                                                                                        aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAATCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCCCGATAGCGAATAGCGACAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAACGAGGCGTTTCCAACTGGACGTGAAAGCCCACCAACAGCGACAGTACAATGGCGGTCGCCCGGTTTTACCg                                                                                                                                          >  2:284776/1‑251 (MQ=255)
                                                                                                              atcaATAATGGTGCGCGGCTCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCACAAAAAGCCACGATAGCGAATAGCGCCAGACTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATAGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCCACAAACAGCGCCAGTACAATGGCGGGCGACACGTTTTACCCAAAGc                                                                                                                                    >  2:61990/1‑250 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGCCCGTa                                                                                                                                                                    >  2:164458/1‑216 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTa                                                                                                                                                                    <  1:164458/216‑1 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAata                                                                                                                               <  1:211646/251‑1 (MQ=255)
                                                                                                                     atGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAATATGt                                                                                                                            >  1:180450/1‑251 (MQ=255)
                                                                                                                      gggTGCGCGGATCATCGAAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGtt                                                                                                                           <  2:427268/250‑1 (MQ=255)
                                                                                                                              ggATCATCGCAATGATGATCACCATTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCa                                                                                                                   <  1:471081/251‑1 (MQ=255)
                                                                                                                              ggATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTTCCGCAAAGCGGATATGTTgggttgc                                                                                                                    >  1:253304/1‑245 (MQ=255)
                                                                                                                                   atcGCAAGGAGGATCACCAGGCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACTATAGCGAATAGCGCCAGAATACGTTTCACTTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCTGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGCTGTTTAGCAATGga                                                                                                              <  2:441544/251‑1 (MQ=255)
                                                                                                                                     cGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGgaga                                                                                                            <  1:506707/251‑1 (MQ=255)
                                                                                                                                                             aaTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCaaa                                                                                    <  2:180450/251‑1 (MQ=255)
                                                                                                                                                                aCAGCCAAGAAATGTGCCGATGATGCGTAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTc                                                                                 <  2:335183/251‑1 (MQ=255)
                                                                                                                                                                         aaaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCg                                                                        <  1:179379/251‑1 (MQ=255)
                                                                                                                                                                                  gatgatGCGCAAAAAGCCACGATAGCGAATATCGCCAGAATACGGTTCACCTACCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTTGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGTTTTTCCCCCAAACGAATATGTTGGTTTGTCATGGGGAAAATACCCATCGTGATTAACCAAACTCACCCACGCGCGGGGCCAt                                                               >  1:141456/1‑251 (MQ=255)
                                                                                                                                                                                      atgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTa                                                           <  1:162434/251‑1 (MQ=255)
                                                                                                                                                                                        gcgcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGTCGATGGGCCATTTTAc                                                          <  2:253304/250‑1 (MQ=255)
                                                                                                                                                                                            aaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaa                                                      >  1:272395/1‑250 (MQ=255)
                                                                                                                                                                                              aaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGc                                                                                                                                                                                      <  1:274164/122‑1 (MQ=255)
                                                                                                                                                                                               aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAc                                                  <  1:112455/251‑1 (MQ=255)
                                                                                                                                                                                                aGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGGATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACACACGCAGGGGAAGGGCCATTTTACGGGAGAAc                                                  >  1:70657/1‑250 (MQ=255)
                                                                                                                                                                                                      gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTACCGCAAAGCGAAAATGTTGGTTAGCAATGGAGAAAAGACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAAAACGAATCCt                                           >  1:142791/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTAGCGCCTGAATACGGTTCCCCCCCCGCAGCAAAGGTCCGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGGTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGTCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                    <  2:479492/251‑1 (MQ=255)
                                                                                                                                                                                                              agagtgcCAGAATTCCGTTCACCCCCCGCAGCAAAGGTCGGACCGGGTGCCATAATCGCCGCTGTCAGTCCCGCCCAGCTGGGCGTGTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGTGCCTTTTTACGGAAGAACGAATCCTGGCTTtcg                                   <  2:12233/246‑1 (MQ=255)
                                                                                                                                                                                                                aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                                 <  2:421669/251‑1 (MQ=255)
                                                                                                                                                                                                                       gAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCACCAGTACAATGGCGGACGCCAGATTTACCGAAAAGCGAATATGGTGGTTAAAAATGGGGAAAATAACCAACGTGATTAACCAAAATACCGCAGGAGATGGGCCATTTTACAGAGAAACCAATCCTGACCTTCGTGGCGGTCt                          >  1:306932/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                     <  1:284776/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAg                <  1:441665/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccccagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGTAAGTGAAAGCCAACAAACAGCGCCAGTGCAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGTCCATTTTACGGAAGAACGAATCCGGGCTTTCGTCGCGATCTTGTTTGCCAGTGac            <  2:465745/246‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccccaccaacggCCGGACCGGCGGCAACAATCGCCGTTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAACCTACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCTAATATGTTGGTTAGCAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATGTTACGGAAGAACGAATCCTTGCTTTCGTCGTGATCTTGTGTGCCAGTGac            <  2:58976/239‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTa  <  1:263617/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTA  >  NC_000913/3387795‑3388284

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