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.7% 62.5 / ‑4.0 23A311T (GCG→ACG) casCtype I‑E CRISPR system Cascade subunit CasC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/9);  total (14/9)
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.20e-01

CGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAAT  >  NC_000913/2880975‑2881395
                                                                                                                                                                                                                                             |                                                                                                                                                                                       
cgcgTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGt                                                                                                                                                                            <  1:75135/251‑1 (MQ=255)
                      ccccccAGGCTTGCATTGGCCACGCAAGTCGCAATATCAACTAAGATCTCATGTTTACGCCTCTCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGGAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATTAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCtt                                                                                                                                                      <  1:389733/250‑1 (MQ=255)
                                       ggCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCaaa                                                                                                                                      <  1:44561/250‑1 (MQ=255)
                                          ccAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGTGTCTACATCAGATAAGCTGAATTGAGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGCTTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGc                                                                                                                                   >  1:467072/1‑250 (MQ=255)
                                                   cGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATa                                                                                                                         <  1:444405/251‑1 (MQ=255)
                                                                          gTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAAcc                                                                                                                                                 >  2:235547/1‑204 (MQ=255)
                                                                          gTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAAcc                                                                                                                                                 <  1:235547/204‑1 (MQ=255)
                                                                                 ccTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGCGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGTTTAAACGTCTGTATAGATGATTGAAAAAAGCCCTCTTACGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGCGAAATTAACAAtt                                                                                           >  2:146534/1‑251 (MQ=255)
                                                                                      ccATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGATGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAACAGCCATCTTTAGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGACATTAACCATTACCat                                                                                      >  2:288868/1‑251 (MQ=255)
                                                                                                       ggATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTaaa                                                                       >  1:1900/1‑249 (MQ=255)
                                                                                                       ggATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTaaa                                                                       >  1:472871/1‑249 (MQ=255)
                                                                                                         aTTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACAATATCCGCAGGGTTAAAAGCg                                                                   >  2:150209/1‑251 (MQ=255)
                                                                                                                  gTTCTAAAGTAGGCATTTGATTAACTTGAGCAGTAATTGGGTCTAAATCAGATAAGCTGAATTGCGCAGCAGCTCTGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCTTAAGt                                                           <  2:22617/250‑1 (MQ=255)
                                                                                                                     cTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGc                                                       >  1:462401/1‑251 (MQ=255)
                                                                                                                        aaGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGt                                                    >  1:182584/1‑251 (MQ=255)
                                                                                                                         aGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCAATAGATAGTGGCATATCGGAGAAATTAACCATTCCCATATACGCAGGGTTAAAAGCGGCATAAGAACTCTGtt                                                   >  2:359271/1‑251 (MQ=255)
                                                                                                                             ggCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGc                                               >  1:130487/1‑251 (MQ=255)
                                                                                                                                   tGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAg                                          <  1:365905/250‑1 (MQ=255)
                                                                                                                                          cTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACctct                                  >  1:201365/1‑251 (MQ=255)
                                                                                                                                                                      tAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCg                                                                   <  1:208821/190‑1 (MQ=255)
                                                                                                                                                                      tAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGAATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTaccagcg                                                                   >  2:208821/1‑184 (MQ=255)
                                                                                                                                                                       aaGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAATTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGc     >  1:282290/1‑251 (MQ=255)
                                                                                                                                                                            gaagtGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGTCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCAGTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTGTTGCTCCAGGGTCCTCTGTTGCCAGCATATGAACAACATGGGTTGCAAt  <  2:247200/245‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                       
CGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCCCAATATTGATTAAACGCCTGTATAGACGGTTGCAAAAAGCCATCTTTCGCTTTAACCGCTTTTTCAAAAGCATTTGCCATAGAAAGTGGCATATCGGAGAAATTAACCATTACCATATCCGCAGGGTTAAAAGCGGCATAAGTACGCTGTTTTGCTCCAGGGACCTCTGTTGCCAGCATATGAACAACATGGGTTGCAAT  >  NC_000913/2880975‑2881395

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