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

AGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151504‑2151924
                                                                                                                                                                                           |                                                                                                                                                                                                                                            
aggaaaaaacaaaaaaaaaaGATATTTCAGAGTGGATGTATACGCGGGGGGAAAAAATTTACGGGGTTTATACCATTTATCCTTTATTGGTCCCCTATCAAAAAGAAGAAGCCGAGGAGTAAAGAGAAAATTTAGTTCGCAAATCAGCCGTTATATTTATTGTTTAACCGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                         <  1:362326/224‑1 (MQ=255)
     aaaaaaaaaaaagaaAATATTTCAGATTGGATGTATACGCGGGGGGAAAAAATTTAGGGCGGTTATCCCATTGATCTTTTATGGGTCACCTACCAAAAAAAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAACCAGCCGTTATATTTATTGTTTAAACGGGCGGTGGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatc                                                                                                                                                                                    <  1:829289/229‑1 (MQ=255)
                      tatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCAATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGtt                                                                                                                                                                  <  1:756705/239‑1 (MQ=255)
                                                                            tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCCGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACCAAAACTGCTTTTACAATTCAGGTAATAAAACGAgttg                                                                                                            >  2:864063/1‑239 (MQ=255)
                                                                            tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAAACGCCATAAGAAAATCACTCTGCCCACAGCGCCAATGCCCAAGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCCGGACACAAAACGCTTTAACAATTTCGGGAATAAATCGAgttg                                                                                                            >  1:1282206/1‑239 (MQ=255)
                                                                            tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttg                                                                                                            >  1:394252/1‑239 (MQ=255)
                                                                                   ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAATTTCGGGAATAAATCGAGTTGTTCCTGc                                                                                                     >  1:885751/1‑239 (MQ=255)
                                                                                         gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTTATTGTTGTCTTCAACACCCCTGAGAAACACACTATGTCCACAACGCCCGTGACAACGGTTTATCTTCCTCCGCTCCCTCGTTTACCGCCCGGCTAGACAAAAACTGCTTTCACCATTCAAGTAATAAAATGAGTTGTTCCTCCgct                                                                                                  >  2:713090/1‑234 (MQ=255)
                                                                                         gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTTTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTTTCATCATTCGTACCCTCGTTTACCCCCCGGCTGGACAAAAACTGATTTAACAATTCAGGTAATAAATTGAGTTGTTTCTGCGTtgct                                                                                               >  1:109452/1‑239 (MQ=255)
                                                                                            aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGGGCCAGTGCCCAAGTTTTATCATTATCCGTACGCTCGTTTACCGCCGGGCTGGACCAAAAATGCTTTAACAATTCAGGGAATAAATCGAGTTTGTTCTGCGTtgtttcc                                                                                            >  1:862143/1‑234 (MQ=255)
                                                                                                aacaaaaaaaaaaaGCCGACGAGTAAAGAAAAAAATTATTTCTATAATCCGCCGTTATTTTTTTTGTTAAACCGGGCGGTGAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCa                                                                                        <  1:184393/235‑1 (MQ=255)
                                                                                                   cAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTTTTTAACCCCCATGAGTAAACCACTCCTCCCACAGAGCCACTTCCCACGTTTTTTCATCATACGCACGCTCGTTTACCGCCGGGGCGGGCAAAAACCACCTTTAAAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTTCTgccgcccct                                                                                      >  2:1359688/1‑235 (MQ=255)
                                                                                                      aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAAAACACTCTGCCCCCAGCGCCAGGGGCAACGTTTTTTCATCATCCCTACGCCCGTTTACCGCCGGGCTGGGCAAAAAATGCTTTAAAAATTCAGGTAATAAAACCAGTTTGTCCTGCGTTTCTGCCGCCAGTccct                                                                                  >  2:1203337/1‑237 (MQ=255)
                                                                                                           gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCGCCCTGAGAAAATCCCTCTTCCCACAACGCCAAGGGCAACGTTTTATTATCATTCCAACGCCCGTTTTCCCCCGGGCTGGGCAAAAAATGCTTTAAAAATTCAGGTCATAAATCGAGTTGTTTCTCCGTTTCTGGCGCCCGTCCATTGtaa                                                                             >  2:85810/1‑239 (MQ=255)
                                                                                                           gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAAATCAGGTAATAAATAGAGTTGTTCCTGCGTTGCTGCCGCCCGTCCATTTaaa                                                                             >  2:1296350/1‑239 (MQ=255)
                                                                                                           gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGAAAAAACTGCTTTAACAAATCAGGTAATAAATCGAGTTGGTCCTGCGTTTCTGCCGCCAGTCCATTTaaa                                                                             >  1:1395009/1‑239 (MQ=255)
                                                                                                                        tAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCTGGGCAAAAACTGCTTTAACAATTCCGGGAATAAATCGAGTTGTTCCTTCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaaga                                                                >  2:45566/1‑239 (MQ=255)
                                                                                                                          agaaCAAAATTGTTCTGCTAATCAGCTTTTTTATTTTTTTTTAAACCGGGCGTTTGATGTGTGTTCTTAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaagaaa                                                              <  2:454590/239‑1 (MQ=255)
                                                                                                                            agCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaagaaag                                                             >  1:182505/1‑238 (MQ=255)
                                                                                                                                aaattttgtccctaaaTCAGCTTTTTTTTTTTTTGTTAAATGGGGGGTTTGATTTTTTTTTTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcac                                                         <  1:48619/225‑1 (MQ=255)
                                                                                                                                 aaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCCCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGca                                                          >  1:306779/1‑236 (MQ=255)
                                                                                                                                    tggTTCGCTAATCAGCCGTTATATTTATTTTTAAAACGGGGGGTGAAATTGTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:879522/237‑1 (MQ=255)
                                                                                                                                    tatttCGTTAATCGGCCGTTATTTTTTTTTTTAAAGCTGGCGTTAGAATGGTTTCCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  1:864063/236‑1 (MQ=255)
                                                                                                                                      tgttcgtcaaaccagctttATTTTTTTTGTTTAAACCGGGGGGTAGATTTGTTGTTCTAACCCGCCATGAGAAACTCCCTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTCCGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:306779/222‑1 (MQ=255)
                                                                                                                                      gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCGCCATTAGGAACTCACTCTGCCCCCAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTTAGGTAATAAATCGAGGTGTTCCTGCGTTGCTGCCGCCCGTCCCTTGAAAAACGGGGGAAAGGAAGGcaaactagga                                                  >  2:830084/1‑231 (MQ=255)
                                                                                                                                                    ctgTTATTTTTATTTTTTAACTGGGGGGTTAAATTGTTTCCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  1:190352/236‑1 (MQ=255)
                                                                                                                                                    cgtttATATTTATTTTTAAAACTGGCGGTTGAATTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  1:1296350/236‑1 (MQ=255)
                                                                                                                                                     ctttATTTTTTTTTTTTAGACTGGGGGTTGTATTTTTGTCTCACACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  1:327670/236‑1 (MQ=255)
                                                                                                                                                       ttATATTTATTGTTTAAACTGGCGGTTTCATTGTTGTCCTTAACCCCCATGAGAAAATCACTCTCCCCCCAGCGCCAGTGCCCACTTTTTATTATCATCCCTACGCTCGCTTTCCGCCCGGCAGGATAAAAACTTTTTTAACAATTCCGGGAATAAATCAAGTTGTTCCTGCGTTTCTGCCGCCAGTCCATTGCAAAAAGGGGTAAAGAAAGGGACATGAGGAGAATTAtaccgacaat                                  >  1:850208/1‑230 (MQ=255)
                                                                                                                                                       ttATATTTATTGTTTAAACTGGCGGTTGTATTGTTGTCTTCAACCGCCCAGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCAACAACCGTACGCTCGTTTAACGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTTCGTTGCTGCCGCCAGGCCATTTAAAAACGGGGTAAAGAAAAGCACATGAGGAGAATTTTCCACCatat                                  >  1:1074491/1‑239 (MQ=255)
                                                                                                                                                                                tttgttttgttgTCTCAAACCGCCATGAAAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCa           <  1:692740/232‑1 (MQ=255)
                                                                                                                                                                                           gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCGGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAt  >  2:992802/1‑236 (MQ=255)
                                                                                                                                                                                           gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAt  >  2:1084957/1‑236 (MQ=255)
                                                                                                                                                                                           |                                                                                                                                                                                                                                            
AGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151504‑2151924

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