Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,342 C→A 21.0% G87V (GGG→GTG)  rutR ← transcriptional regulator, TetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,3420CA21.0% 38.9 / 2.6 19G87V (GGG→GTG) rutRtranscriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base C (10/5);  new base A (2/2);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

TGGGTAGACGCCCAGATCATATACAGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTGGTCAGCACCACACGCGGGTCGTCGAAGTTC  >  NC_002947/4557210‑4557467
                                                                                                                                    |                                                                                                                             
tGGGTAGACGCCCAGATCATATACAGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAAcccccc                                                                                                                            >  1:285548/1‑135 (MQ=255)
               aTCATATACAGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCgg                                                                                                             >  1:38348/1‑135 (MQ=255)
               aTCATATACAGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGt                                                                                                             >  1:427072/1‑136 (MQ=255)
                        aGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCAtct                                                                                                    <  1:358420/136‑1 (MQ=255)
                        aGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCAtct                                                                                                    <  1:441528/136‑1 (MQ=255)
                                cGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGcc                                                                                            >  1:2642/1‑136 (MQ=255)
                                cGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGcc                                                                                            <  2:436467/136‑1 (MQ=255)
                                cGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGcc                                                                                            >  1:117458/1‑136 (MQ=255)
                                cGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCACCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGcc                                                                                            >  1:471408/1‑136 (MQ=255)
                                cGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGCGCACCAGGCTTTCGTCGAGGTTAACCACCAGCACCGGCGCCCCGTGCAGCATCTCTTCggac                                                                                            >  1:335853/1‑134 (MQ=255)
                                     gCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGAt                                                                                       >  1:44083/1‑136 (MQ=255)
                                     gCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGAt                                                                                       >  1:278865/1‑136 (MQ=255)
                                            cGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAg                                                                                >  1:4980/1‑136 (MQ=255)
                                                               tCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGcgcg                                                             >  2:290565/1‑136 (MQ=255)
                                                                                     cttcttTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGAt                                       >  2:490807/1‑136 (MQ=255)
                                                                                      ttcttTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATg                                      <  2:278865/136‑1 (MQ=255)
                                                                                             aTCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCACCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCTGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTgg                               <  2:335853/136‑1 (MQ=255)
                                                                                             aTCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCACCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTgg                               <  2:471408/136‑1 (MQ=255)
                                                                                                            aCCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTGGTCAGCACCACAcgcg                >  2:494570/1‑136 (MQ=255)
                                                                                                                          gAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTGGTCAGCACCACACGCGGGTCGTCGAAGTtc  <  2:44083/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
TGGGTAGACGCCCAGATCATATACAGCAACGCCGAAGGCTCTACCGGCAGAATACGCCCCTCTTCCACCCAGCAGCGGATCTTGGCTTCTTTCATCTTCGCCCAGGGCACCAGGCTTTCGTCGAGGTTAACCCCCAGCACCGGCGCCCCGTGCAGCATCTCTTCGGCCCAGATCTTCGAGCCCAGCGGCCGCGAGCGCGAATGGCCCATCTTCGCCCGGATGTAGCTGGTCAGCACCACACGCGGGTCGTCGAAGTTC  >  NC_002947/4557210‑4557467

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