Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,264,970	G→A	*419* (TAG→TAA)	hemA →	glutamyl tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,264,970	0	G	A	100.0%	70.6 / NA	23	*419* (TAG→TAA)	hemA	glutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/10);  total (13/10)

TCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264850‑1265103                                                                                                                         |                                                                                                                                      tCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCttt                                                                                                                     >  2:164034/1‑139 (MQ=255)               tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGga                                                                                                       >  1:219209‑M1/1‑137 (MQ=255)                                 aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatat                                                                                     >  1:464827‑M1/1‑119 (MQ=255)                                         aCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactg                                                                             >  2:202848‑M1/1‑111 (MQ=255)                                                                       ggggATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggc                                               <  1:202848‑M1/139‑59 (MQ=255)                                                                               cgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     >  1:255142‑M1/1‑73 (MQ=255)                                                                               cgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgagg                                                                                                     <  2:255142‑M1/77‑5 (MQ=255)                                                                                   cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataa                                                                                             >  2:636727‑M1/1‑69 (MQ=255)                                                                                   cgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataa                                                                                             <  1:636727‑M1/81‑13 (MQ=255)                                                                                              tCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggc                                               >  2:97187‑M1/1‑58 (MQ=255)                                                                                              tCTGCGCGACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggc                                               <  1:97187‑M1/116‑59 (MQ=255)                                                                                                  cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctga                    >  2:521330‑M1/1‑54 (MQ=255)                                                                                                  cgcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctga                    >  1:425053‑M1/1‑54 (MQ=255)                                                                                                   gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaata                                                                                      <  2:356841‑M1/72‑20 (MQ=255)                                                                                                   gcgACAGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaata                                                                                      >  1:356841‑M1/1‑53 (MQ=255)                                                                                                        aGCCTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcct              <  2:425053‑M1/139‑92 (MQ=255)                                                                                                          ccTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              >  1:617012/1‑41 (MQ=255)                                                                                                          ccTCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACa                                                                                                              <  2:617012/41‑1 (MQ=255)                                                                                                            tCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                             <  1:757435/40‑1 (MQ=255)                                                                                                            tCGGGCTGGAGTAACAGTACATCATTTTCTTTTTTTACAg                                                                                                             >  2:757435/1‑40 (MQ=255)                                                                                                                    gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatcaacactggttacgtgaagcaataagccaacttcaggcgagcgaaagcccgcggcgtgatgctgaaatcctgctggagcatgt  <  1:521330‑M1/139‑104 (MQ=255)                                                                                                                    gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                    <  1:502547‑M1/57‑22 (MQ=255)                                                                                                                    gAGTAACAGTACATCATTTTCTTTTTTTACAGGGTGgaggaggaataatggaatatc                                                                                    >  2:502547‑M1/1‑36 (MQ=255)                                                                                                                         |                                                                                                                                      TCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGA  >  NC_000913/1264850‑1265103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑