Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,189,980 A→G T156T (ACT→ACC phoP ← DNA‑binding transcriptional dual regulator PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,189,9800AG95.7% 72.8 / NA 24T156T (ACT→ACCphoPDNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (12/10);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189735‑1190220
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ccGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGgtgt                                                                                                                                                                                                                                              >  2:73141/1‑250 (MQ=255)
           gCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCgg                                                                                                                                                                                                                                  >  1:303335/1‑251 (MQ=255)
                aaaaaaaGACGCAGTCATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGAGGTCAGt                                                                                                                                                                                                                             >  2:357183/1‑251 (MQ=255)
                         cGCACTAATGTTTTCATCAGCGAAAGTCGAACATCTAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTTCCCCAATTTTTTTGCTCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCTGATAGAGTTGGAGCATTAACGTATCTTTGCTGACTTCTCTGCCATTCTGGCGTCGCAACTTTTCCCTAATGGTGTCTTCGAACGCGGTCAGTTTGATccct                                                                                                                                                                                                                    <  1:128108/251‑4 (MQ=255)
                            aGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCg                                                                                                                                                                                                                 <  1:357183/251‑1 (MQ=255)
                                             cgcAATTCGAACAGATATCCCTGGCCTCGAATGGTGGTAATCTCTTCTTGCGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACTGCAATGGTATGGCTTTCCCGCAGCTCTGCATCCGGATATAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTGCCATAATTGTGTATTCGAACGTGGTCATTTTGATCACTTCGTCATTAATAGATAATTc                                                                                                                                                                                                <  1:17568/251‑1 (MQ=255)
                                                                  tGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTg                                                                                                                                                                           >  1:50548/1‑251 (MQ=255)
                                                                     tcgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgagagagctacacctgaa                                                                                                                                                                         >  2:219469/2‑239 (MQ=255)
                                                                               tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGc                                                                                                                                                              >  1:201921/1‑251 (MQ=255)
                                                                                                 gATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGcc                                                                                                                                            >  1:400813/1‑251 (MQ=255)
                                                                                                         gCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGc                                                                                                                                     <  1:73141/250‑1 (MQ=255)
                                                                                                             gAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCATCGCCGCATCCGGATAGAGTTGTAGCATTAACGAATCTTTGCTGGCCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATATATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTTAAGCCAGACAGCTATTa                                                                                                                                <  2:297390/251‑1 (MQ=255)
                                                                                                                           aGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGGCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCGGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGc                                                                                                                  <  2:50548/251‑1 (MQ=255)
                                                                                                                                                                            ccGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGACATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTAGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCAATCACCTCTTCCATATGACACGGTTTAGTCACATa                                                                  >  2:247870/1‑250 (MQ=255)
                                                                                                                                                                                         gAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTCCCATAATGGTGTATTCGAACGCGGTAAGTTTGATAACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACATGAAACGGGGGGAGCGAAATGACATGTGAAGCCAGACCGCTATTACGCCGCATTAATGGCTGCATTCGCGACATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCAcggg                                                    >  1:12283/1‑248 (MQ=255)
                                                                                                                                                                                               taacgacacTTTGCTGACCACTTTGCCATTATTGCGTAACAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCAAGGCGAGAGAGATCAACCTGAAACGGGGGGAGAGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTATTGCCTGCATTCGCGCCATc                                                                                                  >  2:168029/9‑199 (MQ=255)
                                                                                                                                                                                               tAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATc                                                                                                  <  1:168029/199‑1 (MQ=255)
                                                                                                                                                                                                       cTTTGCTGACCACTTTGCCATTGTTGCGGATCAACGTTACCCTACTAGTGCATTCGCCCGCGGTAAGTCTGATCACTTCGTCATTAATAGATAATTTACGGCGAGAGAGATCAACCTGAAACGGGAGGAGCGCAATGAACGGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTCAATTCGCGCCATCACCTCTTCACTATGAACCGGTTTCGTCACATAATCATCAGCACCGCCACTTAATACTTCg                                      >  2:377780/1‑251 (MQ=255)
                                                                                                                                                                                                          tGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGAc                                    <  1:366544/250‑1 (MQ=255)
                                                                                                                                                                                                                        ccATTATTGCGTATCCACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGCGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTATTCACTTAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCt                     <  1:157983/251‑1 (MQ=255)
                                                                                                                                                                                                                          tttattGCGTAGCAACGGTTCCATAATGGTGTATTCGAACGCGTTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAACCGGGGGGAGCGAAATGAACTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATCTGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAACttt                   >  2:82948/2‑251 (MQ=255)
                                                                                                                                                                                                                          attattGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCttt                   <  1:2855/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ttCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATAcc  <  1:164580/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ccATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATAcc  >  1:39903/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
CCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189735‑1190220

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