Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,233,609	C→T	Q4* (CAG→TAG)	yeiS →	DUF2542 domain‑containing protein YeiS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,233,609	0	C	T	100.0%	52.7 / NA	16	Q4* (CAG→TAG)	yeiS	DUF2542 domain‑containing protein YeiS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/6);  total (10/6)

GTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTATCTTATTTTCTATCGTTAATAT  >  NC_000913/2233419‑2233842                                                                                                                                                                                               |                                                                                                                                                                                                                                          gTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCAt                                                                                                                                                                               <  1:227691/251‑1 (MQ=255)                                    cTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACgg                                                                                                                                            >  1:153739/1‑251 (MQ=255)                                        atatTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGAGGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTa                                                                                                                                        <  2:108417/251‑1 (MQ=255)                                            ttttAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCAAACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTaaaa                                                                                                                                     >  1:29087/1‑250 (MQ=255)                                                aaaCGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCa                                                                                                                                <  2:81120/251‑1 (MQ=255)                                                                  ttAGGGCCGCTGGTCCCTATTCCGGCTAGGCACCAGGTACCGGATGATGCGCAGGGTTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGATGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGt                                                                                                               <  1:215828/250‑1 (MQ=255)                                                                       gCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGCACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCg                                                                                                         >  1:129986/1‑251 (MQ=255)                                                                               taccTATTCAGGCTATAACCCAGGTATCGGCAGCAGCGCAGGGGTATTCCGCCGTCACACCCGCACAGCTACTTGAATTACATAAGACACAAGGAATGTCGTTATGGATGTATAGCAGTTTTTTGTCGATGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGTAACATGGAAATGCTGGCGACCACGCGCGTTTGATAAACGGGTTAAAACTGCAATGGACCCGACGTATGTATGGCGAGCaaaa                                                                                                  >  2:209774/3‑250 (MQ=255)                                                                                cccTATTCCGGCTATGCTCCAGGTACCTGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACGTGGATTTCAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTGTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGGGCGATTGATAAGCGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCaaaaaa                                                                                                <  1:24215/251‑1 (MQ=255)                                                                                cccTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCaaaaa                                                                                                 >  1:222710/1‑250 (MQ=255)                                                                                           cTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGCTAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGActct                                                                                     >  2:60906/1‑251 (MQ=255)                                                                                                                                                cAGTTACTTGAATTACAACAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAACCGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATg                                >  2:65383/1‑251 (MQ=255)                                                                                                                                                           aTTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGCTTGTTTATCt                     >  2:144286/1‑251 (MQ=255)                                                                                                                                                               cAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTATCTTAtt                 >  1:4776/1‑251 (MQ=255)                                                                                                                                                                aaaaGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTATCTTAttt                >  1:230629/1‑251 (MQ=255)                                                                                                                                                                              aaGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATAGCGGCTTCGGAATTTTATCTATCGGCATGATTGTTGATCGTATTTTCTAGCTTTAatat  <  2:129986/251‑1 (MQ=255)                                                                                                                                                                                               |                                                                                                                                                                                                                                          GTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTATCTTATTTTCTATCGTTAATAT  >  NC_000913/2233419‑2233842

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