Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	363,292	G→T	S1005* (TCG→TAG)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	363,292	0	G	T	100.0%	61.2 / NA	18	S1005* (TCG→TAG)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (12/6);  total (12/6)

TTTTGCTGATATGGTTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTA  >  NC_000913/363053‑363525                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          ttttGCTGATATGGTTGATGTCATGTAGCCAACTCGGGAAAACCGGGAAGTAGGCTCCCATGATAAAAAAGCAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGAGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGACTTCCGCCTATACTGACggg                                                                                                                                                                                                                                >  2:104636/1‑251 (MQ=255)      ctgctatGGTTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAAACGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCAGCCTATACTGAAGGGCTcc                                                                                                                                                                                                                            >  2:225935/5‑250 (MQ=255)                                               gAAGTAGGCTCCCATGATAAAAAAATACAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTATCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCACTCCCCATATGGAAACCGTCGctat                                                                                                                                                                                  >  2:22607/1‑251 (MQ=255)                                                                 aaaaaaGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCcttctt                                                                                                                                                                <  1:259004/251‑1 (MQ=255)                                                                        tAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCCGACATGGCCTGCACGGTTATTATTATTTTAGCCACAAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTACGCCTATACTGCCGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTg                                                                                                                                                         >  1:189091/1‑251 (MQ=255)                                                                                                            gtgtTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCTATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTg                                                                                                                     <  2:175365/251‑1 (MQ=255)                                                                                                             tgtTTTTTCAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGAACGGTTATTATTATTTTTGACACCAGAACAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCCCCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCCGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGt                                                                                                                    >  2:161021/1‑251 (MQ=255)                                                                                                                tttttAAATAGTACATAATGGATATCATTACGCGAAATACGGGCAGTCATGCCCTGCCTGGTTATTATAATTTTTGACAACAGCACAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATCCCGCCTATACTGTCGCGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCAATGCGCCTTCTTCCGCGTGCAGCAGATGGCGCTGGCTGGTTTCCATCAGTTGCTGTTgt                                                                                                                 >  2:26356/1‑250 (MQ=255)                                                                                                                         aGTACATCATGTATTTCCTTACGTGAAATACTGGCATACCTTGCCTGCCCGGTTATTATTATTTTTGACACCAGACCCACTGGTAATGGAAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGCGTCGTCGCCCCCAATCCCCATATGGCAACCGCCGATATTCAGCCATGCGCCTTCTTCCGCGTGCAGAAGATGACACTAGCTGGTTTCCAACAGTTGCTGTTGACTGTAGCgg                                                                                                        >  2:2419/1‑251 (MQ=255)                                                                                                                                 aTGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGtt                                                                                                >  2:242955/1‑251 (MQ=255)                                                                                                                                    gATTTCCGTACTCGAAATACGGGCAGACATGGCCTGCCCGGTTAGTATTATTGTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGTCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGaa                                                                                             <  2:200708/251‑1 (MQ=255)                                                                                                                                                        gggCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGtt                                                                                                <  2:82847/228‑1 (MQ=255)                                                                                                                                                        gggCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGtt                                                                                                >  1:82847/1‑228 (MQ=255)                                                                                                                                                               cATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTgg                                                                  >  1:124631/1‑251 (MQ=255)                                                                                                                                                               cATGGCCTGACCGGTTATTATTATTTTTGACACCAGACCAACTGCTACTGGTAGCGACCGGCGCTCAGCTGGAATTCCGACTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAACCCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGCTGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTAATGTTGAACTGGAAGTCGACGCGCCCCTGGTGTgg                                                                  >  2:193502/1‑251 (MQ=255)                                                                                                                                                                                     attTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTc                                            >  1:101339/1‑251 (MQ=255)                                                                                                                                                                                                                      gACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCTATAATTCAATTCTCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACg           <  1:161021/251‑1 (MQ=255)                                                                                                                                                                                                                               ccatcTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTTCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATTTTGAACTGGAAGGCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGta  <  1:182258/247‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          TTTTGCTGATATGGTTGATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTA  >  NC_000913/363053‑363525

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