Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,880,635	T→C	L138L (TTA→TTG)	casD ←	type I‑E CRISPR system Cascade subunit CasD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,880,635	0	T	C	100.0%	116.8 / NA	36	L138L (TTA→TTG)	casD	type I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/17);  total (19/17)

CCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGA  >  NC_000913/2880393‑2880877                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   cctcctTTAATCACATACCATTCGCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGc                                                                                                                                                                                                                                            <  1:221786/251‑1 (MQ=255)       ttAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGCTGCCTGACATGTCCCCAAAAAACGCGGGTGCGTTAGTGGGCAACTTCTCCGCCCCAGGCAAGGTGTAAACCGAGGCTTCAATACTGCTTTTTc                                                                                                                                                                                                                                      >  2:235110/1‑251 (MQ=255)                tACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTg                                                                                                                                                                                                                              >  1:465594/1‑250 (MQ=255)                           ggAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTATCTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTGAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCg                                                                                                                                                                                                                  <  1:281289/251‑1 (MQ=255)                                                 gagATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGAACCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGAAAGGTGTATATCGAGTCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCatgg                                                                                                                                                                                            >  2:383780/3‑248 (MQ=255)                                                 gTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGCGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCag                                                                                                                                                                                              >  2:35769/1‑249 (MQ=255)                                                 gTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCAAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGAAAGGTGAATACCGAGGCTTCAATAATGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGGGTTAACAaga                                                                                                                                                                                             >  1:462941/1‑250 (MQ=255)                                                  gtatcaacGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATGCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCAGCGTTGCATGGGGTGTTAACCagaga                                                                                                                                                                                           <  2:439195/249‑1 (MQ=255)                                                      catcGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGa                                                                                                                                                                                       <  1:133284/251‑1 (MQ=255)                                                            ttCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTACGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTa                                                                                                                                                                                  >  1:18390/1‑250 (MQ=255)                                                                gcgcttgcCGTAAATTTTAAATGCTTCCCTATAACTGGTTCCTTACTATATATATCGCCGCCAACGGTGTCATAATTTAATAGCGTCTTCTGAGGATACGATGACTGACATGTTCCCAAAAAAAGCGGTTGTTTTAGTGGGAAACGTCTCCGTCCCAGGTAAGGTGTATGCTTAGGCTTCATTACTGCTTTTTCAAGGTATGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAagga                                                                                                                                                                             <  1:303715/245‑1 (MQ=255)                                                                gcgcgcgCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAagga                                                                                                                                                                             <  1:35769/250‑1 (MQ=255)                                                                 cgcgcgCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATATCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAagga                                                                                                                                                                             <  1:383780/250‑1 (MQ=255)                                                                                    atgatgCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTcgc                                                                                                                                                         >  1:6939/1‑251 (MQ=255)                                                                                                   tGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATgg                                                                                                                                                                                                          <  2:334184/187‑1 (MQ=255)                                                                                                   tGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATgg                                                                                                                                                                                                          >  1:334184/1‑187 (MQ=255)                                                                                                                    tatCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCaa                                                                                                                          >  1:194878/1‑250 (MQ=255)                                                                                                                       cgccgcCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCa                                                                                                                      <  1:49210/251‑1 (MQ=255)                                                                                                                              aaCGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAt                                                                                                               >  1:416070/1‑251 (MQ=255)                                                                                                                                  ggCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAtcttc                                                                                                           >  1:424433/1‑251 (MQ=255)                                                                                                                                    cTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTtctc                                                                                                         >  1:22617/1‑251 (MQ=255)                                                                                                                                     tCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCa                                                                                                                                                                                                                                     <  2:98502/126‑1 (MQ=255)                                                                                                                                     tCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCa                                                                                                                                                                                                                                     >  1:98502/1‑126 (MQ=255)                                                                                                                                        tAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCgcgc                                                                                                     >  2:125095/1‑251 (MQ=255)                                                                                                                                         aaTTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCt                                                                                                    <  2:424433/251‑1 (MQ=255)                                                                                                                                         aaTTTAATAGCGCCGTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCt                                                                                                    <  1:49303/251‑1 (MQ=255)                                                                                                                                                           gAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCGTCTCGCGCTCCAAGGACTGTATGGTaa                                                                                  <  1:256813/251‑1 (MQ=255)                                                                                                                                                                      gCCTGACATGTCCCCAAAAAAAGCGGGGGTGGTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCGTTTTCAAGTTCTGAGATAACCATCGTTTCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAAcc                                                                        <  2:465594/250‑1 (MQ=255)                                                                                                                                                                                   ccAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGacaca                                                          <  2:6939/251‑1 (MQ=255)                                                                                                                                                                                        aaaaGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGCCACacgacg                                                     >  2:346449/1‑251 (MQ=255)                                                                                                                                                                                             cGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGt                                                <  2:416070/251‑1 (MQ=255)                                                                                                                                                                                                 tgtgtTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCa                                              >  1:205687/1‑249 (MQ=255)                                                                                                                                                                                                              aaCTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCg                               >  1:77239/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTACTCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCCCTGCAAATTGca            >  2:356864/1‑250 (MQ=255)                                                                                                                                                                                                                                         aCCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACActctct    >  1:144415/1‑251 (MQ=255)                                                                                                                                                                                                                                           cGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAAACAACGGTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATAATCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGGCACACGACGAACGTTAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGa  <  1:125095/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGA  >  NC_000913/2880393‑2880877

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