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

ATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACG  >  NC_000913/363070‑363520                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      aTGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGCATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGcc                                                                                                                                                                                                           >  1:344144/1‑251 (MQ=255)       tGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACAAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAAt                                                                                                                                                                                                     >  1:1806/1‑251 (MQ=255)                     gAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAAc                                                                                                                                                                                       >  1:58413/1‑251 (MQ=255)                         aaCGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAAATGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTc                                                                                                                                                                                   >  2:427603/1‑251 (MQ=255)                                        cccATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCAt                                                                                                                                                                      <  1:414430/249‑1 (MQ=255)                                                aaaaaaGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCcttctt                                                                                                                                                            >  1:46285/1‑251 (MQ=255)                                                            gaaaaagaaTAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTgcagc                                                                                                                                                 >  1:58105/1‑250 (MQ=255)                                                                   aaTAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAgatggcg                                                                                                                                         >  1:123656/1‑251 (MQ=255)                                                                      aaaCCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTGCCGCGTGCAGCAgatggcgat                                                                                                                                       <  1:76243/250‑1 (MQ=255)                                                                       aaCCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAgatggcgatgg                                                                                                                                     <  1:427603/251‑1 (MQ=255)                                                                            aaCATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGAc                                                                                                                                                                                                                               <  2:423804/156‑1 (MQ=255)                                                                            aaCATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGAc                                                                                                                                                                                                                               >  1:423804/1‑156 (MQ=255)                                                                                                 tttAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCACGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCCTATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTgcct                                                                                                           >  2:147018/1‑248 (MQ=255)                                                                                                               aaTGGATTTCCTTACTCGAAATTCGGGCAGACAGGGCCTGCCCTGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGGCCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTGCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGGTGGTTTCCATCAGTTTCTGTTG‑ACTGCAGCggaggatgt                                                                                            <  2:442596/251‑8 (MQ=255)                                                                                                                                           aGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGgtg                                                                >  1:121394/1‑251 (MQ=255)                                                                                                                                                 ggCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGc                                                           <  2:123656/250‑1 (MQ=255)                                                                                                                                                          gtttattattaTTTTTGACACCAGACCAACTGGTAAGGATACCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTGTAGGATTCGTCGCCACCAAGCTCCATATGGAAACCGTCGATATGCAGCCATTTGCCTTCGTCCGCGTTCAGCATATGGCGATGGCTGGGTTCCATCAGTGGCTTTTG‑ACGGTATCTTTTGATGTTTAGCTGTAAGTCGCGGCGACACTTGTGTTGGCCATatttcaa                                                 <  2:180677/249‑1 (MQ=255)                                                                                                                                                          ggTTATGATTATTGTTGACACCAGACCAACGGTTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGACTATACTGACGGGCGCCAGGAGTCGTCGCCACCTATCCCCATATGGAAACCGGCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGGGATGGCGGGTTTCCATCAATTG‑TGGTGTACGGTAGGGGCAGATGTTGAACTAGAATTCGGCTCGCCACTTGTGTGGGCTATatttcaa                                                 <  2:422739/251‑1 (MQ=255)                                                                                                                                                           gTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCGCCAATCCCCATATGGAAACCGGCGATATTCAGCCATGTGCCTTCTTCCGCGTGCATCAGATGGCGATTGCTGGTGTCCATCAGTTGCTGTTG‑ACTGTAGCGTCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATaattcaat                                                <  1:185673/251‑1 (MQ=255)                                                                                                                                                            ttattattaTTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATaattcaatt                                               >  1:331931/1‑251 (MQ=255)                                                                                                                                                                   tattTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATaattca                                                  >  1:81657/1‑241 (MQ=255)                                                                                                                                                                   tattTTTGACACCAGACCAACGGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATaattca                                                  <  2:81657/241‑1 (MQ=255)                                                                                                                                                                        ttGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCcgc                                   >  1:310139/1‑251 (MQ=255)                                                                                                                                                                                 gACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGAc                          >  1:339567/1‑251 (MQ=255)                                                                                                                                                                                    cAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGTCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCGGATGGTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAGTTCGCGCGACCCGCAGCGCAGACCGt                       <  2:60206/251‑1 (MQ=255)                                                                                                                                                                                     aaCTGGAAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGtt                      <  2:310139/251‑1 (MQ=255)                                                                                                                                                                                           tAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCt                >  1:96236/1‑251 (MQ=255)                                                                                                                                                                                                        cGGCGCTCAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTAc   >  2:105196/1‑251 (MQ=255)                                                                                                                                                                                                         ggCGCACAGCTGGAATTCCGCCTATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGCTATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCCGTTGCTGTTG‑AATGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGAGCGTCCCGCAGCGCAGAACGTTTTAGCTCGGGAAGACGAACg  >  2:442832/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      ATGTCATGTAGCCAAATCGGGAAAAACGGGAAGTAGGCTCCCATGATAAAAAAGTAAAAGAAAAAGAATAAACCGAACATCCAAAAGTTTGTGTTTTTTAAATAGTACATAATGGATTTCCTTACGCGAAATACGGGCAGACATGGCCTGCCCGGTTATTATTATTTTTGACACCAGACCAACTGGTAATGGTAGCGACCGGCGCTCAGCTGGAATTCCGCCGATACTGACGGGCTCCAGGAGTCGTCGCCACCAATCCCCATATGGAAACCGTCGATATTCAGCCATGTGCCTTCTTCCGCGTGCAGCAGATGGCGATGGCTGGTTTCCATCAGTTGCTGTTG‑ACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACG  >  NC_000913/363070‑363520

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