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%	84.5 / NA	26	S1005* (TCG→TAG)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/8);  total (18/8)

TAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGG  >  NC_000913/363078‑363523                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        tAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATcc                                                                                                                                                                                                     >  1:274984/1‑251 (MQ=255)                aaaCGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGt                                                                                                                                                                                      <  2:18526/251‑1 (MQ=255)                  aCGGGAAGTAGGCGCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACAGCCAAAAGTTTGTGTGTTTAAAATAGTACAGAATGGATGTCCGTACGCGAAATACGGGCAGACATGGCCTGCCCGTTTATTAGTATTTTTGACACCAGACCAACTGGGAAGGGTAGCGACCGGCGCCCAGCTGGAATTCCGCCTATACTGACTGGCTCCAGGAGGCTCCGCCACCAATCCCCAGATGGAAGCCGTCg                                                                                                                                                                                    <  2:188533/251‑1 (MQ=255)                                        aaaaaaGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCcttctt                                                                                                                                                              >  1:401595/1‑251 (MQ=255)                                                               aaCCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGATATTATTATTTTTGACACCCGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAgatggcgatgg                                                                                                                                       >  1:250987/1‑251 (MQ=255)                                                                               gTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGtt                                                                                                                       <  1:15898/251‑1 (MQ=255)                                                                               gTTTGTGTTTTTTAAATAGGACATAATGGATTTCCTTACGCGAAATTCGGGCAGACATGGCCTGCCCGGTGATTATTATTTGTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTTCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGTCGATGGCTGGTTTCCATCAGtt                                                                                                                       <  1:137168/251‑1 (MQ=255)                                                                                            aaaTAGTACATCATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGCCCGTCGCTCAGCTGGAACTCCGCCTATACTGACGGGCTCCAGCAGTCGTCGCCACAAATCCCCATATGGACACCGTCGATAATCAGCCATGTGCCTTCTTCCGCGTGCAGCAGCTGGCGATAGCTGGTCTCCATCAGTTGCTGTTGACTGt                                                                                                           >  2:62092/1‑250 (MQ=255)                                                                                            aaaTAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGCCACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGt                                                                                                           >  1:178251/1‑250 (MQ=255)                                                                                                    caaaatagATTTCCTGACGCGAAAAACGGGAAGACATGACCTGCCCGGTTATTATTCTTTTTGCCACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCAGGAATTCCGCCTATACTGACAGGATCCTGGGATCGCCGCCACCAATCCCCATATGTAAACCAGCGATATTCAGCCAATTGCCTTCTTCCGCGTGCAGCCAATGGAGTTGGCTGGTTTCCATCAGTGGCTGTTggtacttgcggagtg                                                                                                  >  2:347926/8‑237 (MQ=255)                                                                                                                ccTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTgg                                                                                       >  1:256725/1‑250 (MQ=255)                                                                                                                     cgcgAAATACGTGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCTTCGATATTCAGCCATGTGCCTTCTTCCGTGTGCAGCAGATGGCGATGGCTGGTTTCGATCAGTTGCTGtt                                                                                                                 <  1:306362/219‑1 (MQ=255)                                                                                                                        gAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAgatggcgatg                                                                                                                                        >  2:216525/1‑193 (MQ=255)                                                                                                                        gAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAgatggcgatg                                                                                                                                        <  1:216525/193‑1 (MQ=255)                                                                                                                                            cTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCAt                                                           >  1:44377/1‑250 (MQ=255)                                                                                                                                                                         gACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGAc                             >  1:159691/1‑251 (MQ=255)                                                                                                                                                                          aCCAACGGGTAATGGTAGCTACCGGCGTCCATCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATAAGGAAACCGTCGATATGCAGCCATGTGCCTTCTTCCGCGTGCAGCAGAGGTCGATGGCCTGTTTCCATCAGTTGCTGTGGACTGTCGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTTTGGGCCATAATTCAATTCGCGCGTCCCCCATCGCAGAcc                            <  2:335970/251‑1 (MQ=255)                                                                                                                                                                           ccAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCg                           >  1:344644/1‑251 (MQ=255)                                                                                                                                                                                  ggagtggtggCGACCGGCGCGCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGGCGTCGCCACCAATCCCCATAGGGAAACCGTCGATATTCAGCCATGTGCCTGCGTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCAGCAGGTGCTGTTGACTGTAGCGGCGGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAAGTCAATTCGCGCGTCCCGCAGCGCAGACCTTGTTCGc                    <  2:18830/247‑1 (MQ=255)                                                                                                                                                                                    gaggggatctatcagcgcTCAGATGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATCTTCAGCCATGTGCCTTATTCCGCGTGCAGCAGATGGCGATGGCTGATTTCCATCAGTTGCTGTTGACTGTCACGGCTGATGTTGCACTGGAAGCCGCCGCGCCACTGGCGTGGGCCCAAATTCAATTCGCGCGTCCCGCAGCGAAGCCCGTTTTCGCTc                  >  2:389658/15‑251 (MQ=255)                                                                                                                                                                                      tGGTAGCGACCGGCGCTCAGCTGGAATTCCGCATATACTGACGGGCTCCAGGAGTCGTCGCCACCACTCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGGGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGGCCGTTTTCGCTCgg                >  2:262228/1‑251 (MQ=255)                                                                                                                                                                                       ggTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCggg               >  1:167548/1‑251 (MQ=255)                                                                                                                                                                                         tAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGCCTGTAGCGGCTGATGTTGAACTGGACGTCGCCGCGCCACTGGTGTGAGCCGTAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCgggc              >  2:111402/1‑249 (MQ=255)                                                                                                                                                                                         tAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGa              >  1:201571/1‑250 (MQ=255)                                                                                                                                                                                               ccGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGt        >  1:107791/1‑250 (MQ=255)                                                                                                                                                                                                    gcTCAGCTCGAATTCCGCATATACTGCCGGGCTCCAGGAGTCGTCGACCCCACTCCACATATGGAAACCGTCGATATTCAGCCATGTGACTTCTTCAGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATACTTCAATTCACGCGTCCCGCAACGCAGACCGTTTTCGCTCGGGACGACGTAAgggg  >  2:262918/1‑251 (MQ=255)                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        TAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGG  >  NC_000913/363078‑363523

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