Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,322,325 G→A L208F (CTC→TTC)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3250GA95.0% 46.7 / ‑3.3 20L208F (CTC→TTC) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/8);  total (12/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 = 5.29e-01

CACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGT  >  NC_000913/1322098‑1322539
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                      
cACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCg                                                                                                                                                                                                 >  1:354383/1‑251 (MQ=255)
                 gCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAAc                                                                                                                                                                                         <  1:273060/242‑1 (MQ=255)
                                                            cTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCCTCGCTCTGATTACCTTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTAACGTAGTCCGTTCAGGCGAGCAGTGAAACgttgt                                                                                                                                                                                                                 >  1:500428/1‑175 (MQ=255)
                                                                aaaaTTTCTCCAGCGCGAATCGCGTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATGCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTGTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTgn                                                                                                                                                                             <  1:316208/207‑2 (MQ=255)
                                                                                     gCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCtt                                                                                                                                                                 >  1:216110/1‑198 (MQ=255)
                                                                                                                  gAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGCCACGGAAACACCTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTATTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCtttttt                                                                                                                                                             >  1:270194/1‑173 (MQ=255)
                                                                                                                          naTCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCactc                                                                                                                                               >  1:96399/2‑176 (MQ=255)
                                                                                                                                                          ggCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAATCAGGCTGGc                                                                                                                                                                                 <  1:411797/113‑1 (MQ=255)
                                                                                                                                                          ggCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACa                                                                                                                                                                                        >  1:14543/1‑106 (MQ=255)
                                                                                                                                                              cGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTGGCTGACGTAGTTCGTTCGGGCGGGAAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGATGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCg                                                                               >  1:265902/1‑207 (MQ=255)
                                                                                                                                                                       cTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAAGCGCGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACACGGTCATCAGAGAACAGGCCGCCGAAGAAc                          >  1:250212/1‑251 (MQ=255)
                                                                                                                                                                            aGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACANNCTGGCCNNa                                                                                                                                                                             >  2:26625/1‑99 (MQ=255)
                                                                                                                                                                                         cTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAggn                                                                                                        <  1:55826/155‑2 (MQ=255)
                                                                                                                                                                                           tcgttCAGTTGCGGACGTAGTTCGTTCAGGCGAGCAGTGAAATGTTTTTTTTCTGCTTCATTCGGAACAAACGGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGTTCAATCACCATCAGCGTTTCAGCGAGAAAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACATTTGGGGAAAATTTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcc        <  1:26625/245‑1 (MQ=255)
                                                                                                                                                                                                cAGTTGCTGACGTAGTTCGTTCAGGCGAGAAGTGAAACGTTGTTTTTCTCATTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGCAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGt  >  1:376277/1‑250 (MQ=255)
                                                                                                                                                                                                      cTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATANGGGNGCTTNTTNTCTGATGGTCAATCANCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATNNTNNNCNGACAGTTg                                                                                 <  2:356119/165‑1 (MQ=255)
                                                                                                                                                                                                      cTGACGTAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAAAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                 >  1:356119/1‑165 (MQ=255)
                                                                                                                                                                                                            tAGTTCGTTCAGGCGAGCAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTGTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                                                                          <  1:451805/150‑1 (MQ=255)
                                                                                                                                                                                                                tcgttcAGGCGAACAGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTNNATCACCNNc                                                                                                                                         >  2:27225/1‑99 (MQ=255)
                                                                                                                                                                                                                              aGTGAAACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGc                                                                                                                                                                                     <  1:439687/41‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                      
CACGTAATAGGCCGCCAGCGGTGACGGGCAGGGCAGAGAGAAACGGCGAGATGGCACCACCTGGAAAATTTCTCCAGCGCGAATCGCTTTTTGCAACAAACGCACTACGCCACCGAACTCTTCATCGCTCTGATTACATTCACAACGCATATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGT  >  NC_000913/1322098‑1322539

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