Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496 A→G E232E (GAA→GAG thrB → homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,4960AG100.0% 26.7 / NA 9E232E (GAA→GAGthrBhomoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/4);  total (5/4)

GCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCC  >  NC_000913/3377‑3618
                                                                                                                       |                                                                                                                          
gcgcAGTATCGCCGCCAGCATTCCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAAc                                                                                                               <  1:132337/133‑1 (MQ=255)
 cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACg                                                                                                              >  2:380113/1‑133 (MQ=255)
              ccAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGc                                                                                                 >  1:62328/1‑133 (MQ=255)
                  gATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGCTATCGCTGCGCCCTACAGTGAACGGATACTCCCAGGCTTcc                                                                                             >  2:171535/1‑133 (MQ=255)
                         ttgcGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCcggcaggc                                                                                      >  1:502428/1‑133 (MQ=255)
                                                           cGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTagcg                                                    <  2:502428/133‑1 (MQ=255)
                                                             ccTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTagcgag                                                  <  1:380113/133‑1 (MQ=255)
                                                                                       gCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCtccggctccggc                               <  1:52431/126‑1 (MQ=255)
                                                                                       gCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCtccggctccggc                               >  2:52431/1‑126 (MQ=255)
                                                                                                             ttATCGCTGAGCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGcc  >  2:483642/1‑133 (MQ=255)
                                                                                                                       |                                                                                                                          
GCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCC  >  NC_000913/3377‑3618

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