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_0009132,880,6350TC96.7% 84.3 / ‑3.9 30L138L (TTA→TTGcasDtype I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/14);  minor base A (1/0);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

AATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTG  >  NC_000913/2880401‑2880876
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
aaTCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCa                                                                                                                                                                                                                                    >  2:89127/1‑250 (MQ=255)
          catgctcGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGGAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTTACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTTCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGAt                                                                                                                                                                                                                         <  2:38462/247‑1 (MQ=255)
                  ggAAGCAAATTGTCGAGGCAAGGAGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAAATGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTATGAGGATCCGATGCCTGACATGTCCCCCAAACAAGCGGGTGTGTTAGTGAGCAACTTCTCCGCCCCAGGTACGGTGTCTACCGAGGCTTAAATACTGCTTTTTCAAGTTCTGAGATAACCATCg                                                                                                                                                                                                                 >  1:129563/1‑251 (MQ=255)
                                                     cgtcgCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGATAGCGACGGTAAag                                                                                                                                                                              <  2:232508/251‑1 (MQ=255)
                                                     cgtcgCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTaaa                                                                                                                                                                               <  1:72512/250‑1 (MQ=255)
                                                             gcCGTAAATTTTAAATGATGCCCTGTAACTGATTACTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGACGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATc                                                                                                                                                                      >  1:253645/1‑251 (MQ=255)
                                                              ccGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATc                                                                                                                                                                      >  1:98618/1‑250 (MQ=255)
                                                                         aaaTGTTGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATGTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGTAACTTCTCCTCCCCAGGTAAGTTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTc                                                                                                                                                          <  1:20159/251‑1 (MQ=255)
                                                                                      tAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGCTGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCACCTTCTCCGCCCCAGGTAAGGTGAATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTg                                                                                                                                              >  1:196818/1‑250 (MQ=255)
                                                                                      tAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTAATACTTTAATAGCGCCTTCTGAGGCTCCGCTGCCTCACATGTCCCCAAAAAAAGCGGGTCTGTTAGTGGGCAACTTCTCCACCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCCAGTTCTGAGATACCCATCGTTGCATGGGGTGTTAACCAGAGAGAGACGGTAAAGGAGGCATCACATAAATATTCGAGCCATGTTTg                                                                                                                                              >  2:16332/1‑250 (MQ=255)
                                                                                          tGATTACTCACTATATATAACGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGTATCCGATGCCAGACATGTCCCCAAAAACAGCGGGTGTGTTAGTAGGCAACTTCTCCGCCCCAGGTAAGGTGGATACCGAGGCTTCAATACTCCTTTTTCAAGTTCTGAGATAACCATCGGTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCCTCACATAAATATTCGCGACATGTTTGAATCg                                                                                                                                         >  2:301618/1‑251 (MQ=255)
                                                                                           gATTCCTCACTATATATAACGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGTGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGc                                                                                                                                                                         <  1:288961/218‑1 (MQ=255)
                                                                                           gATTCCTCAATCTATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGc                                                                                                                                                                         >  2:288961/1‑218 (MQ=255)
                                                                                                      atatatatCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGTATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCTCCCCAGGTAAGGTTTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGATCGACGGTAAAGGAGGCATCACATAAATATTCGCTCCATGTTTGAATCGTTTCATGACtt                                                                                                                              <  1:50601/250‑1 (MQ=255)
                                                                                                          atatCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCa                                                                                                                          >  1:115689/1‑250 (MQ=255)
                                                                                                                                aaTTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTCAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCt                                                                                                   >  1:105331/1‑251 (MQ=255)
                                                                                                                                      aTAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAgg                                                                                             >  1:239647/1‑251 (MQ=255)
                                                                                                                                       aaGCGCCTTCTGAGGATCCGATTCCTGACACGTCCCCAAAAAAAGCGGGTGTGTTAGTTGACAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTGCAATACTGCTTTTTCAAGTTCTGAGATAACCAGCGTTGCATGGGGTGTTAACCAGAGAGAGACGGTAAAGGAGGCAGCACATAAATAGTCGCGCCATGTTTGAATCGTTTCATGACTTGTCAAAGCACGGTAATCTTCTCGCGCTCCAAGGa                                                                                            <  2:160976/250‑1 (MQ=255)
                                                                                                                                            ccTTTTGAGGGTCCGATGCCTGACAGGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGGTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACAGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTa                                                                                       <  2:173381/251‑1 (MQ=255)
                                                                                                                                                cggAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGATCCAAGGACTGTATgtt                                                                                   <  2:226181/249‑3 (MQ=255)
                                                                                                                                                  gAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTa                                                                                  >  1:230362/1‑250 (MQ=255)
                                                                                                                                                     aaTCCAATGCCTGAAATGCCCCCAAAACAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCACGCTCTGAGATAACCATCGTTGCATGGGGCGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATc                                                                               >  2:135163/2‑250 (MQ=255)
                                                                                                                                                                  aCATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATAGTTGCCTGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGtt                                                                  >  1:219666/1‑250 (MQ=255)
                                                                                                                                                                                 aaGCGGGTGTGTTAGTGGGCAACTGCTCCGCCCCTGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTTGGAGAAAACAATCGTTGCATGGGGTGTTAAACAGAGAGCGACGGTAAAGGAGGCATCACATAAATATGCGCGCCATGTTTGAATCGTTTAATGACTGTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAAGCACGCAACCCCGTTTCAGACACacgacgtt                                                  <  2:143836/251‑3 (MQ=255)
                                                                                                                                                                                                          cgtcGCCCCAGGTTAGGTTTATACCGAGGCTTCCATAATGCTTTTTCAAGTTCTGAGATAACCAGCGCTGAATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTGTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACTCAACCCCCGTACATACACACGTCGAGCGTCAAGAATGAGtgcagcgctgc                          <  2:222642/247‑11 (MQ=255)
                                                                                                                                                                                                            ccGCCCCAGGTAAGGTGTATACCGAGGCGTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTCAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAgc                          <  2:105331/248‑1 (MQ=255)
                                                                                                                                                                                                               cactagttAAGGTGTATACCGAGGAGTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTGAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCAc                      <  2:274971/242‑1 (MQ=255)
                                                                                                                                                                                                                   aGGTAAGGTGTCTACAGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGtt                                                                                                                                                >  2:42811/1‑123 (MQ=255)
                                                                                                                                                                                                                   aGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGtt                                                                                                                                                <  1:42811/123‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTg  >  1:174368/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
AATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTG  >  NC_000913/2880401‑2880876

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