Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,419,277 T→C 23.2% H68H (CAT→CAC PP_3905 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,419,2770TC23.2% 30.6 / 6.1 18H68H (CAT→CACPP_3905LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base C (2/2);  total (11/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 = 4.05e-01

CACCAGGCTGCATCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATG  >  NC_002947/4419143‑4419409
                                                                                                                                      |                                                                                                                                    
cACCAGGCTGCATCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCAta                                                                                                                                     >  1:490973/1‑135 (MQ=255)
cACCAGGCTGCATCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGACCGTCACGACCAATGCAGCCCACCCCCCAGAACAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCAta                                                                                                                                     >  1:46723/1‑135 (MQ=255)
         gCATCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCCCGACCAAATCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAg                                                                                                                            >  2:411078/1‑136 (MQ=255)
           aTCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGCCACGATCACTATAGGCGCATGGATAGGCCAACTCGGCCGAGGc                                                                                                                          >  2:40784/1‑136 (MQ=255)
                  gTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGgcgc                                                                                                                   >  1:347262/1‑136 (MQ=255)
                  gTGGCGTTCGCCCTCCCCAGGTCCGGGAGGTTCATGGCAGCGCGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGCCCCGATCCCTATCGGCGCATGGTTTGGCCATCTCGGCCGCGGCCTGgcgc                                                                                                                   >  1:207712/1‑136 (MQ=18)
                                              ggTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGt                                                                                       >  1:109329/1‑136 (MQ=255)
                                               gTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTa                                                                                      <  2:490973/136‑1 (MQ=255)
                                                   aTGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGccc                                                                                  >  2:326840/1‑136 (MQ=255)
                                                                           gCAGCCCACACCGAGATCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGcc                                                        >  1:277724/1‑136 (MQ=14)
                                                                           gCAGCCCACACCGAGATCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGcc                                                        >  1:512614/1‑136 (MQ=14)
                                                                                       cGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCg                                              <  1:411078/136‑1 (MQ=255)
                                                                                        gagaTCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGc                                            <  2:512614/136‑1 (MQ=18)
                                                                                        gagaTCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGc                                            <  2:277724/136‑1 (MQ=18)
                                                                                                                 aaTCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGaaa                    <  1:192484/135‑1 (MQ=255)
                                                                                                                          aTGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGt           >  2:27662/1‑136 (MQ=255)
                                                                                                                              aTAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACg       >  2:253519/1‑136 (MQ=255)
                                                                                                                              aTAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACg       >  1:565355/1‑136 (MQ=255)
                                                                                                                                  gCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAAt   <  2:347262/136‑1 (MQ=255)
                                                                                                                                   ccATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATg  <  1:383639/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CACCAGGCTGCATCGGTCGTGGCGTTCGCCCTCCCCTGGTCCGGGAGGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATG  >  NC_002947/4419143‑4419409

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