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

ACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGC‑GCGGGGATAAAATTTACGGCGG‑TTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TGTTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGC  >  NC_000913/2151486‑2151926
                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
atgtaaaatggggCAAGGAGTGAAAAAACAAAAAAGAAAAT‑TTTCAGGCTGGTGGTTAACACGGGGGGGAAAAATTTTGCGGGGT‑TTTAACCTTTAATCTTTTATTGGTCCCCTTTCAAAAAAAAGAAGCCGCCGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCGTTTATTTTTATTGTTAAAACG‑GGCGTTTAAAT‑TGTTGTCTT‑AAACCGCCATGAGAAACTCACTCTGCCcaca                                                                                                                                                                                                             <  2:184374/228‑1 (MQ=255)
               agaggtgtAAAAAAAAAAAAAAAAATTTTTTAGGCTGTATTTTTTCGC‑GGGGGGAAAAAATTTGCCGTGGTTTATACCTTTGATATTTTGGTGTTCACCAACAAAAAAGAATAACCCGAGGAGCAAAGGAAAAAATTGTTGCTCTAA‑CAACCGTTTTTTTTTTTTTTTAAAGTGGGGGGTTGATT‑GTTTGTCTC‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                                                                                                                                                                                                <  1:693269/234‑1 (MQ=255)
                                                                                               ttGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TGTTTTCTT‑CAACCGCCCTTGGAAAATCACCCTTGCCCCAGCGCCCAGGCCCACGGTTTTTACACATCCCTACCCCCGGTTT‑CCGCCCGGCTGGGCAAAAACATCTTTAACCATTTAGGGaaaaaaatcatgt                                                                                                               >  2:614384/1‑232 (MQ=255)
                                                                                                tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTTAAT‑TGTTGTCTT‑TCACCGCCATGAGAAAACCACTCTTTCCACAGCGCCACTGCCCAAGCTTTATCATCTTCCGTACGCTCCCTTA‑CCTCCCGGGCGGGCACAAATACATTTAACAATTCAGGTAACAAATCgc                                                                                                                  >  2:827086/1‑234 (MQ=255)
                                                                                                  aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGTGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGGAT‑TGTTTTCTT‑TAAACGCCATGGGAAACTCACTCTTCCCACAGCGCCAGTTCCAACGTCTTATCATTATTCTGACCCCCGTTTT‑TCGCCGGGGTGTGCAAAAACAGCTTTTACAATTCAGGCAATAAAtatcctttg                                                                                                             >  2:832377/1‑230 (MQ=255)
                                                                                                       ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTTGAA‑TGTTGTCTT‑CAACCCCCCTGAGAAACTCACTCTTCCCCCAACGCCAGTGCCAACAGTTTTTTATCATCCGTACACCCGTTTACCCCCCGGGTGGGACAAAAACTCTTTTAACAATTCGGGAaaaaaaaacgagtttcccgc                                                                                                        >  1:226889/1‑225 (MQ=255)
                                                                                                                aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TGTTTTCTT‑CAACCGCCATTAGAAAATCACTCTGCCCCCAGCGCCCGTGCCAACGGTTTATTCTCATCCCTACGCCCGTTTA‑CCGCCCGGCCGGGCAAAAACTGGTTTTACAATTCAGGTAATAAAACGAGGTGTTTCTCGGTtgctgc                                                                                               >  2:718140/1‑238 (MQ=255)
                                                                                                                 ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTTAAT‑TTTTTTCTT‑CAACCGCCATTAGGAAATCACTCTGCCCCCAACGCCAGTGCCCAAACTTTATCATCATCCCTCCGCCCGCTTT‑TCGCCACGGCGGCAAAAAAATACTTTTACAAATACTGGAAAAAAACGAGATGTTCATACGTTGTTgccg                                                                                             >  1:130752/1‑239 (MQ=255)
                                                                                                                  cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TGTTGTCTT‑CAACCGCCCTTAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGGTTA‑CCGCCGGGCTGGGCAAAAACTGCTTTTACAATTTAGGTAATAAATAGAGTTGTTTCTGCGGTTCTgccg                                                                                             >  2:135487/1‑238 (MQ=255)
                                                                                                                                                         tgttcgcaaaatcgcctttttttttattTTTTAAAGT‑GGGGGTTTAAT‑TGTTTTTTTAAAACCCCCATGAAAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                      <  2:778136/217‑1 (MQ=255)
                                                                                                                                                            tCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTTAAT‑TGTTTTCTT‑CAACCTCCCTGAGAAACACACACTGCCCCCAGCGCCCGTGCCCACGTTTTTTCATCATCCGTACCCCCGCTTA‑CCTCCGGGCTGGCCAAAAAATTCTTTAACAAATCACGTAATAAAAAGCGTTGTTCCTGCGTTTCTGCGCCCAGTCCCTTGAAAAAAGGGGGaaaaaaaaaaaaaaagggaga                                                  >  1:856347/1‑225 (MQ=255)
                                                                                                                                                               cTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TTTTGTCTT‑CCACCGCCATGAGAAAATCACTCTTCCCACAGCCCCCAGGCCAAACTTTTTTCATCATCCGCAAGCTCGTTTT‑CCGCCCGGCTGGCCAAAAAATGCTTTTAAAAATCACGTAATAAAACCGGTTGTTCCTGCGTTGGTGCCCCCCGCCCATTGAAAAAAGGGGGAAAAAAAGGAACAATAGGGGaaaa                                               >  1:778136/1‑237 (MQ=255)
                                                                                                                                                                                     tgttAAACT‑GGGCGTTTATTGTGTTTTTTG‑CAACCGGCATGGGAAACTCACTCTGCCCACAGCCTCAGTTCCAACGTTTTATCCACATCCGTACGCTTGGTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                          <  2:856347/237‑1 (MQ=255)
                                                                                                                                                                                                  ggtttaattttttgtGTC‑CACCCGCCATGGAAAACTCCCTCGGCCCCCAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATAc               <  2:596470/227‑1 (MQ=255)
                                                                                                                                                                                                               gtCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGCGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGTTACCATCCGTAAAtgc  >  2:413218/1‑238 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
ACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGC‑GCGGGGATAAAATTTACGGCGG‑TTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACT‑GGCGGTTGAAT‑TGTTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTA‑CCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGC  >  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: ‑