Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 503,429 C→T intergenic (‑191/+47) ybaL ← / ← fsr putative transporter YbaL/fosmidomycin efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913503,4290CT93.8% 42.4 / ‑5.8 16intergenic (‑191/+47)ybaL/fsrputative transporter YbaL/fosmidomycin efflux pump
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/7);  minor base G (1/0);  total (9/7)
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.92e-01

AATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAG  >  NC_000913/503206‑503670
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
aaTGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGttt                                                                                                                                                                                                                        <  1:302491/251‑1 (MQ=255)
                                                 ggtggcTATTGTCCATTTTTGGCCGGGGAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGACTTATGGTGACTGTAAAACGATACAGCCAAAGTTTGCCTTTGGCTGCCATGAAATCAGTCTTAATGCCGGTTATCAGACAGGa                                                                                                                                                                        >  2:363253/1‑250 (MQ=255)
                                                   tggtTATTGTCCATTTTTGGACGGGAAAACCAAAATTACAGGTAATAGTTTATGACACTTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat                                                                                                                                                                      >  2:361690/1‑250 (MQ=255)
                                                     gtTATTGTCCATTTTTGGCCGGGAAAGCCAAAATTACAGGTAATAGTTTATGACAATTTCATTTATGATGGTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCTCGGATCGGAAAACAGTAGGCGGTTTAACAGAGGAAAGGCGTATGGTGACTGTAAAACGATGCAGCCAAAGTGTGCCTTTGGTTGCAACGAATTCAGTCTTTATGCCGGTTATCAGGCAGGAATATg                                                                                                                                                                   <  1:160886/251‑1 (MQ=255)
                                                                        cGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGTGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCTTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTg                                                                                                                                                <  1:363253/251‑1 (MQ=255)
                                                                                                 aGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGAATTTATAGAc                                                                                                                       >  2:63921/1‑251 (MQ=255)
                                                                                                     tATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAAc                                                                                                                   >  1:4160/1‑251 (MQ=255)
                                                                                                                cATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGGAAGCACAGATTTTATAGACTAACTCgctactggg                                                                                                        >  2:207522/1‑243 (MQ=255)
                                                                                                                        atgTTCATGACTAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCACAAAAAGAACGGATCGGAACACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGAGTTTATAGCCTAACTTGATGCAGGTGGGATCgg                                                                                                >  2:338577/1‑251 (MQ=255)
                                                                                                                                        attGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATTCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGGGGCAGGAAAGCACAGATTTTAGAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAAc                                                                                <  1:361690/251‑1 (MQ=255)
                                                                                                                                                   cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCa                                                                     <  1:113537/251‑1 (MQ=255)
                                                                                                                                                          gtGGAGTTAGCTGCAAAAAAAGCACGGATCGGCAGCCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGGGACGGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATAAGTCTTTATGCCGGTTATCAGGCCGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGAAGAGATTATAAAGAATAACTCGGTGACGGGGGGAGCGGCGACAAGACCCAGAAcaggcgcgccgaagcatcc                                                              >  2:95496/2‑233 (MQ=255)
                                                                                                                                                                                   ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAAATAGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaa                                       >  1:13510/1‑249 (MQ=255)
                                                                                                                                                                                         gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa                               <  1:207522/251‑1 (MQ=255)
                                                                                                                                                                                           aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGCTGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTcc                                >  1:342909/1‑248 (MQ=255)
                                                                                                                                                                                                                      ggCTTATGGTGACGGTAAAACGATGCAGCCAAAGTGTGCCTTTGGCTGCAATGAAATCAGTCTTTAGGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCAGACCGGTACGTCCTTGAagcag  <  2:18128/251‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
AATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAG  >  NC_000913/503206‑503670

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