Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,719 C→T P196P (CCG→CCA)  arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,7190CT88.2% 35.4 / ‑1.1 19P196P (CCG→CCA) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (9/6);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

AGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTC  >  NC_000913/4639491‑4639941
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
agaagTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAATCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACATAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGTCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCTCGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCa                                                                                                                                                                                                          <  1:369281/251‑1 (MQ=255)
       aCAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTttctt                                                                                                                                                                                                   >  1:101707/1‑251 (MQ=255)
             gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTcagcag                                                                                                                                                                                             <  1:101945/251‑1 (MQ=255)
                                    aTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGCGACTGAATTTTGc                                                                                                                                                                      >  1:68289/1‑251 (MQ=255)
                                      aaGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCt                                                                                                                                                                    >  1:301040/1‑251 (MQ=255)
                                        ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTgg                                                                                                                                                                  <  1:146153/251‑1 (MQ=255)
                                                                         gTTTTTTTTGACGGGGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGTTCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTTCCTGGGTTTTCACAGAAGGGAAGAATGGAACGGAACTc                                                                                                                                 >  1:333898/1‑251 (MQ=255)
                                                                                      gtggtAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGGGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGAATGAATTTTGCCTGGGTTTTCACAGTAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGc                                                                                                                    >  1:357967/1‑251 (MQ=255)
                                                                                        ggtAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGGCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCtt                                                                                                                  >  1:224012/1‑251 (MQ=255)
                                                                                           aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTa                                                                                                               <  1:98443/251‑1 (MQ=255)
                                                                                                        cTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTAGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGGCTGAATTTTGCCTGGGTTTTCACCGACGTGAAGCAAGGCGCGGAACTCGCTGCGCGGCAGCATGTACTGCTCGCCATc                                                                                                   >  1:57159/1‑250 (MQ=255)
                                                                                                                        gAAGCGATAACCTGTGCCGGGGATGGGGGCTATGATTTCCGGCGTATCCGGCTTAGATTCGAAAAGTTTACGAATCCGGCGTATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGCCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGTCGCGGAACTCGCGGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacgaac                                                                                  <  2:308999/251‑1 (MQ=255)
                                                                                                                                      cACCGTGTATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTc                                                                     <  1:339697/250‑1 (MQ=255)
                                                                                                                                                        tgatTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAAc                                                  >  1:118283/1‑251 (MQ=255)
                                                                                                                                                               cggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGGCACGTCTACAGTACGGTCGTGTGGTCTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGGTGTCCAGTTCCCATCCATTGAACTTGTAGc                                           <  1:114367/251‑1 (MQ=255)
                                                                                                                                                                                            tACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGCAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGAc               >  1:247634/1‑250 (MQ=255)
                                                                                                                                                                                                    cggcgaATCGTCACGTCTACAGTACGGTCGTGAGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGCGACTGAATTTTGCCTGGGTTTTCACAGAAGTGCAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGCACTTGTCGCTTTAAACGCTACGGCGTTCTTCGCTGACAGTAccc        >  2:19597/1‑249 (MQ=255)
                                                                                                                                                                                                     ggcggATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCCGTTCCCAACCATTGAACTTGTAGCTTTCAACGCAACGACGTTCTTCGCTGACAGTAccc        >  1:238159/1‑248 (MQ=255)
                                                                                                                                                                                                         gATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTAACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAATGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGACCGATCAACGAACGGCCGTTGATGTCCAGTTACCAACCATTGAACATGTAGCTTACAACGCTACGCCGTTCTTCGCTGACAGtcacaagatta  >  2:151945/1‑240 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
AGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTC  >  NC_000913/4639491‑4639941

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