Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,397,745 A→G 23.5% N236S (AAT→AGT)  PP_1214 → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,397,7450AG23.5% 25.1 / 6.2 17N236S (AAT→AGT) PP_1214transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (5/8);  new base G (2/2);  total (7/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.65e-01

ACGCGGAAATCGTCATGCAGCCGACCACCAGTGCCGAGCTCGACCAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCAT  >  NC_002947/1397625‑1397879
                                                                                                                        |                                                                                                                                      
aCGCGGAAATCGTCATGCAGCCGACCACCAGTGCCGAGCTCGACCAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCcgc                                                                                                                         >  1:176178/1‑135 (MQ=255)
aCGCGGAAATCGTCATGCAGCCGACCACCAGTGCCGAGCTCGACCAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTCGACGACGAGCAGAACGTCTATTCCAATGCCCAGATTTCcga                                                                                                                         >  1:5199/1‑136 (MQ=255)
                                       tCGACCAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCAtt                                                                                  >  2:242086/1‑136 (MQ=255)
                                          cccAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTtg                                                                               <  1:250181/135‑1 (MQ=255)
                                            cAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTtgtg                                                                             >  1:200466/1‑136 (MQ=255)
                                                   gTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAActc                                                                      <  2:215635/136‑1 (MQ=255)
                                                       cGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAg                                                                  <  2:212225/136‑1 (MQ=255)
                                                         aaaaGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCg                                                                >  2:379086/1‑136 (MQ=255)
                                                                 cTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAg                                                        <  1:19036/136‑1 (MQ=255)
                                                                 cTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAg                                                        <  1:168753/136‑1 (MQ=255)
                                                                                            cTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTaa                             <  2:200466/136‑1 (MQ=255)
                                                                                                    cgTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTc                     <  1:242086/136‑1 (MQ=255)
                                                                                                                  aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTactgca       <  2:207766/136‑6 (MQ=255)
                                                                                                                   ttCCAGTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCt                   >  1:351894/1‑123 (MQ=255)
                                                                                                                   ttCCAGTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCt                   >  1:349965/1‑123 (MQ=255)
                                                                                                                   ttCCAGTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCt                   <  2:349965/123‑1 (MQ=255)
                                                                                                                   ttCCAGTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCt                   <  2:351894/123‑1 (MQ=255)
                                                                                                                       aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCAt  <  2:372923/136‑1 (MQ=255)
                                                                                                                        |                                                                                                                                      
ACGCGGAAATCGTCATGCAGCCGACCACCAGTGCCGAGCTCGACCAGGACGGTGCCGAAAAGGTGCTCAAGCTGATCGACATGCTCGAAGACCTGGACGACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCAT  >  NC_002947/1397625‑1397879

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