Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,885,998	A→G	S408P (TCA→CCA)	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,998	0	A	G	100.0%	102.8 / NA	32	S408P (TCA→CCA)	cas3	CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/21);  total (11/21)

ATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAA  >  NC_000913/2885871‑2886126                                                                                                                                |                                                                                                                                 aTTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAAcc                                                                                                                             <  1:423218/133‑1 (MQ=255)     aaaCACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTg                                                                                                                         >  1:43308/1‑133 (MQ=255)      aaCACTTCGAGCAATGCCCAAACCCCGGATAACGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGa                                                                                                                        <  2:205864/133‑1 (MQ=255)                ccAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGaa                                                                                                              <  2:308431/133‑1 (MQ=255)                   aTTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATGCCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGGTTTGTGGAAACTACTGACAACTCTGAAccc                                                                                                           >  1:218488/1‑133 (MQ=255)                     tCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCAc                                                                                                         >  1:367402/1‑133 (MQ=255)                         aaaCCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGctt                                                                                                     <  1:297741/133‑1 (MQ=255)                              aCGGATAAAGCGGTGTTTAACTGGCAATACCGCTACCAACACCTGATCAATCGCGCAAACGCCGATTTGCCCAAGAAACACTTTCCTATTGTTTTGCGGCAACCACTGACAACGCTGAACCCACGcctcttct                                                                                                >  2:128655/1‑133 (MQ=255)                                   tAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGccc                                                                                           <  2:367402/133‑1 (MQ=255)                                    aaaGCGGTGTTTAACGGGAAATCCCGATATGAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGcccc                                                                                          <  1:162411/133‑1 (MQ=255)                                           tgtTTAACTGGCAAGACCGATATCAACCCCTGATCAATCGTGCAAACGCCGCTTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAg                                                                                   >  2:337713/1‑133 (MQ=255)                                            gtTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGt                                                                                  >  1:496480/1‑133 (MQ=255)                                                   tGGCAATACCAATATTAACACCTGGTCAATCGTGCAAACGCGGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATcgcg                                                                           <  2:218488/133‑1 (MQ=255)                                                     gaaaTACCGAAATCAACACCTGAGCAATCGTGTAAACGCCGGTTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATcgcgcg                                                                         <  2:240915/131‑1 (MQ=255)                                                     gCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATcgcgcg                                                                         <  1:457734/133‑1 (MQ=255)                                                     gCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTCGTGGCAACCACTGACAACAGTGAACCCAAGCTGCATGTTGCCCCTGTTCATGAATcgcgcg                                                                         >  2:226517/1‑133 (MQ=255)                                                     gCAATACCGACATCAACACCTGACGAATTGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATcgcgcg                                                                         <  2:248823/133‑1 (MQ=255)                                                               tatCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATAGCGCGTGATTTTatt                                                               >  1:174654/1‑133 (MQ=255)                                                                atCAACACCCAATCATTCGTGCAATCGCAGATTTGCCCAAGTAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTattg                                                              <  2:232690/133‑1 (MQ=255)                                                                         tGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAaga                                                     <  1:424817/133‑1 (MQ=255)                                                                          gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAagag                                                    >  2:348850/1‑133 (MQ=255)                                                                           aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACGCTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAgg                                                   <  1:337713/133‑1 (MQ=255)                                                                           aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTTACAACACTGAACCCACGCTTCTTCTTGCCCATGTTCAGGAATCGCGCGTGATTTTATTGATTGAAAGAgg                                                   >  2:379986/1‑133 (MQ=255)                                                                                  gTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTaa                                            >  1:394694/1‑133 (MQ=255)                                                                                  gTGCAAACGCCGATTTGCCCAAGAAACACTTTATTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTaa                                            >  2:281964/1‑133 (MQ=255)                                                                                                    ccAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAg                          <  1:348850/133‑1 (MQ=255)                                                                                                    ccAAGAAACACTTTCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAg                          <  2:496480/133‑1 (MQ=255)                                                                                                     cAAGAAACTCTTTCTTATTGCTGTGTGGCATCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGc                         <  1:226517/133‑1 (MQ=255)                                                                                                                 tCTTATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGAtt             >  2:445447/1‑133 (MQ=255)                                                                                                                  attATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGAttt            <  1:281964/132‑1 (MQ=255)                                                                                                                   ttATTGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTg           <  2:385621/133‑1 (MQ=255)                                                                                                                      ttGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTgggg        <  2:70015/133‑1 (MQ=255)                                                                                                                      ttGCTTTGTGGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTgggg        <  1:322559/133‑1 (MQ=255)                                                                                                                      ttGCTTTGTGGCAAACACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTCCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTgggg        <  2:504439/133‑1 (MQ=255)                                                                                                                           ttGTGGCAACCACTGGCAACACTGCACCCACGCTTCTTCTTCCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGaa   <  2:174654/133‑1 (MQ=255)                                                                                                                            tgtgGCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa  >  1:321115/1‑133 (MQ=255)                                                                                                                                |                                                                                                                                 ATTAAAACACTTCGACCAATTCCCAAACCACGGATAAAGCGGTGTTTAACTGGCAATACCGATATCAACACCTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAA  >  NC_000913/2885871‑2886126

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑