Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,883,115	G→A	P342S (CCT→TCT)	casA ←	type I‑E CRISPR system Cascade subunit CasA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,883,115	0	G	A	93.5%	87.2 / ‑6.0	32	P342S (CCT→TCT)	casA	type I‑E CRISPR system Cascade subunit CasA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (15/14);  minor base C (1/1);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCTTAAGAAAACCGG  >  NC_000913/2882890‑2883354                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 cATGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGCACTTTGCGGCGCAAT‑ATTTCTGAAt                                                                                                                                                                                                                        >  1:344374/1‑251 (MQ=255)                         tGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGTCTGTTGTATATCCCAAACCAACAGTTACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTGCAATAAGAGATGCGGGATTATTACGAGGTCCCCCCATAATCAATTCAAGAGAACTGTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                                                                                                                                                 <  1:221965/250‑1 (MQ=255)                                  tttttttAACCCCTTCTGCAAAGGTATATAACGCCTTGCGGAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTGTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACACCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATCCACAACCGCCGCCACGCGATTTCcatt                                                                                                                                                                                        <  2:403250/248‑1 (MQ=255)                                       ttaaaacattCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCATGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATACCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGCTTTCcattttaatt                                                                                                                                                                                  >  2:24989/9‑251 (MQ=255)                                                     aaaGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTAt                                                                                                                                                                    >  2:90507/1‑251 (MQ=255)                                                               aaCGCCTTGCGGAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTTCCACTTTGCCGTATTGTTGCCACCCCTGATTAAACACCAACACACCATTACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGCACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCGTATCTACCACaa                                                                                                                                                           <  2:90134/250‑1 (MQ=255)                                                                aCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACt                                                                                                                                                         <  1:90507/251‑1 (MQ=255)                                                                           aaGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCTTTTATCACATTGCCGTATTGTTGCCACCCCTTATTAAACATCAACCCATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCTCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAttt                                                                                                                                               <  2:275343/250‑1 (MQ=255)                                                                           aaGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTtg                                                                                                                                              <  1:472988/251‑1 (MQ=255)                                                                              gCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCt                                                                                                                                                    >  2:206620/1‑242 (MQ=255)                                                                              gCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCt                                                                                                                                                    <  1:206620/242‑1 (MQ=255)                                                                               cttttttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTc                                                                                                                                                                                            <  2:319735/198‑1 (MQ=255)                                                                               cTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTc                                                                                                                                                                                            >  1:319735/1‑201 (MQ=255)                                                                               cTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCCTTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTc                                                                                                                                          >  2:452822/1‑251 (MQ=255)                                                                                   tttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATAAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGACCTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCat                                                                                                                                       <  1:401628/250‑1 (MQ=255)                                                                                                ccAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGCTGGTGCGGAggt                                                                                                                         >  2:145930/1‑251 (MQ=255)                                                                                                     aGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGgtggtg                                                                                                                     >  2:116882/1‑250 (MQ=255)                                                                                                      gtctcTATTTCGTTTATCACATGGCCGTATTGTTGCCACCCCTGATTAAACATCAACACAGCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGaa                                                                                                                   <  1:399559/247‑1 (MQ=255)                                                                                                       ccACTATTTCTTTTATCACATGGCCGTTGTGTTGCCTCCCCTGATTAAACATCACCTCAGCATGACGCCGTTCAAGAATAGATGTGTGATTATTACGATAGCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCGGAATTGATTCACAACCGCCGCCACGCGATGTCCATTTTCATTTTGAATAATCTTATCTACCTCAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGTTGaaa                                                                                                                  <  2:464587/250‑1 (MQ=255)                                                                                                                    tATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAAttttt                                                                                                     >  2:358614/1‑251 (MQ=255)                                                                                                                         caTTGCCGTATTGTTGCCCCCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGATCTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCTGAGGTGGTGAAAGCAAGAAATTTTTCCTCa                                                                                                <  1:116882/251‑1 (MQ=255)                                                                                                                                                 gataaaacatctaaccaccatGACGCCGTTCAATAAAAGAGGTTTGATGATTACTATATCCCCCCAGAATCAATGCAAGAGATCGTGGCG‑CGGAATGATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATGTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCGTTCTTGACTGTTAcc                                                                        <  2:183507/233‑1 (MQ=255)                                                                                                                                                    tAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCACTTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGCTTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTAcc                                                                        >  2:64693/1‑248 (MQ=255)                                                                                                                                                    tAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGTTGGTGCGGAAGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTAcc                                                                        >  2:145896/1‑248 (MQ=255)                                                                                                                                                     aaaCATCAACACATCATGACGCCGTTCAAGAATAGCTGCTTGATTAATACGATATACCCCCATAATCAATTCCAGAGAACTTTGCGGCGCAAT‑ATTTCTGACTTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCAATGATGGTGCGGAGGTGGTGACAGCAAGAAATTTTTCCTCAACATCCCCTTTCTTGACTGTTACCAGAc                                                                    >  2:229668/1‑251 (MQ=255)                                                                                                                                                           aaacaccccaggACGCCGCTCAAGAATAGATGCTTGATTTCTACGTTATCCCCCCCTAATCAATCCAAGAGAACCTTGCGGCGAAAT‑ATTTCGGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAGTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAgggg                                                              <  2:390180/240‑1 (MQ=255)                                                                                                                                                             acacATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGCCAAGGGGa                                                             >  1:460901/1‑250 (MQ=255)                                                                                                                                                                 atcatGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCg                                                        >  1:153304/1‑251 (MQ=255)                                                                                                                                                                          ccGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGAATGTTACCAGACAAGGGGAATGCGGATGGGGcc                                               >  2:406926/1‑251 (MQ=255)                                                                                                                                                                              tCAAGAATAGATTCTTGATTATTACGATACCCCCCCATAACCAATTCAAGAGACCTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATa                                            <  2:153304/250‑1 (MQ=255)                                                                                                                                                                                                         tatCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTCAATTTTTc                >  1:178516/1‑251 (MQ=255)                                                                                                                                                                                                                        aaTTCAAGAGGACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACAGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGCCTGATTTGTGTCCATGTTGGTGCGGTGGTGGTGAAAGCAAGAAATTTTTCTTCACCCACCCCTTTCTTGACTGTAACAAGACAAGGGGAATGCGGTTGGGGACATAGCCCATTAACTGTAACGGTAAATTTTTCCTTAAGAAAACCga  >  2:330910/1‑249 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 CATGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAAT‑ATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCTTAAGAAAACCGG  >  NC_000913/2882890‑2883354

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