Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,881,212 C→T A311T (GCG→ACG)  casC ← type I‑E CRISPR system Cascade subunit CasC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,881,2120CT95.2% 53.3 / ‑4.1 21A311T (GCG→ACG) casCtype I‑E CRISPR system Cascade subunit CasC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

GGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAATTTCCAGAGCCTGCTCCCTGGAGGCACCACCTAAATTTTCCTGAAGTTGAGCGAG  >  NC_000913/2880972‑2881449
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
ggTCGCGTTCCTTCAAAGGCCGGCTGCCCCAAGGCTGTCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGAc                                                                                                                                                                                                                                     <  1:52753/251‑1 (MQ=255)
             caagggTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCACTTATTACGAACCCAGGATTTTAACTGTTCTCAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCAGTTCAGACTATATCCATTGGCAACGCGATCCCACTATTGATTAAACGTCTGTATAGACGATCGCAAAAAGc                                                                                                                                                                                                                        >  2:335549/5‑251 (MQ=255)
                    cGGCTGCCCCCAGGCTTGCATTGGCCCCGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCttt                                                                                                                                                                                                                 <  1:382246/251‑1 (MQ=255)
                    cGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCttt                                                                                                                                                                                                                 <  1:303322/251‑1 (MQ=255)
                     gtcctcgcccaggatgtcATTGTCCCAGCAAGACGCAAGATCAAAATATATCTCAGGTACACGCCTCTCCATTATGACGAATCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTGGAGCAGTAATTGGTTCGACATCAGAGAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATATACGGTTGCAAAAAGCCATCTTTc                                                                                                                                                                                                                <  2:336247/234‑1 (MQ=255)
                     ggCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTc                                                                                                                                                                                                                <  1:108167/251‑1 (MQ=255)
                                              cAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCAt                                                                                                                                                                                       <  1:128757/251‑1 (MQ=255)
                                                                       tctcATTTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGTCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTGTTCAAAAGCATGTGCCATAGAAAGTTGCATATCgg                                                                                                                                                               <  2:345677/250‑1 (MQ=255)
                                                                                             tattaCGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATc                                                                                                                                         <  1:309559/250‑1 (MQ=255)
                                                                                                                                      tGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAg                                                                                                <  1:138893/250‑1 (MQ=255)
                                                                                                                                             cTTGAGCAGTAATTGGGGCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGAGGCATATCGGAGTAATTAACCATTACCATATCCGCAGGGTTCAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGCCctc                                                                                         >  2:4874/1‑250 (MQ=255)
                                                                                                                                                   cAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTAAGCTGTTTTGCTCCAGGGACCTCTGTTGcc                                                                                  >  1:90098/1‑251 (MQ=255)
                                                                                                                                                                 aCATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGaacaac                                                                    <  1:354831/251‑1 (MQ=255)
                                                                                                                                                                  cataaGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGaacaac                                                                    >  1:300231/5‑250 (MQ=255)
                                                                                                                                                                            gCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTg                                                                                                                                                                        <  1:329568/140‑1 (MQ=255)
                                                                                                                                                                            gCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGATGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCAATAGACAGTg                                                                                                                                                                        >  2:329568/1‑140 (MQ=255)
                                                                                                                                                                                      gcagcaGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAACGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAATTTCCAGAGc                                               >  1:14861/1‑251 (MQ=255)
                                                                                                                                                                                            gcTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAATTTCCAGAGCCTGCTc                                         >  1:94792/1‑251 (MQ=255)
                                                                                                                                                                                                                                    atTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAATTTCCAGAGCCTGCTCCCTGGAGGCACCACCTAAATTTTCCTGAAGTTGAGCGAg  <  1:253918/250‑1 (MQ=255)
                                                                                                                                                                                                                                      tGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAg                                                                                <  2:104684/170‑1 (MQ=255)
                                                                                                                                                                                                                                      tGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATATTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAg                                                                                >  1:104684/1‑170 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAATTTCCAGAGCCTGCTCCCTGGAGGCACCACCTAAATTTTCCTGAAGTTGAGCGAG  >  NC_000913/2880972‑2881449

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