Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,979,132 C→T 20.0% I50I (ATC→ATT serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,979,1320CT20.0% 40.1 / 5.0 20I50I (ATC→ATTserC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (6/10);  new base T (2/2);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

TTAACTTCTGCGCAGGCCCTGCCGCGCTTCCTGACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCCAGCCAGCAGTTCGCTGAAATCCCGCTGAACCTGCTGCCGG  >  NC_002947/1978999‑1979253
                                                                                                                                     |                                                                                                                         
ttAACTTCTGCGCAGGCCCTGCCGCGCTTCCTGACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGc                                                                                                                         >  1:336368/1‑136 (MQ=255)
          cgcAGGCCCTGCCGCGCTTCCTGACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGcc                                                                                                               >  1:19273/1‑136 (MQ=255)
               gCCCTGCCGCGCTTCCTGACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCa                                                                                                          <  2:336368/136‑1 (MQ=255)
                       gcgcTTCCTGACGCTGTCCAGCAGCGCGCACAGGCCGAGATGCAGGACAGGCGTGGCAAGGGGTTGTCGGTGATGGAAAATAGCCAACCCAGCGACGATTACCTGGCCAATGCCCAAAACGCCCAGCACGACCAGc                                                                                                  >  2:265336/1‑136 (MQ=255)
                       gcgcTTCCTGACCCTGTACTGCAGCGCGCACAGGCCCAGAAGCTGGACCGGCGTGGCAAGGGGTTGTCGGTGATGGAAAAGATCCAACCCAGCGACGATTAACTGGCCAATGACCAAAAGGCCGAGCAGGACCAGc                                                                                                  >  2:262750/1‑136 (MQ=255)
                                gACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACctgc                                                                                         >  1:210303/1‑136 (MQ=255)
                                     tGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTgtcc                                                                                    <  2:527153/136‑1 (MQ=255)
                                      gTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTgtccg                                                                                   >  2:381030/1‑136 (MQ=255)
                                                  gcACAGGCCGAGATGCTGGACGGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATTGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACt                                                                       <  1:265336/136‑1 (MQ=255)
                                                  gcACAGGCCGAGATGCTGGACGGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATTGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACt                                                                       <  1:262750/136‑1 (MQ=255)
                                                        gCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAgg                                                                 <  2:210303/136‑1 (MQ=255)
                                                           gagaTGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGc                                                              <  1:419162/136‑1 (MQ=255)
                                                           gagaTGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGc                                                              <  1:456450/136‑1 (MQ=255)
                                                           gagaTGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGc                                                              <  1:2648/136‑1 (MQ=255)
                                                                          ggTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAggcg                                               <  2:128882/135‑1 (MQ=255)
                                                                          cGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAggcg                                               <  2:484948/136‑1 (MQ=255)
                                                                                      ttGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCcagccagc                                   >  2:10705/1‑136 (MQ=255)
                                                                                                           cATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCCAGCCAGCAGTTCGCTGAAATCCCGCTGa              <  2:307198/136‑1 (MQ=255)
                                                                                                                  gCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCATCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCCAGCCAGCAGTTCGCTGAAATCCCGCTGAACctgct       <  2:19273/136‑1 (MQ=255)
                                                                                                                       gaTTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCCAGCCAGCAGTTCGCTGAAATCCCGCTGAACCTGCTGCCgg  >  1:37594/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                         
TTAACTTCTGCGCAGGCCCTGCCGCGCTTCCTGACGCTGTCCTGCAGCGCGCACAGGCCGAGATGCTGGACTGGCGTGGCAAGGGGTTGTCGGTGATGGAAATGAGCCATCGCAGCGACGATTACGTGGCCATCGCCGAAAAGGCCGAGCAGGACCTGCGTGACCTGCTGTCCGTCCCCTCCAACTACAAGGTGCTGTTCCTGCAAGGCGGCGCCAGCCAGCAGTTCGCTGAAATCCCGCTGAACCTGCTGCCGG  >  NC_002947/1978999‑1979253

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: