Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,880,559	C→T	A164T (GCA→ACA)	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,559	0	C	T	95.7%	65.4 / ‑6.0	23	A164T (GCA→ACA)	casD	type I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/11);  minor base A (0/1);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CCCTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCAC  >  NC_000913/2880323‑2880801                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   cccTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCa                                                                                                                                                                                                                                      <  1:141708/251‑1 (MQ=255)                                           atgatgACTTTACAGAGATACATCTATACCTCCTTTAATCACATCCCATTATCGGGCGGCAAATTGTCCAGGCAAGATGATCATCAGTCAGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACCGATTCCTCACTATATATATCGCCGCCAACGAGCTCATAATTTACTCGCGCCTTCTGAGGCTCCGCTGTCTGACATGTCCCAAAAAACAGAGGGCGTGTTCGTGGGCAACTTATCCGCCCCAGGt                                                                                                                                                                                            >  2:264484/1‑251 (MQ=255)                                              atgaGTTTACTGAGATACATCCATACCTCCTTTAATCACATGCCATTCTCTGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCCCGCCGTAAATTTAAAATGATGCCCAGTAACTGATTCCTCACTATATATATCTACTCCAACGGGCTCATAATTTAATAGAGCCTTCTGACGATCCGATGTCTGACATGTCCCAAAAAAAAGAGGGTGTGTTAGTGGGCAACTTATCCGCCCCAGGTaat                                                                                                                                                                                         <  2:92116/251‑2 (MQ=255)                                                       cTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATAcc                                                                                                                                                                                >  2:243401/1‑251 (MQ=255)                                                                aTCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGAAAGGTGATCATCGGTTCGTCGAGCGCCGTAAACTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCttt                                                                                                                                                                        >  1:266186/1‑250 (MQ=255)                                                                        ctcctTTAATCACCTACCATTCTCGGGCAGCAAATTTTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTCTATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGACCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGAGGGCAACTTCTCCGCCCCAGGCAAGGTGTCTACCGAGGCTCTAATACTGCt                                                                                                                                                               >  2:37939/1‑251 (MQ=255)                                                                             ttAATCACATACCATTCTCTGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTGCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCAGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCtttt                                                                                                                                                            <  1:271922/249‑1 (MQ=255)                                                                                 tCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGt                                                                                                                                                      >  1:186924/1‑251 (MQ=255)                                                                                   acaTACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCCGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCACGTTc                                                                                                                                                    >  2:150636/1‑251 (MQ=255)                                                                                       aCCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTgag                                                                                                                                                <  1:321/251‑1 (MQ=255)                                                                                                        aaTTGTCGAGGCAAGGCGATCGTCAGTTCGTCACGCGTCGGAAATTTAAAATGATGCCCTGTAACTGACTGCTCACTAGATAGATCGCCGACAAGGGGCTCATTATTTAATTGCGCGTTCGGAGGAGGCGATGACTGACCGGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCGCAGGTGACGTGTATACCGGGGCGTTAATACTGGTTTTTCAATTTCTGAGATAGCCAtgttggcatg                                                                                                                               <  2:143668/251‑10 (MQ=255)                                                                                                                                              cGTAAATGTTAAATGATGCCCTGTAACTGAGTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTGTAATACTGCTTTTTCAAGTTCTGAGAGAACCATCGTTGCATGGGGTGTTATCCAGAGAGCGACGGTAAAGGAGGCATca                                                                                          <  2:245050/250‑1 (MQ=255)                                                                                                                                                 aaaTTTTAAATGATGCCCTGTATCTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTTTGAGGATCCGATGTCTGTCATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTAGTCCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCAGCGTTGCATGGGGTGTTAACCAGAGAGCGACGTTAAAGGAGGCATCACATa                                                                                      <  2:264383/251‑1 (MQ=255)                                                                                                                                                     tttAAATGATGCCGTGTAAATGCTTCCTCACTATATATATTGCCGACAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCAGCCCAAGGCAAGGTGTATACCGAGGCTTTAATACTTCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAat                                                                                   >  2:240232/1‑250 (MQ=255)                                                                                                                                                      ttAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGTTCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTGCGCCCGAGGTAAGGTGTATACCGATGCTTTAATACTGCTTGTTCAAGTTCTGAGATAACCAGCGTGGCAGGGGGTGTTAACCATAGAGCGACGGTAAAGGAGGCATCACATAAata                                                                                  <  1:165512/250‑1 (MQ=255)                                                                                                                                                                       aCTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCCCATAAATATTCGCGCCATGTTTGAAt                                                                >  2:98580/1‑251 (MQ=255)                                                                                                                                                                        cTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATc                                                               >  1:52521/1‑251 (MQ=255)                                                                                                                                                                               cTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCAt                                                        >  1:168485/1‑251 (MQ=255)                                                                                                                                                                                                    aaCGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAt                                   >  2:240930/1‑251 (MQ=255)                                                                                                                                                                                                            cATAATTTAATAGCGCCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTtctc                             <  1:74481/249‑1 (MQ=255)                                                                                                                                                                                                                       aGCGCCTTCTGAGAATACGATGTCTGACATGTCCCCAAAAAAAGCGAGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAAGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCAGGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGAc                <  2:55161/251‑1 (MQ=255)                                                                                                                                                                                                                        gcgcCTTCTGAGGATCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTGCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTATCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACt               <  2:17219/251‑1 (MQ=255)                                                                                                                                                                                                                                     aTCCGATGTCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCAc  <  1:98580/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   CCCTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCAC  >  NC_000913/2880323‑2880801

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