Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,574,250 G→A 18.2% G419D (GGC→GAC)  tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,574,2500GA18.2% 39.4 / 2.4 22G419D (GGC→GAC) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (9/9);  new base A (2/2);  total (11/11)
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.34e-01

TTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACGACTCAATTGGCCTGG  >  CP009273/2574114‑2574390
                                                                                                                                         |                                                                                                                                            
ttGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGGGCGTGAATTTGGCATGACCGCTCTCGCCAGCGGCATcgctc                                                                                                                                      >  2:151658/1‑143 (MQ=255)
          tttgttcggATCTTGCTCCCATCCACCTGACCATCTCGATAGAGATATATTCGCTGACGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcggcgg                                                                                                                            <  2:80835/139‑1 (MQ=255)
          gTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCACCGCGCACACcggcgg                                                                                                                            >  1:76486/1‑145 (MQ=255)
                  gATCTGGCTCCCAGCAACCTGACCATCTAGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCCCTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGc                                                                                                                    <  1:237131/145‑1 (MQ=255)
                     cTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACAGCTATCGCACACGGCATCGCGTACCACTGCGGTTTTGTGCCGt                                                                                                                 <  1:354945/145‑1 (MQ=255)
                              aGCAATCTGATCATTTGGAAATGTTATGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGa                                                                                                        <  2:131342/145‑1 (MQ=255)
                                  aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTCCGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGACATCGCGCACCACGGCGGCTTTGTGCCGAATACCGCGACGtt                                                                                                    >  2:60046/1‑145 (MQ=255)
                                  aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGACATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGtt                                                                                                    >  2:268964/1‑145 (MQ=255)
                                                           tttCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGGTCCTGATGGTTGTTGAATACGACCg                                                                            >  1:352150/1‑145 (MQ=255)
                                                              cGCTGAAGGAAGATCCAGCGGGCAACTCTATTCACTACGGGGTGCGTGCACTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACTTTCCTGATGTTTGTTGAATACGCCCGTaa                                                                         <  1:233154/145‑1 (MQ=255)
                                                                          aTCCAGCGGGCAACTACAATCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGAt                                                             >  2:311603/1‑145 (MQ=255)
                                                                                cGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATa                                                                                  <  1:17404/118‑1 (MQ=255)
                                                                                cGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATa                                                                                  >  2:17404/1‑118 (MQ=255)
                                                                                                     gtgtgtgtGAATTTGGCATGACCGCTATGGCCAACGGCCTCGCGCACCACGGCGGCTTTGTGCCGTATCCCGCGACGTTTCTGACGTTTGTTGCATCCGCCCGTGACACCGCGCGGCTGCCGCCATTGTGGAACGCGCCGCAGAt                                  >  1:40287/3‑145 (MQ=255)
                                                                                                         gCGTGGATTTTGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCGGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATg                              <  1:60180/145‑1 (MQ=255)
                                                                                                                tttGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTtata                       >  1:338757/1‑145 (MQ=255)
                                                                                                                     cATGCCTGCTATCGGCAACGGCATCGCGCACCACTGCGGCTTTGTGCCGTATACCGCGACGTTCCTGCTGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCAc                  <  2:239985/145‑1 (MQ=255)
                                                                                                                      aTGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACg                 >  2:338715/1‑145 (MQ=255)
                                                                                                                      aTGACCGCTATCGCCAACGACCTCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACg                 <  1:60046/145‑1 (MQ=255)
                                                                                                                      aTGACCGCTATCGCCAACGACATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACg                 <  1:268964/145‑1 (MQ=255)
                                                                                                                                   tcAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACGACTCAATTGGCCt    >  2:138065/2‑145 (MQ=255)
                                                                                                                                     aaCGGCATCGCGCACCACGGCGGCTTGGTGCCGTATACCGCGACGTCCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACGACTCAAGTGGCCTgg  <  1:336578/145‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TTGCTCGGCGGTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTTGAATACGCCCGTAACGCCGCGCGGATGGCGGCACTGATGAAAGCGCGGCAGATTATGGTTTATACCCACGACTCAATTGGCCTGG  >  CP009273/2574114‑2574390

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