Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 367,852 (GC)6→5 coding (568‑569/834 nt) mhpR ← mhp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913367,8420G.87.5% 47.7 / 3.8 16coding (579/834 nt)mhpRmhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (7/7);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01
*NC_000913367,8430C.87.5% 48.1 / 4.0 16coding (578/834 nt)mhpRmhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (7/7);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGC  >  NC_000913/367708‑367965
                                                                                                                                      ||                                                                                                                          
aTGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcg                                                                                                                         <  2:145751/139‑1 (MQ=255)
               gCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCt                                                                                                        >  1:503342/1‑139 (MQ=255)
               gCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCt                                                                                                        >  2:51473/1‑139 (MQ=255)
                ccATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCtt                                                                                                       >  2:106473/1‑139 (MQ=255)
                            gATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTaatg                                                                                           >  1:50007/1‑137 (MQ=255)
                             aTTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgc                                                                                                                      <  1:544896/113‑1 (MQ=255)
                             aTTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgc                                                                                                                      >  2:544896/1‑113 (MQ=255)
                                            cACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTg                                                                           <  1:312778/139‑1 (MQ=255)
                                             aCCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                                                                                  <  1:704673/131‑1 (MQ=255)
                                             aCCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATATTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                                                                                  >  2:704673/1‑131 (MQ=255)
                                                             gcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCGCTTT‑‑GCGCGCGCGCCAGAAGGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGc                                                          <  2:68949/139‑1 (MQ=255)
                                                                 cGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGa                                                      <  2:246253/139‑1 (MQ=255)
                                                                                          tcctccTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGTCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTg                             <  2:373641/139‑1 (MQ=255)
                                                                                          tcctccTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAATGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTg                             <  2:50007/139‑1 (MQ=255)
                                                                                              cctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGt                         >  2:220635/1‑139 (MQ=255)
                                                                                                       aGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTc                  >  1:609777/1‑139 (MQ=255)
                                                                                                                       tCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgccagc  <  2:609777/139‑1 (MQ=255)
                                                                                                                           tAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAgg       <  2:31440/128‑1 (MQ=255)
                                                                                                                           tAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAgg       >  1:31440/1‑128 (MQ=255)
                                                                                                                                      ||                                                                                                                          
ATGGTCATTGCGCTCGCCATATACACCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGC  >  NC_000913/367708‑367965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: