Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,886,731	T→C	I163M (ATA→ATG)	cas3 ←	CRISPR‑associated endonuclease/helicase Cas3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,886,731	0	T	C	100.0%	43.9 / NA	14	I163M (ATA→ATG)	cas3	CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/8);  total (6/8)

CAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATC  >  NC_000913/2886535‑2886975                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATa                                                                                                                                                                                                 >  1:140055/1‑250 (MQ=255)                 aCCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCAt                                                                                                                                                                                <  1:44962/251‑1 (MQ=255)                                agaACAATCAGGGGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTGGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGaaaa                                                                                                                                                                 <  1:53598/250‑1 (MQ=255)                                                                aCCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATcc                                                                                                                                 >  1:346200/1‑251 (MQ=255)                                                                               aaaaTAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGCTCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATAc                                                                                                                  >  1:288205/1‑251 (MQ=255)                                                                                                tcAGATATCCGCTCCTCGCGGGCCTGTCTATCTTGCGCATCATAAGATTACAGAGAAGCGGGGATTTCCCAACGCGACTTATCTTGATTCTGGGAGTGTACCATAAATCCATGATGTCCTGTTTCGGCCTCTACCCATTTGAACCAGGTCGCATAATGATGAGCGGCGTCGGCAAAAAAGCTGAAAAAATCCCCGAGAGGTTGCTCAGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCAt                                                                                                  <  1:388382/249‑1 (MQ=255)                                                                                                aCAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATg                                                                                                 >  1:325556/1‑251 (MQ=255)                                                                                                  aGATATCCACTCCTCACGAGCCTGTTTATCCTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAt                                                                                               <  1:151899/251‑1 (MQ=255)                                                                                                          aCTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAc                                                                                                                        <  1:161554/218‑1 (MQ=255)                                                                                                          aCTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAc                                                                                                                        >  2:161554/1‑218 (MQ=255)                                                                                                                   gAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAtt                                                                              >  1:300530/1‑251 (MQ=255)                                                                                                                                                                                           gggAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTg       <  2:325556/250‑1 (MQ=255)                                                                                                                                                                                           gggAATGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCGGCTGATTTATATTg       <  2:61970/250‑1 (MQ=255)                                                                                                                                                                                               aTGTAACATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGTAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAGTTACGTCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTGAATTGCAGCCTACTTTCTGCTGATTTATATTGGAATc  <  2:302379/251‑1 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCACGAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATC  >  NC_000913/2886535‑2886975

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