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

CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTA  >  NC_000913/3387795‑3388284
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
caATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACcg                                                                                                                                                                                                                                                   <  1:187624/251‑1 (MQ=255)
   tGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATACATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGG‑CCGGCCcggcg                                                                                                                                                                                                                                                >  2:211646/1‑251 (MQ=255)
     aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGCTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCggcggc                                                                                                                                                                                                                                              >  1:191368/1‑251 (MQ=255)
     aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAAAGG‑CCGGACCggcggc                                                                                                                                                                                                                                              >  1:321356/1‑251 (MQ=255)
             gggataACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAACGG‑CCGGACCGGCGGcaacaaacc                                                                                                                                                                                                                                      >  2:136094/4‑248 (MQ=255)
                        caacctcccacgcactctacttttCTAGTCGTCCTAGCAAGGAATTCCATGTAAAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGGGCGCTGATCATCGCAATGATGATCCCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGCCCGGCGGCATCAATCGCCGCTGTCAg                                                                                                                                                                                                                            <  2:278427/230‑1 (MQ=255)
                                   gCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCg                                                                                                                                                                                                                >  1:292560/1‑251 (MQ=255)
                                     aTACGAGTTTTCTATGCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCATCAATAATGGTGCGCGGATCATCGCAAGGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTg                                                                                                                                                                                                              <  2:292560/251‑1 (MQ=255)
                                                      gTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAg                                                                                                                                                                                             <  2:238210/251‑1 (MQ=255)
                                                             cGAGTATATCCAGGTACAAAATCCGGCCCTGATCCAGCACACCAGAATCATCATGACTGGTGCGCGGATCATCGCAATTATGAGCACCAGTCCGGCAATACAGCCAATAAATGGGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCGGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGCAAGTGAAAGc                                                                                                                                                                                      <  2:219000/251‑1 (MQ=255)
                                                                gagtatccatctACAAAAACCGTCCCAGAGACAGCAAACCAGAACCATCAAAAAGGGGTCGCGGATTGTCGAAATGATGACCACCAGTCCGTCAATACAGCCAATAAATGTGCCGATTATGCGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGTCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaa                                                                                                                                                                                   <  2:279575/240‑1 (MQ=255)
                                                                       cAGGTAAAAAAACCGGCCCAGAGAAAGCACACCAGAATCATCAAAAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAg                                                                                                                                                                            <  2:354841/251‑1 (MQ=255)
                                                                           tACAAAAACCGGCCCAGATACAGCACACCAGAATCATCGATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCTGGAAGTGAAAGCCAACAAACAACGcc                                                                                                                                                                        >  2:418460/1‑251 (MQ=255)
                                                                               aaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAACCAGCGCCAGt                                                                                                                                                                     >  1:212114/1‑250 (MQ=255)
                                                                                   ccGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                                <  1:496393/251‑1 (MQ=255)
                                                                                    cGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGGCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGAACAATg                                                                                                                                                               >  2:441665/1‑251 (MQ=255)
                                                                                     ggCCCAGATACAGCACACCAAAAACATAAATAATGGTGCGTGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGACGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                              <  2:131396/251‑1 (MQ=255)
                                                                                       cccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCACCAAACAGCGCCAGTACAATGGCg                                                                                                                                                            >  1:213682/1‑251 (MQ=255)
                                                                                                        aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAATCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCCCGATAGCGAATAGCGACAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAACGAGGCGTTTCCAACTGGACGTGAAAGCCCACCAACAGCGACAGTACAATGGCGGTCGCCCG‑GTTTTACCg                                                                                                                                          >  2:284776/1‑251 (MQ=255)
                                                                                                              atcaATAATGGTGCGCGGCTCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCACAAAAAGCCACGATAGCGAATAGCGCCAGACTACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATAGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCCACAAACAGCGCCAGTACAATGGCGGGCGAC‑ACGTTTTACC‑CAAAGc                                                                                                                                    >  2:61990/1‑250 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGCCCGTa                                                                                                                                                                    >  2:164458/1‑216 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTa                                                                                                                                                                    <  1:164458/216‑1 (MQ=255)
                                                                                                                  ataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAata                                                                                                                               <  1:211646/251‑1 (MQ=255)
                                                                                                                     atGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GGTTTACCGCAAAGCGAATATGt                                                                                                                            >  1:180450/1‑251 (MQ=255)
                                                                                                                      gggTGCGCGGATCATCGAAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGtt                                                                                                                           <  2:427268/250‑1 (MQ=255)
                                                                                                                              ggATCATCGCAATGATGATCACCATTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCa                                                                                                                   <  1:471081/251‑1 (MQ=255)
                                                                                                                              ggATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑C‑GGTTTTCCGCAAAGCGGATATGTTgggttgc                                                                                                                    >  1:253304/1‑245 (MQ=255)
                                                                                                                                   atcGCAAGGAGGATCACCAGGCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACTATAGCGAATAGCGCCAGAATACGTTTCACTTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCTGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGCTGTTTAGCAATGga                                                                                                              <  2:441544/251‑1 (MQ=255)
                                                                                                                                     cGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGgaga                                                                                                            <  1:506707/251‑1 (MQ=255)
                                                                                                                                                             aaTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCaaa                                                                                    <  2:180450/251‑1 (MQ=255)
                                                                                                                                                                aCAGCCAAGAAATGTGCCGATGATGCGTAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTc                                                                                 <  2:335183/251‑1 (MQ=255)
                                                                                                                                                                         aaaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCg                                                                        <  1:179379/251‑1 (MQ=255)
                                                                                                                                                                                  gatgatGCGCAAAAAGCCACGATAGCGAATATCGCCAGAATACGGTTCACCTACCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTTGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCC‑C‑GTTTTTCCCCCAAACGAATATGTTGGTTTGTCATGGGGAAAATACCCATCGTGATTAACCAAACTCACCCACGCGCGGGGCCAt                                                               >  1:141456/1‑251 (MQ=255)
                                                                                                                                                                                      atgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTa                                                           <  1:162434/251‑1 (MQ=255)
                                                                                                                                                                                        gcgcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGTCGATGGGCCATTTTAc                                                          <  2:253304/250‑1 (MQ=255)
                                                                                                                                                                                            aaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑C‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaa                                                      >  1:272395/1‑250 (MQ=255)
                                                                                                                                                                                              aaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGc                                                                                                                                                                                      <  1:274164/122‑1 (MQ=255)
                                                                                                                                                                                               aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAc                                                  <  1:112455/251‑1 (MQ=255)
                                                                                                                                                                                                aGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGGATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACACACGCAGGGGAAGGGCCATTTTACGGGAGAAc                                                  >  1:70657/1‑250 (MQ=255)
                                                                                                                                                                                                      gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCC‑A‑GTTTTACCGCAAAGCGAAAATGTTGGTTAGCAATGGAGAAAAGACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAAAACGAATCCt                                           >  1:142791/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTAGCGCCTGAATACGGTTCCCCCCCCGCAGCAAAGGTCCG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGGTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGTCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                    <  2:479492/251‑1 (MQ=255)
                                                                                                                                                                                                              agagtgcCAGAATTCCGTTCACCCCCCGCAGCAAAGG‑TCGGACCGGGTGCCATAATCGCCGCTGTCAGTCCCGCCCAGCTGGGCGTGTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGTGCCTTTTTACGGAAGAACGAATCCTGGCTTtcg                                   <  2:12233/246‑1 (MQ=255)
                                                                                                                                                                                                                aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                                 <  2:421669/251‑1 (MQ=255)
                                                                                                                                                                                                                       gAATACGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCACCAGTACAATGGCGGACGCC‑A‑GATTTACCGAAAAGCGAATATGGTGGTTAAAAATGGGGAAAATAACCAACGTGATTAACCAAAATACCGCAGGAGATGGGCCATTTTACAGAGAAACCAATCCTGACCTTCGTGGCGGTCt                          >  1:306932/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGGTTCACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                     <  1:284776/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cACCTCCCGCAGCAAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAg                <  1:441665/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccccagcaAAGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGTAAGTGAAAGCCAACAAACAGCGCCAGTGCAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGTCCATTTTACGGAAGAACGAATCCGGGCTTTCGTCGCGATCTTGTTTGCCAGTGac            <  2:465745/246‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccccaccaacgg‑CCGGACCGGCGGCAACAATCGCCGTTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAACCTACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCTAATATGTTGGTTAGCAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATGTTACGGAAGAACGAATCCTTGCTTTCGTCGTGATCTTGTGTGCCAGTGac            <  2:58976/239‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGG‑CCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTa  <  1:263617/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGG‑CCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑A‑GTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTA  >  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: ‑