Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,689 G→T 53.1% intergenic (+43/+22) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6890GT53.1% ‑5.1 / 3.8 13intergenic (+43/+22)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (5/2);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

ACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGC  >  NC_000913/2151486‑2151926
                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
atgtaaaatggggCAAGGAGTGAAAAAACAAAAAAGAAAATTTTCAGGCTGGTGGTTAACACGGGGGGGAAAAATTTTGCGGGGTTTTAACCTTTAATCTTTTATTGGTCCCCTTTCAAAAAAAAGAAGCCGCCGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCGTTTATTTTTATTGTTAAAACGGGCGTTTAAATTGTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCcaca                                                                                                                                                                                                             <  2:184374/228‑1 (MQ=255)
               agaggtgtAAAAAAAAAAAAAAAAATTTTTTAGGCTGTATTTTTTCGCGGGGGGAAAAAATTTGCCGTGGTTTATACCTTTGATATTTTGGTGTTCACCAACAAAAAAGAATAACCCGAGGAGCAAAGGAAAAAATTGTTGCTCTAACAACCGTTTTTTTTTTTTTTTAAAGTGGGGGGTTGATTGTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                                                                                                                                                                                                <  1:693269/234‑1 (MQ=255)
                                                                                               ttGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCCTTGGAAAATCACCCTTGCCCCAGCGCCCAGGCCCACGGTTTTTACACATCCCTACCCCCGGTTTCCGCCCGGCTGGGCAAAAACATCTTTAACCATTTAGGGaaaaaaatcatgt                                                                                                               >  2:614384/1‑232 (MQ=255)
                                                                                                tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTTCACCGCCATGAGAAAACCACTCTTTCCACAGCGCCACTGCCCAAGCTTTATCATCTTCCGTACGCTCCCTTACCTCCCGGGCGGGCACAAATACATTTAACAATTCAGGTAACAAATCgc                                                                                                                  >  2:827086/1‑234 (MQ=255)
                                                                                                  aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGTGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTTTCTTTAAACGCCATGGGAAACTCACTCTTCCCACAGCGCCAGTTCCAACGTCTTATCATTATTCTGACCCCCGTTTTTCGCCGGGGTGTGCAAAAACAGCTTTTACAATTCAGGCAATAAAtatcctttg                                                                                                             >  2:832377/1‑230 (MQ=255)
                                                                                                       ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGAATGTTGTCTTCAACCCCCCTGAGAAACTCACTCTTCCCCCAACGCCAGTGCCAACAGTTTTTTATCATCCGTACACCCGTTTACCCCCCGGGTGGGACAAAAACTCTTTTAACAATTCGGGAaaaaaaaacgagtttcccgc                                                                                                        >  1:226889/1‑225 (MQ=255)
                                                                                                                aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATTAGAAAATCACTCTGCCCCCAGCGCCCGTGCCAACGGTTTATTCTCATCCCTACGCCCGTTTACCGCCCGGCCGGGCAAAAACTGGTTTTACAATTCAGGTAATAAAACGAGGTGTTTCTCGGTtgctgc                                                                                               >  2:718140/1‑238 (MQ=255)
                                                                                                                 ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTCAACCGCCATTAGGAAATCACTCTGCCCCCAACGCCAGTGCCCAAACTTTATCATCATCCCTCCGCCCGCTTTTCGCCACGGCGGCAAAAAAATACTTTTACAAATACTGGAAAAAAACGAGATGTTCATACGTTGTTgccg                                                                                             >  1:130752/1‑239 (MQ=255)
                                                                                                                  cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTTAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGGTTACCGCCGGGCTGGGCAAAAACTGCTTTTACAATTTAGGTAATAAATAGAGTTGTTTCTGCGGTTCTgccg                                                                                             >  2:135487/1‑238 (MQ=255)
                                                                                                                                                         tgttcgcaaaatcgcctttttttttattTTTTAAAGTGGGGGTTTAATTGTTTTTTTAAAACCCCCATGAAAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                      <  2:778136/217‑1 (MQ=255)
                                                                                                                                                            tCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTCTTCAACCTCCCTGAGAAACACACACTGCCCCCAGCGCCCGTGCCCACGTTTTTTCATCATCCGTACCCCCGCTTACCTCCGGGCTGGCCAAAAAATTCTTTAACAAATCACGTAATAAAAAGCGTTGTTCCTGCGTTTCTGCGCCCAGTCCCTTGAAAAAAGGGGGaaaaaaaaaaaaaaagggaga                                                  >  1:856347/1‑225 (MQ=255)
                                                                                                                                                               cTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCCACCGCCATGAGAAAATCACTCTTCCCACAGCCCCCAGGCCAAACTTTTTTCATCATCCGCAAGCTCGTTTTCCGCCCGGCTGGCCAAAAAATGCTTTTAAAAATCACGTAATAAAACCGGTTGTTCCTGCGTTGGTGCCCCCCGCCCATTGAAAAAAGGGGGAAAAAAAGGAACAATAGGGGaaaa                                               >  1:778136/1‑237 (MQ=255)
                                                                                                                                                                                     tgttAAACTGGGCGTTTATTGTGTTTTTTGCAACCGGCATGGGAAACTCACTCTGCCCACAGCCTCAGTTCCAACGTTTTATCCACATCCGTACGCTTGGTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                          <  2:856347/237‑1 (MQ=255)
                                                                                                                                                                                                  ggtttaattttttgtGTCCACCCGCCATGGAAAACTCCCTCGGCCCCCAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATAc               <  2:596470/227‑1 (MQ=255)
                                                                                                                                                                                                               gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGCGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGTTACCATCCGTAAAtgc  >  2:413218/1‑238 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
ACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGC  >  NC_000913/2151486‑2151926

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