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	95.2%	57.1 / ‑5.8	21	L138L (TTA→TTG)	casD	type I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/9);  minor base A (0/1);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTG  >  NC_000913/2880463‑2880876                                                                                                                                                                             |                                                                                                                                                                                                                                                  ccGTAAATTTTAAATGATGCCCTGTAACTGCTTCCTTATTATATATATCGCCGCCCACGGGCTCATAATTTAATCGCGCATTCTGAGGATCCGATGCATGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATc                                                                                                                                                                      >  1:130846/1‑250 (MQ=255) ccGTAAATTTTAAATGATGCCCTGTAACTGCTTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATc                                                                                                                                                                      >  1:280545/1‑250 (MQ=255)              atgatgCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTcgc                                                                                                                                                        <  1:94716/251‑1 (MQ=255)                     ccTGTAACTGATTCCTCACTATATATATCGCCGCTAACGGGCTCATAATTTAATAGCGCCTTCTGAGAATCCGATACCTGACATGTCCCCAAAAAAAGCGGGAGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGTAGGCATCACATAAATATTCGCGCCATGt                                                                                                                                                 >  1:142705/1‑251 (MQ=255)                                      aCTATATATATCGACGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCACTACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAAc                                                                                                                                                                                                >  1:117188/1‑187 (MQ=255)                                                  gccgccAACGGGCTCATAATTTAATANCGCCTTNNGANGATCCGATGCCTGACATGTNCCNNNANNNANCGGGTNTGNNAGTGGGCAANTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTNCTGAGATAANNATCGTTGCATGGGGTGTTAACNAGANNGNNANGGTnaan                                                                                                                                                                              >  2:83001/1‑192 (MQ=255)                                                  gccgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGCTCCGATGCATGACATGTCCCCAAAAAACGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCAc                                                                                                                    >  1:238695/1‑251 (MQ=255)                                                                                                 ccTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACAGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGAc                                                                                           >  1:329601/1‑229 (MQ=255)                                                                                                           ccccAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTg                                                                                                                                                                                                              <  1:13062/104‑1 (MQ=255)                                                                                                                          gtgtgtTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGCGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGAc                                                                                                                                >  1:149073/1‑167 (MQ=255)                                                                                                                              nttAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTg                                                                                         >  1:223978/2‑202 (MQ=255)                                                                                                                                tAGTGTGCAACTGCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTAAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTGGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATGCGCGCCATGTTTGAGTCGGTTCATGCCTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGt                                                                                   <  1:288776/206‑1 (MQ=255)                                                                                                                                tAGTGGGCAATTTCGCCGCCTCAGGTAATGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGATATAACCATCGTGGCATGGGGTGTTAACCAGATAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTGTCATGACTTTGCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCGCGCAACCCCGTTACAGACACACGACGATCGTCAAGAATg                                       <  1:207410/250‑1 (MQ=255)                                                                                                                                tAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATg                                       <  1:83001/250‑1 (MQ=255)                                                                                                                                    gggCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTGTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCAGGTTTGAATCGTGTCATGACTTGTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCAcac                                                                           <  1:327601/210‑3 (MQ=255)                                                                                                                                     ggCAACTTCTCCGCCCCAGGTAAGGTGTATACCGATGCTTCAATACTGCTTGTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATAGTCGCGCCATGTGTGAATCGTTTCATGACTTTGCAAACCACGTTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACg                                                                            <  1:334782/208‑1 (MQ=255)                                                                                                                                         aCTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAg                                                              <  1:94166/218‑1 (MQ=255)                                                                                                                                            tctcCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAAGCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGGATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCan                           <  1:33634/250‑2 (MQ=255)                                                                                                                                               ccGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGAAACACGACGATCGTCAAGAATGAGTTCATCGCAgcgc                        >  1:289788/1‑250 (MQ=255)                                                                                                                                                                    ccGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGtt                                                                                                                                       >  1:13878/1‑118 (MQ=255)                                                                                                                                                                    ccGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAGCGTCGCGCGCTCCAAGGACTGTATGATAAGCACGCAGCCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTg  <  1:46536/251‑1 (MQ=255)                                                                                                                                                                             |                                                                                                                                                                                                                                                  CCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTG  >  NC_000913/2880463‑2880876

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