Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,339 C→T 37.8% T34T (ACC→ACT talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,3390CT37.8% 3.9 / 13.6 16T34T (ACC→ACTtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base T (4/2);  total (11/5)
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.12e-01

AATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAA  >  CP009273/8229‑8476
                                                                                                              |                                                                                                                                         
aaTACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGaa                                                                                                    >  2:41101/1‑150 (MQ=255)
     tATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCg                                                                                               >  1:149455/1‑150 (MQ=255)
        aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                                                                                              >  1:183315/2‑148 (MQ=255)
          tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                                                                                          <  2:59430/150‑2 (MQ=255)
            aCGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttg                                                                                        >  1:97702/1‑150 (MQ=255)
            aCGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGtt                                                                                         >  1:166783/1‑149 (MQ=255)
              ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                                                                                      <  1:155216/150‑1 (MQ=255)
                     ttGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTgatgat                                                                               <  2:166783/150‑1 (MQ=255)
                              cTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGcc                                                                      >  1:24705/1‑150 (MQ=255)
                                 cGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCt                                                                     <  1:41101/148‑1 (MQ=255)
                                                                        gcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctc                            >  2:42749/1‑147 (MQ=255)
                                                                            cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                        >  1:46430/1‑143 (MQ=255)
                                                                             aaTGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGAAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAcgc                       >  2:126440/1‑150 (MQ=255)
                                                                                               gCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAg      >  1:145049/1‑149 (MQ=255)
                                                                                                cAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGt     >  1:102498/1‑149 (MQ=255)
                                                                                                  ggATGCCACAACCAACCCTTCTCTCATGCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa  <  2:97702/150‑1 (MQ=255)
                                                                                                              |                                                                                                                                         
AATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAA  >  CP009273/8229‑8476

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