Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,689 G→T 50.3% 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,6890GT50.3% ‑4.7 / 3.2 15intergenic (+43/+22)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (4/3);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑GAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGAT  >  NC_000913/2151578‑2151911
                                                                                                                  |                                                                                                                                                                                                                                   
aTTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTGTCTTCAACCCCCCTGAGAAACTCACTCTTCCCCCAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCCG‑GCTGGGCCAAAACTGCTTTT‑ACAATTCAGGTAATAAAACG‑Agt                                                                                                   >  1:842521/1‑239 (MQ=255)
   gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTTTTTTTAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCCGTGCCAAAGTTTTTTCATTATTCGTACGCCCGGTTAACCCCCG‑GCCGGACAAAAAATTCGTTT‑ACAATTTAGGGGATAAAtacaggtgt                                                                                                 >  2:248226/1‑231 (MQ=255)
    ttcttttATGGGTCACCTATCAAAAAAAAGAAGCCGACGAGAAAAGAGCAAAATTTGTT‑CGCAAACCACCTTTTATTTTTTTTTTTTAAACGGGG‑GGTTGAT‑TTTTTCTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAgcgc                                                                                                                                                                                              <  2:881304/143‑1 (MQ=255)
              ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTTTTTTTAAACGCCCTGAGAAAATCACTCCTCCCCCCGCCCCAGTGCCAACGGTTTTTTATCATCCCTACGCTCGTTTATCGCCGGTGCGGGAGAAAAA‑TGCCTTATAAAATACAGGT‑ATAAATCG‑AGTTGTTCCTGccgtgg                                                                                     >  2:1209329/1‑234 (MQ=255)
              ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCAAc                                                                                                                                                                                    >  1:881304/1‑147 (MQ=255)
                    cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTTAA‑TTGTTTTCTTCAAACCCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACg                                                                                                                                                                                   >  1:695524/1‑142 (MQ=255)
                                   cgcagaagaGCAAAGAACAATAGTTGTT‑CGATAATCCGCTGTTATATTTATTGTAAAAACGGGG‑GGTTGAA‑TTGTTTTTTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑Gaaaaa                                                              <  1:477991/232‑1 (MQ=255)
                                     cggacgaGAAAAGAAAAAATTTAGTT‑CTCAAATCAGCTTTTATTTTTTTTG‑TTAAACCGGG‑GGTTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑Gaaaaa                                                              <  1:64596/235‑1 (MQ=255)
                                        acgaGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑ATGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCCCAGCGCCAGGGCCAACGGTTTTTCATCAACCGTACGCTCGGTTACCGCCCG‑GCTGGACAAAAAATGCTTTT‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTTCGTTTCTGCCGCCAGTCCACTT‑TAAAAACggg                                                          >  2:738822/1‑239 (MQ=255)
                                         cgaGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTTTTTTCAACCGCCATGAGAAAATCACCCTGCCCCCAGCGCCACTGCCAACGGTTTTTCATCATCCGTACGCCCCTTTACCGCCGG‑GCTGGGCAAAAACTGCTTTT‑ACAATTCAGGGAATAAATCGCAGTT‑TTTCTTCGTTGCTGCCGCCAGTCCA‑TTGGAAAAAAgggg                                                         >  2:64596/1‑239 (MQ=255)
                                                      aaaTTTGTTCCGCTAACCAGCTTTTATTTTTATTGTTAAACGGGGC‑GGTAAAA‑TTGTGTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑GAAAAACGGGGTAAAGAAAGGcac                                            <  1:706623/239‑1 (MQ=255)
                                                                          cttttatttttatttTTAAACCGGGC‑GTTTAAT‑TTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑Gaaaaa                                                              <  1:915936/198‑1 (MQ=255)
                                                                                        ttttAAACGGGGGGGTTGAT‑TTGTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑GAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg           <  1:738822/238‑1 (MQ=255)
                                                                                                 ggC‑GGTTGAA‑TTTTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCG‑GCTGGACAAAAACTGCTTTT‑ACAATTCAGGGAATAAATAC‑AGTTGTTCCTTCGTTGCTGCCGCCAGTCCC‑TT‑TTAAAAAGGGGGAAAGAAATGCACATTAGGGGAAATATACAACACATTTTAAGCCCGTCGCttgga  >  2:641927/1‑234 (MQ=255)
                                                                                                  gC‑GGTTGAA‑TTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCCACGGTTTTTCATCATCCGTACGCTCGTTTACCGCCCG‑GCTGGACAAAAAATGCTTTT‑ACAATTCAGGGAATAAAACG‑AGTTGTTCCTGCGTTTCTGCCGCCAGTCCA‑TT‑GAAAAACAGGGGAAAGAAAGGCACATTAGGGGAATTATCCAGCATATTTaa                 >  2:394269/1‑222 (MQ=255)
                                                                                                                  |                                                                                                                                                                                                                                   
ATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTT‑CGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAA‑TTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGG‑GCTGGACAAAAACTGCTTTA‑ACAATTCAGGTAATAAATCG‑AGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TT‑GAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGAT  >  NC_000913/2151578‑2151911

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