Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,654,922 G→A 18.2% A190T (GCC→ACC)  cysB → transcriptional regulator CysB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,654,9220GA18.2% 41.5 / 4.8 22A190T (GCC→ACC) cysBtranscriptional regulator CysB
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base A (2/2);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGG  >  NC_002947/2654787‑2655069
                                                                                                                                       |                                                                                                                                                   
ttCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATc                                                                                                                                       <  2:254788/150‑1 (MQ=255)
         gccaCCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTAc                                                                                                                              <  1:306870/150‑1 (MQ=255)
                                          aTCATGATGCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGAt                                                                                             >  2:57760/1‑150 (MQ=255)
                                              tgatgCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCg                                                                                                                       <  2:297068/120‑1 (MQ=255)
                                              tgatgCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCg                                                                                                                       >  1:297068/1‑120 (MQ=255)
                                                  gCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCttt                                                                                     <  2:121430/150‑1 (MQ=255)
                                                           cAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCg                                                                            >  1:420773/1‑150 (MQ=255)
                                                               tGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTg                                                                                                >  1:470997/1‑126 (MQ=255)
                                                               tGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTg                                                                                                <  2:470997/126‑1 (MQ=255)
                                                                     cgttgcgttgTAGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACg                                                                  >  2:152894/1‑150 (MQ=255)
                                                                       ttgcgttgTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGcc                                                                >  1:192782/1‑150 (MQ=255)
                                                                                ggtgCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGt                                                       <  1:29856/150‑1 (MQ=255)
                                                                                          gggCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACcgcg                                             >  1:348174/1‑150 (MQ=255)
                                                                                            gCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGc                                           >  1:485315/1‑150 (MQ=255)
                                                                                                                 gAAACTGACCCTGGAAGCCGTCACCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGAc                      >  2:223891/1‑150 (MQ=255)
                                                                                                                 gAAACTGACCCTGGAAGCCGTCACCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCGCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGAc                      >  2:224216/1‑150 (MQ=255)
                                                                                                                     cTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTAt                  >  1:486960/1‑150 (MQ=255)
                                                                                                                              gAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTg         >  2:379479/1‑150 (MQ=255)
                                                                                                                               aaGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTgg        <  2:351232/150‑1 (MQ=255)
                                                                                                                                  ccGTCACCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGt     <  1:224216/150‑1 (MQ=255)
                                                                                                                                  ccGTCACCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGt     <  1:223891/150‑1 (MQ=255)
                                                                                                                                     tCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTgg  <  2:486960/150‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
TTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGCTACAAGTGGAACCGTTGCGTTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGG  >  NC_002947/2654787‑2655069

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