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

AATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCAT  >  NC_000913/2151511‑2151916
                                                                                                                                                                                          |                                                                                                                                                                                                                                         
aactaagaaaggAAGATGTTACAGTCGGGATGTATGCGCGGGGGAATAAAATTGAGGTGTGTAATACATTTGTCCTTTAATTGCCCACCTAACAAAAAAAAAACCCCGCAGAGTTAAAAAAAAATTTGTTTCGCAAACCGGCCTTTATTTTTTTTTTTAAAACGGGGGTTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                                                                                                                                                                                           <  2:1304544/223‑1 (MQ=255)
                        gcgtggtatataCGCGCGGGTAAAAAAATTACGGGCGGTTATACATTTCATCTTTTTTTGGCCACCAAAACAAAAAAAAAAGCCGACGATAAAAAAAAAAATTTAGTCCCTAAACCACCCTTTATTTTTATTTTTAAACCGGGCGGTTAATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                      <  2:142996/201‑1 (MQ=255)
                                                                        gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTTAACCGCCCTGAGAAACTCAATCTTCCCCCAACGCCCGTGCCAACGGTTTATCATCATcccttc                                                                                                                                                                      >  2:178162/1‑173 (MQ=255)
                                                                        gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCCTGTGAAAACCACCCTCCCCACAGCGCCAGGGGCCACGTTTTATCACTATTCCTAACCCCGCTTAACCCCCGGCTGGGGAAAAAATGTCTTAACAATTTCGGtatga                                                                                                                   >  2:712428/1‑225 (MQ=255)
                                                                         aTCCTTTATTGGTCCCCTATCAAAAAGAAGAAGCCGACGAGATAAGAGCAAAATTAGTTCCTTAATCAGCCTTTATATTTTTTTTTTAAACGGGGGTTTAATTTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACg                                                                                                                                                                     <  1:178468/177‑1 (MQ=255)
                                                                                    ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgt                                                                                                      >  1:812693/1‑228 (MQ=255)
                                                                                     gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTTTTCAACCGCCCTTAGAAAATCACCCTTCCCCCAACCCCAGGGCCAACGTTTTTTCATAATCCCTACCCCCGTTTTCCGCCCGGCCGGCCCAAAACTGCTTTTATAAATCAGGTAATAAATCGAgtt                                                                                                        >  2:251585/1‑225 (MQ=255)
                                                                                      tCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTCACCGCCCTGAGAAACTCACTCTTCCCCCAGCGCCCGTGCCCAAGTTTTTTTCTCATCCGTACGCCCGGTTTCCGCCGGGGTGGGCAAAAAATTGTTTTACAATTCACGTCATAAATCgcggtttt                                                                                                      >  2:1755559/1‑220 (MQ=255)
                                                                                            aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTTTTTTCCACCGCCCTGAGAAACTCACTCTGCCCACAACGCCAGTGCCCACGTTTTATCATCATCCGCAAGCTCCTTTACCGCCCGGCTGGACCAAAACTGCTTTTACCAATCAAGGAATTAATCGAGTTGTTCCTGccg                                                                                              >  2:909302/1‑226 (MQ=255)
                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCCCCATGAGGAACTCACCCTGCCCACAACGCCCGTGGCAACGTTTTATCATCATCCGTACGCCCGGTTTTCGCCGGGCTGGGCAAAAAATGCCTTAAAAATTCCAGGAATAAATCGAGGTGTTTCTTCgg                                                                                              >  2:994434/1‑226 (MQ=255)
                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGGTTACCGCCggg                                                                                                                                                     >  1:698200/1‑174 (MQ=255)
                                                                                                       gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCAACCGc                                                                                                                                                                                                                               >  1:543144/1‑90 (MQ=255)
                                                                                                       gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTTCCCCCAAGGCCAGGGCCAACCTTTTATCATCATCCGTACGCCCGGTTACCGCCCGGCTGGGTAAAAAAAGGTTTTAACATTTAGGGAATTAATCAGAGTGTGTCTGCGTTGCAgccgccc                                                                                   >  2:922614/1‑227 (MQ=255)
                                                                                                       gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCCCGTTTTCCCCCGGGGCGGGCAAAAAATGCTTTTACAATTCAGGGGATAAAATCAGTTGTTTCTTCGGTGCTtccgccc                                                                                   >  2:300348/1‑227 (MQ=255)
                                                                                                        aagCCGAGGAGCAAAGACAAATATTAGTTCTCTATCCAGCTTTTATATTTTTTTTTTAAGCTGGGGGTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                           <  1:1323231/189‑1 (MQ=255)
                                                                                                                           aaaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGTTATAAAACGAGATGTTCATGCGTTTCTGCCGGCCGTGCCTTGaaaaaaaggggt                                                            >  2:1341281/1‑221 (MQ=255)
                                                                                                                            aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATTAGAAACTCACTCTTCCCACAGCGCCCGTGGCAACGTTTTTTCATCATCCCTACGCCCCTTTACCCCCCGGCGGGGCAAAAAATTCTTTAAAAATTCAGGGAATAAATCGATGCGGTTCCTCCTTTCTTGCGGCCGTCCACTTGAAAAATgggg                                                              >  2:1487159/1‑226 (MQ=255)
                                                                                                                            aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTTAATTTTTGTCTTCCACCGCCATTAGAAACTCACTCTACCCCCAGCGCCCGTGCCAACGGTTTTTCATCATCCGTACGCCCGTTTTCCGCCCGGGTGGACAAAAACAGCTTTAACAAATCTGGGAATAAATTGAAGTGTTTCCGCTGTGGTTTCGCCCGTTCATTGGAAAACcgggga                                                            >  2:1619972/1‑225 (MQ=255)
                                                                                                                             aaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAAATCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTTCCGCCAGTCCATTGGAAAACAGGGt                                                             >  2:624734/1‑225 (MQ=255)
                                                                                                                                  gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAAACCCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCAAGGTTTATCATCATCCGTAAGCCCGTTTACCGCCGGGCCGGGCCAAAAATGCTTTTACAATTCAGGt                                                                                                                       >  2:566090/1‑167 (MQ=255)
                                                                                                                                   ttCTCTAACCACCCTTTATTTTTTTTTTTAAACCGGGGGGTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                                                                                                                      <  1:567012/167‑1 (MQ=255)
                                                                                                                                   ttCGCTAATCCCCCTTTATTTTTTTTTTTAAACCTGGGGGTTAAATTTTTTTCTTCAACCGCCATGAGAAactcact                                                                                                                                                                                                                <  1:11147/77‑1 (MQ=255)
                                                                                                                                   ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCCTTAGAAACTTACTCTTCCCCCCGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTTCCGCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGGAATTAAACGAGGTGTTTCTGCGGTTCTTCCCCCAGTCCATTGGAAAACAGGGGAAAGAAAgg                                                    >  2:1710121/1‑228 (MQ=255)
                                                                                                                                       gataaccggctttATTTTTTTTTTTTAAACCGGGGGGTTAATTTGTTTTCTAAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcac                                                 <  1:1343686/217‑1 (MQ=255)
                                                                                                                                         gaTCAGCCTTTTATTTTTTTTTTAAACCGGGGCGTTTAATTTGTTTCCCACACCCGCCATGAGAAACTCCCTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACAt                                               <  2:811298/227‑1 (MQ=255)
                                                                                                                                            tcaagtttATTTTTTTTGTATAAGCGGGGGGTTAAATTTTTCTCTTCAACCGCCATGAGAAACTCCCTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                            <  1:234680/221‑1 (MQ=255)
                                                                                                                                            gccgccttAAATTTATGTGTTAAAGCGGGGGTTGAAATTTTTCTCTACAACCGCCATGAGAAACTCCCTCTGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                            <  1:1489817/222‑1 (MQ=255)
                                                                                                                                             aGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGGACGCTCGTTTACCGCCGGGGTGGGCAAAAAATGCTTTAACAATTc                                                                                                                           >  1:1175442/1‑152 (MQ=255)
                                                                                                                                              gCCTTTATTTTTTTTTTTAAAACTGGCGTTTAAATTTTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCa                                                                                                                          <  2:1173337/152‑1 (MQ=255)
                                                                                                                                                  ttATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCTGGGCAAAAAATGCTTTAACAATTTAGGTAATAAATCGAGTTGTTTCTGCGTTGCTGTCGCCAGGCCCTTGAAAAAACGGGTTAAGAAAGGCACATGAGGAGaaat                                     >  2:664043/1‑226 (MQ=255)
                                                                                                                                                  ttATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAAAGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAAAAtt                                                                                                                            >  1:205367/1‑146 (MQ=255)
                                                                                                                                                            tttttAAACCGGGGGTATAATTTGTTGTCTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTCCGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                              <  1:1758661/224‑1 (MQ=255)
                                                                                                                                                            tGTTTAAACTGGCGGTTGAATTTTTGTCTTTAACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGGTTTTTCATCATCCGTACGCTCCTTTACCCCCGGGCTGGACAAAAAATGCTTTAACAATTCAGGGAATAAAACGAGGTGTTCCTGCGGTGGTTCCGCCAGTCCATTGGAAAACCGGGTAAAGAAAGGGACATGAGGAGAATTAACCAGcctc                           >  2:1684113/1‑225 (MQ=255)
                                                                                                                                                            tGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCCACGTTTTTTTATCATCCGTACGCTCGTTTTCCGCCCGGGTGGGCCAAAAATGCTTTTACAATTTCGGGAATAAAACCGGGTGTTCCCTCGTTTCTTCCGCCCGGCCATTTAAAAAAGGGGGAAAGAAAGGCACATTAGGGGAAATTTcccatg                             >  2:234284/1‑222 (MQ=255)
                                                                                                                                                            tGTTTAAACTGGCGGGTGAATTGTTGTCTTCAACCGCCCTGAGAAAATCACTCTTCCCCCAGCGCCAGTGGCAAAGGTTTTTCATCATACCTAAGCTCGTTTTACGCCGGGCTGGACAAAAACTGCTTTAAAAATTCAGGGAATAAAACGAGTTGTTTCTGCGTTGCTGCCGCCAGTCCATTGGAAAAACGGGGAAAGAAAAGAGCATGAGGAGAACTATCACGCata                           >  2:1808128/1‑228 (MQ=255)
                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTTTCTTCAACCCCCATGAGAAAATCACTCTTCCCACAACGCCAGTGCCCACGTTTTATCATCATCCCTACGCCCGTTTACCGCCGGGCCGGACAAAAACTGCTTTTACAATTCCGGTAATAAATTGGGGTGTTTCCGCGTTGTTGCCCCCCGTCCCTTTTAAAAAAGGGGTAAGAAAAGCGACATAGGGGGATATATCAGCATAtt                         >  2:1840091/1‑227 (MQ=255)
                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAGGGGCAGTGcc                                                                                                                                                                                            >  1:459908/1‑69 (MQ=255)
                                                                                                                                                                tAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTTCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTAAGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTTAAAAAAGGGGTAAAGAAAGGCACATCAGGGGAATTTTCCAGCCAATTTa                       >  2:154254/1‑228 (MQ=255)
                                                                                                                                                                         cGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGCGTTGTTCCTTCGTTGCTGCCGCCAGTCCATTGAAAAAAGGGGTAAAGAAAAGACCATGAGGGGAATTATTCAGCATAATTactgttctc               >  2:561032/1‑220 (MQ=255)
                                                                                                                                                                            gTTTAATTTTTGCTCTCCACCCGCCATGAGAACCTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATc                                  <  1:561947/208‑1 (MQ=255)
                                                                                                                                                                                          gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATAATTAACGCCAGTCGCTGGATACCAt  >  1:892131/1‑228 (MQ=255)
                                                                                                                                                                                          |                                                                                                                                                                                                                                         
AATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCAT  >  NC_000913/2151511‑2151916

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