Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,266,086	C→G	Q359E (CAG→GAG)	prfA →	peptide chain release factor RF‑1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,266,086	0	C	G	100.0%	109.8 / NA	31	Q359E (CAG→GAG)	prfA	peptide chain release factor RF‑1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (14/17);  total (14/17)

CGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGACCAACTGGCGGCGTTGTCCGAGCAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACT  >  NC_000913/1265957‑1266204                                                                                                                                  |                                                                                                                       ttgatccatgcgccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGaataat                                                                                                               <  1:748529‑M2/13‑1 (MQ=255)   gatccatgcgccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATgg                                                                                                             <  2:19677‑M2/15‑1 (MQ=255)         tgcgccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATc                                                                                                       <  1:592558‑M2/21‑1 (MQ=255)            gccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAac                                                                                                    <  1:1007253‑M2/24‑1 (MQ=255)             ccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAaca                                                                                                   <  1:169117‑M2/25‑1 (MQ=255)             ccaacgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAaca                                                                                                   <  2:226536‑M2/25‑1 (MQ=255)             ccaacgaaatcacgtcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAaca                                                                                                   <  1:324186‑M2/25‑1 (MQ=255)                    aatcacttcaacaggccgcccgtgacggtgataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGaa                                                                                                                   >  2:456780‑M2/108‑116 (MQ=255)                    aatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacttcattttctttttttacagggtgGAGGAGGaa                                                                                                                   <  1:456780‑M2/9‑1 (MQ=255)                          ttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGaa                                                                                      >  1:114748‑M2/102‑139 (MQ=255)                                  gccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATgg                                                                                                             <  1:187793‑M2/15‑1 (MQ=255)                                  gccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATgg                                                                                                             >  2:187793‑M2/94‑108 (MQ=255)                                             cggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAgg                                                                   <  2:200823‑M2/57‑1 (MQ=255)                                                                 atattctgcgcgacagtctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGt                                               >  2:488838‑M2/63‑139 (MQ=255)                                                                       tgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCt                                         >  2:899667‑M2/57‑139 (MQ=255)                                                                        gcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTg                                        >  2:501597‑M2/56‑139 (MQ=255)                                                                         cgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGa                                       >  1:626231‑M2/55‑139 (MQ=255)                                                                                 tctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCctgc                               <  1:488838‑M2/93‑1 (MQ=255)                                                                                 cctcgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCctgc                               <  1:918203‑M2/93‑1 (MQ=255)                                                                                    cgggctggagtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTgg                            <  1:501597‑M2/96‑1 (MQ=255)                                                                                             gtaacagtacatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTa                   <  2:228204‑M2/105‑1 (MQ=255)                                                                                                     acatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCgcggcg                                                >  1:883487‑M2/27‑102 (MQ=255)                                                                                                     acatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCgcggcg                                                <  2:883487‑M2/76‑1 (MQ=255)                                                                                                     acatcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCaga           >  2:919993‑M2/27‑139 (MQ=255)                                                                                                      catcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCagag          >  1:75643‑M2/26‑139 (MQ=255)                                                                                                      catcattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCagag          >  1:326844‑M2/26‑139 (MQ=255)                                                                                                         cattttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGc       <  2:114748‑M2/117‑1 (MQ=255)                                                                                                           ttttctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGt     >  2:12405‑M2/21‑139 (MQ=255)                                                                                                              tctttttttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACt  >  1:93901‑M2/18‑139 (MQ=255)                                                                                                                     ttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTa                   >  1:697132‑M2/11‑115 (MQ=255)                                                                                                                     ttacagggtgGAGGAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTa                   <  2:697132‑M2/105‑1 (MQ=255)                                                                                                                                  |                                                                                                                       CGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGACCAACTGGCGGCGTTGTCCGAGCAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACT  >  NC_000913/1265957‑1266204

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