Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,748 G→A S404S (AGC→AGT)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,7480GA91.7% 51.1 / ‑3.4 24S404S (AGC→AGT) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (10/12);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT  >  NC_000913/6513‑6985
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
tGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGACGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCCTGCAGGCCATTATGATATCTGCCAGTTGCCAAATCAGCGGAAGACTTAGCAAGGGgccg                                                                                                                                                                                                                                >  1:314347/1‑251 (MQ=255)
  aTCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATAGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTtccgcc                                                                                                                                                                                                                              >  1:108151/1‑251 (MQ=255)
                           aaaCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAAAGCAATGCAGGCCATTCTGATATCTGCCAGTTGCCACATCAGAGGAAAACTTAGCAAGGTGCCGCCGATGACCGTTGCgcaagggc                                                                                                                                                                                                          >  1:219265/1‑239 (MQ=255)
                                  atcCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCCTGCAGGCCCTTATGATATCTGCCAGTTGCCACATCAACGGAAGACTTAGCAAGGAGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAac                                                                                                                                                                                              >  1:370291/1‑251 (MQ=255)
                                                 agaAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAg                                                                                                                                                                               <  1:166333/251‑1 (MQ=255)
                                                            ggccGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAgg                                                                                                                                                                    <  1:382790/249‑1 (MQ=255)
                                                                           aaCGCAACGGATCAAACACCTGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATGGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACAGCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAagag                                                                                                                                                      <  1:375262/250‑1 (MQ=255)
                                                                                        aatcgcagggagctcGCCGAGTTTACGCTGGAGTAGA‑AATGCGATGTCAATGGTATTAACCATATGTGAGAAAAGTAAAATTGCGTTCAAAATGGTAATAGCCATGCAGGCCATGATGAGATCTGTCAGTTTCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATTACCTTGGCGAAGGTGAAAATCCGCAAACACTTGACCGCTTTAGGGTTGTTCAGTCGTAAAAAGAATAGATTCCTTTCGGCaa                                                                                                                                        <  2:94067/237‑2 (MQ=255)
                                                                                                  cgaaCGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGAAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGtt                                                                                                                              <  1:280985/248‑1 (MQ=255)
                                                                                                                aCGCTGGCGTAGATAAA‑CACTGGCAAGGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGc                                                                                                                <  1:101925/251‑1 (MQ=255)
                                                                                                                      gCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGCGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCAAGCCCTTATGATAACTGCCAGTTGCCCAAGCAGGGGAAGACTTAGACAGGCGCCGCCGGTGACCGTTTCGAAGGTGCAGAGCCGCAAAAACCAGAACGCCTTAGGGTTGTTTAGGGGTACAATGTAGAAACTGTTTTCAGCGTAAATTTCGATGGCAACGATgagggtgtg                                                                                                            >  2:31168/1‑241 (MQ=255)
                                                                                                                        gTAGATAAT‑CACTGGCATGGGTATGAACCACAGGCTAGAGCATTAAAAGGGAGGTCAAATGGGTAATATCGATGCAGGGCATTATGATATCTGCCAGTGGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAGTGTGTTCGGCATAAATTTAGTTGGCAACGATGGAGCTGAAGGCa                                                                                                        <  2:141842/251‑1 (MQ=255)
                                                                                                                            aTAAT‑CACTGGCAATTGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGTTCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTGCTGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACa                                                                                                    <  1:141323/251‑1 (MQ=255)
                                                                                                                                                    ccACAGGCGAGAGCAGTAAAATGGCGGTCAAAATGGTAATAGCCATTCAGGACATTATGAAATCTGCCAGGTGCCACATCAAAGGGAGACTTAGCAAGGTTCCGCCGATGACCGTTGCGACGGAGCAGAACCGCAAACACCAACTCGCTTTTGGGGTGTTCAGTCGAAAAAAGAAGAGATTGATTTCGGCGTTAATCTAGTTTGCGAAGATGAAGATGAAGGCAAACAAAATAACCCCGAGGGTAACAAAc                                                                             >  2:49993/1‑251 (MQ=255)
                                                                                                                                                      acaGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTc                                                                           <  1:21518/251‑1 (MQ=255)
                                                                                                                                                          gCGAGAGCTGTAAAATGGCGGTCAAATTGGTAATGGCCATGCAGGCCATTAGGATATCTGCCAGTTGCCACATCAGCGGAAGACGTAGCAAGGTGCCGCCGATGACCGTTGCGAATGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGGAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCa                                                                       <  2:116872/251‑1 (MQ=255)
                                                                                                                                                          gCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGAGATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCa                                                                       <  1:83169/251‑1 (MQ=255)
                                                                                                                                                                   gTAAAATGGCTGGCAAATTTGGAATATCCTGGCAGGCCAGTATGATATCCGCGAGTTGCCACAGCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGTGTGTGCAGGCGTAAAAAGAAGAGATTGCTGTCGGCGTAAATGTAGTTGGCAACGAGGGAGCTGAAGGCAGACATAATAACCACGAGGGTGACAAAATCAGCAcctgaggaa                                                              <  2:162637/251‑8 (MQ=255)
                                                                                                                                                                          ggCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCAtt                                                        >  1:240974/1‑250 (MQ=255)
                                                                                                                                                                             ggTCAACTTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCCAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGc                                                                                                                                                                   >  2:294873/1‑140 (MQ=255)
                                                                                                                                                                             ggTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGc                                                                                                                                                                   <  1:294873/140‑1 (MQ=255)
                                                                                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATa                                <  1:365025/250‑1 (MQ=255)
                                                                                                                                                                                                    aGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCATAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGc                             >  1:2346/1‑251 (MQ=255)
                                                                                                                                                                                                                               acaTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGtt  >  1:351439/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             
TGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAAT‑CACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT  >  NC_000913/6513‑6985

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