Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,885,049	G→A	P724L (CCT→CTT)	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,885,049	0	G	A	100.0%	65.4 / NA	19	P724L (CCT→CTT)	cas3	CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/9);  total (10/9)

ATTGCCGACCCATTCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAGGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCTATAAA  >  NC_000913/2884976‑2885124                                                                          |                                                                            aTTGCCGACCCATTCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGaaaaa                                                                      >  1:523392/1‑81 (MQ=255)  ttGCCGACCCATTCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAAt                                                                     <  1:523354/81‑1 (MQ=255)           cATTCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCt                                                            <  1:523355/81‑1 (MQ=255)             ttCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTcc                                                          <  1:523356/81‑1 (MQ=255)                      aTCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGct                                                  <  1:523357/80‑1 (MQ=255)                              ccGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAatat                                         >  1:523399/1‑81 (MQ=255)                                   tcatcGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTtgc                                    >  2:523394/1‑81 (MQ=255)                                          aaaTGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTcc                             >  2:523389/1‑81 (MQ=255)                                          aaaTGCTATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTcc                             >  2:523395/1‑81 (MQ=255)                                                tATCCAGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCcaca                       >  2:523387/1‑81 (MQ=255)                                                     aGCCATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGAct                  >  2:523391/1‑81 (MQ=255)                                                       ccATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGActct                <  2:523358/81‑1 (MQ=255)                                                       ccATTGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGActc                 >  2:523406/1‑80 (MQ=255)                                                          ttGCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAAc             <  2:523359/81‑1 (MQ=255)                                                            gCCGGTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACg            <  2:523360/80‑1 (MQ=255)                                                                gTAAGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCt       <  2:523361/81‑1 (MQ=255)                                                                  aaGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCta      >  2:523398/1‑80 (MQ=255)                                                                   aGCATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCtat     >  2:523408/1‑80 (MQ=255)                                                                     cATCAAGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCTATaaa  <  2:523362/81‑1 (MQ=255)                                                                          |                                                                            ATTGCCGACCCATTCTGGCTCATCCATTTCCGCATCATCGTAAATGCTATCCAGCCATTGCCGGTAAGCATCAGGGAAAAATAAGGATGCTCCATTAAGCTCCTCAATATGTTGCTGCGTCCGCCACATGACTCTAACGTTGCTATAAA  >  NC_000913/2884976‑2885124

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