Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,437,890	C→T	Q62* (CAG→TAG)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,437,890	0	C	T	100.0%	52.9 / NA	18	Q62* (CAG→TAG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/8);  total (10/8)

TTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTAC  >  NC_000913/4437664‑4438136                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       tttAATGGAGACATAAAGTAACAAAATCAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAAAAATTTTATTGGTGCGGAGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAgg                                                                                                                                                                                                                                >  2:234288/1‑251 (MQ=255)  ttAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAACAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGt                                                                                                                                                                                                                               >  2:57205/1‑251 (MQ=255)                         aataaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATATTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTg                                                                                                                                                                                                        >  2:198421/1‑251 (MQ=255)                            aaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTc                                                                                                                                                                                                     <  1:65316/251‑1 (MQ=255)                               ataaGGTTGTCATATGTTACCCAGGATCAGACACAAGAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCATACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTAGTTGCGAAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTAATCCAAACATGAAACTGTTCACg                                                                                                                                                                                                  >  2:142739/1‑251 (MQ=255)                                                 aCCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTa                                                                                                                                                                                <  1:39443/251‑1 (MQ=255)                                                            gACATAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAATCTTCGTATACAAAACCATATTAAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGTATTTAGAATTATGATAGCTTATTAACTTTGAGGGAAGTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTAGCAGGTGCTATAATCAAAAAATGAAACTGTTCACGTCATGGATAAGTTTCTTATCAGTTTCAGc                                                                                                                                                                     <  2:149981/251‑1 (MQ=255)                                                                          aTTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACgg                                                                                                                                                       <  1:25248/251‑1 (MQ=255)                                                                                    tGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATAATGCATTTAAGCAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAatat                                                                                                                                             >  1:155850/1‑251 (MQ=255)                                                                                        gTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACaa                                                                                                                                         <  1:213432/251‑1 (MQ=255)                                                                                          tATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGTAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAATTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAAt                                                                                                                                        <  1:156829/250‑1 (MQ=255)                                                                                                                                                             tAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACtcattc                                                                     >  2:174871/1‑250 (MQ=255)                                                                                                                                                                          tGGATTTCGAATTAGGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGAAACTTAGATTGATAAAAGTGGTTAGCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCaaaa                                                        <  2:185041/250‑1 (MQ=255)                                                                                                                                                                            gATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAAcc                                                     >  2:71622/1‑251 (MQ=255)                                                                                                                                                                               ttCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACATAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTg                                                  >  1:93958/1‑251 (MQ=255)                                                                                                                                                                                    aTTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGGTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAATCTGTGCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAAttt                                             <  1:37678/251‑1 (MQ=255)                                                                                                                                                                                                        ttGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                         >  1:14468/1‑251 (MQ=255)                                                                                                                                                                                                                               aaCTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTAc  >  1:88580/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       TTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTAC  >  NC_000913/4437664‑4438136

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