Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,886,389	C→T	K277K (AAG→AAA)	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,389	0	C	T	100.0%	116.7 / NA	32	K277K (AAG→AAA)	cas3	CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/20);  total (12/20)

GAGCATCAACTAACACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGCTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATACGTTCTCAG  >  NC_000913/2886313‑2886463                                                                             |                                                                           gagCATCAACTAACACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTAt                                                                        >  1:520986/1‑81 (MQ=255)    cATCAACTAACACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTg                                                                     <  1:520940/81‑1 (MQ=255)          cTAACACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACaa                                                               <  1:520941/81‑1 (MQ=255)            aaCACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGt                                                             >  2:520988/1‑81 (MQ=255)              cacCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTcc                                                           <  1:520942/81‑1 (MQ=255)                  tGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACt                                                        >  2:707198/1‑80 (MQ=255)                      aTTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAAc                                                    <  2:520943/80‑1 (MQ=255)                       ttGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTc                                                  >  1:520999/1‑81 (MQ=255)                       ttGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTc                                                  <  2:520944/81‑1 (MQ=255)                         gTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTcc                                                 <  2:520945/80‑1 (MQ=255)                           cTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAAt                                              <  2:520946/81‑1 (MQ=255)                             gggTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATAc                                            >  2:520992/1‑81 (MQ=255)                                ttGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCg                                         <  2:520947/81‑1 (MQ=255)                                    tAGCCATTGTCCAGTAGTGCGTGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTc                                     <  2:520949/81‑1 (MQ=255)                                    tAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTc                                     <  1:520948/81‑1 (MQ=255)                                         aTTGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATc                                >  1:520991/1‑81 (MQ=255)                                          ttGTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATcc                               <  1:520950/81‑1 (MQ=255)                                            gTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCctg                             >  2:520993/1‑81 (MQ=255)                                            gTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCct                              <  2:520951/80‑1 (MQ=255)                                            gTCCAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCct                              <  2:520952/80‑1 (MQ=255)                                               cAGTAGTGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCctgct                           <  1:520953/80‑1 (MQ=255)                                                 gtagtGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTc                        >  1:520997/1‑81 (MQ=255)                                                  tagtGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCg                       <  2:520954/81‑1 (MQ=255)                                                    gtGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGt                     <  1:520956/81‑1 (MQ=255)                                                    gtGCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGt                     <  2:520955/81‑1 (MQ=255)                                                      gCATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTcc                   >  1:520990/1‑81 (MQ=255)                                                       cATGAACACCTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCTAATACCCGGCTCGCATCCTGCTGTCGGTCCt                  >  2:520987/1‑81 (MQ=255)                                                               ccTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATAc          <  2:520957/81‑1 (MQ=255)                                                               ccTTCATAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATAc          >  1:520984/1‑81 (MQ=255)                                                                    aTAACATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATACGTtct     >  2:521001/1‑81 (MQ=255)                                                                       aCATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATACGTTCTCag  <  2:520958/81‑1 (MQ=255)                                                                       aCATCGTTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATACGTTCTCa   <  2:520959/80‑1 (MQ=255)                                                                             |                                                                           GAGCATCAACTAACACCTGTAATTGTCTGGGTTGATAGCCATTGTCCAGTAGTGCATGAACACCTTCATAACATCGCTTATTTGATACAAGTCCACTCAACTCCAATACCCGGCTCGCATCCTGCTGTCGGTCCTGGAAATACGTTCTCAG  >  NC_000913/2886313‑2886463

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