Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,256,708 G→T 21.0% V104V (GTC→GTA tdcD ← propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,256,7080GT21.0% 36.4 / 3.2 19V104V (GTC→GTAtdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base T (2/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

GCCACCTGAGTTACGCCCGGAAATAATTGCTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGA  >  CP009273/3256601‑3256962
                                                                                                           |                                                                                                                                                                                                                                                              
gCCACCTGAGTTACGCCCGGAAATAATTACTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCa                                                                                                                                                       >  1:369891/1‑213 (MQ=255)
 ccACCTGAGTTACGCCCGGAAATAATTGCTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATTACTTCATCg                                                                                                                                                                                                                                                       >  2:46314/1‑116 (MQ=255)
    ccTGAGTTACGCCCGGAAATAATTGCTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATg                                                                                                                                                   <  2:369878/213‑1 (MQ=255)
     cTGAGTTACGCCCGGAAATAATTGCTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATTACTTCATCGGtaa                                                                                                                                                                                                                                                   <  1:46316/116‑1 (MQ=255)
                                cgcCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAg                                                                                                               >  1:71140/1‑221 (MQ=255)
                                    gATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAg                                                                                                           <  2:71138/221‑1 (MQ=255)
                                                cTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGAc                                                                                                                                                                                                                                                              >  2:218419/1‑62 (MQ=255)
                                                   aaaTTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATTACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTGTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAg                                                                                                               >  1:94678/1‑202 (MQ=255)
                                                   aaaTTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGGAATAATGGCGGACTAGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTTTGGTGAGCAAGTGGTGCTGGCTCTCCCCCATTAACGGAAAAGAATGAATTTCCCGAGTTAATCCCTTCGGccacatct                                        >  2:415333/1‑266 (MQ=255)
                                                    aaTTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCa                                                                                                                                                                                                                                                          <  1:218426/62‑1 (MQ=255)
                                                      ttGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACAt                                   <  1:61014/275‑1 (MQ=255)
                                                       tGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATTACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAg                                                                                                           <  2:94676/202‑1 (MQ=255)
                                                           gTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATa                              >  1:75324/1‑275 (MQ=255)
                                                                         gTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCg                                                                                                                                                                                                                                           >  1:213690/1‑56 (MQ=255)
                                                                         gTGCCAGTGGAGAAACGCGACGGATAATATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCgg                                                                                                                                                                                                       >  1:324132/1‑92 (MQ=255)
                                                                             cAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACt                                                                                                                                                                                                                                       <  2:213683/56‑1 (MQ=255)
                                                                             cAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGc                                                                                                                                                                                                   <  2:324122/92‑1 (MQ=255)
                                                                                       aCGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTCTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGTATAATAATGCATTTTACGAGTTAATACCGTCGCAAATACATGACATTAAAACTTCACAGTCGCTGGCATCGAgaccggc  >  2:447930/1‑269 (MQ=255)
                                                                                                          tGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCa                                                                                                                                                                                                 >  1:1072/1‑65 (MQ=255)
                                                                                                           |                                                                                                                                                                                                                                                              
GCCACCTGAGTTACGCCCGGAAATAATTGCTGCGCCGATTCAATACCACTTAAATTGGCGTAATTATGCAGGGGTGCCAGTGGAGAAACGCGACGGATATTATCAATGACTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGA  >  CP009273/3256601‑3256962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: