Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,313,853 T→A 57.5% I65F (ATC→TTC)  trpC ← fused indole‑3‑glycerolphosphate synthetase/N‑(5‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,313,8530TA57.5% 2.1 / 10.2 14I65F (ATC→TTC) trpCfused indole‑3‑glycerolphosphate synthetase/N‑(5‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base A (4/4);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

AAATCGGCTGCGGGGCGATTTGGCTGACGATGGGGAGGAAATTAAAGCTCCCCTGAAAATATTTCTCATCAGTCAGCACCGAAATTGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTTGCTGCTGTTTGCGGGCTT  >  CP009273/1313709‑1313993
                                                                                                                                                |                                                                                                                                            
aaaTCGGCGGCGGGGCGATCCGGCTGACGATGGGGAGGTTACTTCTACGACCGTTAAAATATTTCTCATCAGTCAGCACCGGAACTGCCGAAGCGGAAGGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGAACACg                                                                                                                                          <  1:84665/149‑1 (MQ=255)
   tCGGCTGCGGGGTGATTTGGCTGACGATGGGGACGAAATTAAAGCTCCCCTGAAAATATTTCTCATCAGTCAGCACCGAAATTGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGATCACGCCt                                                                                                                                       <  1:1156/149‑1 (MQ=255)
                                            aaGCGCCCCTGGAAATATTGCTCACCAGTCAGCACCGAAATGGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGAGCGAAATCATCACGGAACACGCCTTTCGACGGCGACGCTTTCTTGCACTCCAGAATaaa                                                                                                    <  2:130099/143‑1 (MQ=255)
                                                       aaaaTATTTCTCATCAGTCAGCACCGAAATTGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCt                                                                                   >  1:157690/1‑149 (MQ=255)
                                                           tatTTCTCATCAGTCAGCACCGAAATTGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCAcc                                                                                     <  2:172677/143‑1 (MQ=255)
                                                                                      gCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCgg                                                    <  1:341554/149‑1 (MQ=255)
                                                                                             cGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTCGCGCGTCATAAAAATGTCGCGTGCTCGGCTGAAcc                                             >  1:9057/1‑149 (MQ=255)
                                                                                             cGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCACCCCGAATCAGGACCATTGATGGGGACGTATTCTAGCATCCGAGAATACATGAACTTCGTTAACCCTGTAGCGCATCATAAAAATCTcgcggccccggc                                                   >  2:78902/1‑135 (MQ=255)
                                                                                                  tGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCAAGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAAc                                              >  2:76725/1‑143 (MQ=255)
                                                                                                                  cAATGCGTGCTGGATCGAAATCATCACGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGGGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTg                              <  2:952/143‑1 (MQ=255)
                                                                                                                             ggATCGAAATCATCACGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGAGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTtg                   >  2:8918/1‑143 (MQ=255)
                                                                                                                                      tcatcaCGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAATAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTTGCTGCTGTTTGCGGGc    <  1:223684/149‑1 (MQ=255)
                                                                                                                                        atcaCGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTTGCTGCTGTTTGCGGGCtt  >  1:50732/1‑149 (MQ=255)
                                                                                                                                         tcaCGGAACACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTTGCTGCTGTTTGCg       <  2:330103/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AAATCGGCTGCGGGGCGATTTGGCTGACGATGGGGAGGAAATTAAAGCTCCCCTGAAAATATTTCTCATCAGTCAGCACCGAAATTGCCGAAGCGTAATGTTTATAAATGGCGGCAATGCGTGCTGGATCGAAATCATCACGGATCACGCCTTTTGACGGCGACGCTTTCTTGCACTCCAGAATAAACGCCGTGCGCGCACCCTGTAGCGCATCATAAAAATGTCGCGTGCTCGGCTGAACCTCATTCTGAAAACTGGCCAGCGGTTGCTGCTGTTTGCGGGCTT  >  CP009273/1313709‑1313993

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