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

ATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGAT  >  NC_000913/2151578‑2151911
                                                                                                                  |                                                                                                                                                                                                                                   
aTTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCCCCCTGAGAAACTCACTCTTCCCCCAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCCGGCTGGGCCAAAACTGCTTTTACAATTCAGGTAATAAAACGAgt                                                                                                   >  1:842521/1‑239 (MQ=255)
   gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTTAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCCGTGCCAAAGTTTTTTCATTATTCGTACGCCCGGTTAACCCCCGGCCGGACAAAAAATTCGTTTACAATTTAGGGGATAAAtacaggtgt                                                                                                 >  2:248226/1‑231 (MQ=255)
    ttcttttATGGGTCACCTATCAAAAAAAAGAAGCCGACGAGAAAAGAGCAAAATTTGTTCGCAAACCACCTTTTATTTTTTTTTTTTAAACGGGGGGTTGATTTTTTCTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAgcgc                                                                                                                                                                                              <  2:881304/143‑1 (MQ=255)
              ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTTAAACGCCCTGAGAAAATCACTCCTCCCCCCGCCCCAGTGCCAACGGTTTTTTATCATCCCTACGCTCGTTTATCGCCGGTGCGGGAGAAAAATGCCTTATAAAATACAGGTATAAATCGAGTTGTTCCTGccgtgg                                                                                     >  2:1209329/1‑234 (MQ=255)
              ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCAAc                                                                                                                                                                                    >  1:881304/1‑147 (MQ=255)
                    cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTCTTCAAACCCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACg                                                                                                                                                                                   >  1:695524/1‑142 (MQ=255)
                                   cgcagaagaGCAAAGAACAATAGTTGTTCGATAATCCGCTGTTATATTTATTGTAAAAACGGGGGGTTGAATTGTTTTTTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                              <  1:477991/232‑1 (MQ=255)
                                     cggacgaGAAAAGAAAAAATTTAGTTCTCAAATCAGCTTTTATTTTTTTTGTTAAACCGGGGGTTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                              <  1:64596/235‑1 (MQ=255)
                                        acgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCCCAGCGCCAGGGCCAACGGTTTTTCATCAACCGTACGCTCGGTTACCGCCCGGCTGGACAAAAAATGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTTCGTTTCTGCCGCCAGTCCACTTTAAAAACggg                                                          >  2:738822/1‑239 (MQ=255)
                                         cgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAAATCACCCTGCCCCCAGCGCCACTGCCAACGGTTTTTCATCATCCGTACGCCCCTTTACCGCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGGAATAAATCGCAGTTTTTCTTCGTTGCTGCCGCCAGTCCATTGGAAAAAAgggg                                                         >  2:64596/1‑239 (MQ=255)
                                                      aaaTTTGTTCCGCTAACCAGCTTTTATTTTTATTGTTAAACGGGGCGGTAAAATTGTGTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcac                                            <  1:706623/239‑1 (MQ=255)
                                                                          cttttatttttatttTTAAACCGGGCGTTTAATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                              <  1:915936/198‑1 (MQ=255)
                                                                                        ttttAAACGGGGGGGTTGATTTGTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg           <  1:738822/238‑1 (MQ=255)
                                                                                                 ggCGGTTGAATTTTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCAGGGAATAAATACAGTTGTTCCTTCGTTGCTGCCGCCAGTCCCTTTTAAAAAGGGGGAAAGAAATGCACATTAGGGGAAATATACAACACATTTTAAGCCCGTCGCttgga  >  2:641927/1‑234 (MQ=255)
                                                                                                  gCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCCACGGTTTTTCATCATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAAATGCTTTTACAATTCAGGGAATAAAACGAGTTGTTCCTGCGTTTCTGCCGCCAGTCCATTGAAAAACAGGGGAAAGAAAGGCACATTAGGGGAATTATCCAGCATATTTaa                 >  2:394269/1‑222 (MQ=255)
                                                                                                                  |                                                                                                                                                                                                                                   
ATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGAT  >  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: