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_0009132,883,1150GA90.0% 45.9 / ‑3.9 20P342S (CCT→TCT) casAtype I‑E CRISPR system Cascade subunit CasA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (9/9);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCTTAAGAAAACCGGTATA  >  NC_000913/2882893‑2883358
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
gAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTg                                                                                                                                                                                                                          >  1:117890/1‑250 (MQ=255)
                   tttGAATTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcg                                                                                                                                                                                                       <  1:190432/250‑1 (MQ=255)
                                 ttttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGAAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCcattttc                                                                                                                                                                                         <  1:159613/250‑1 (MQ=255)
                                      aaCCCTTCTGAAAAGGTATATACCGCTTTTGGTAAGGCTTTTTTTTATCCAAAACTACCATTCACTAGTGCGTTTAACCCATTGCCGTATTGTTGCCACCCCTGATAAAACATCAACAAATCATAACACCGTTCAAGAATAGATGCTTGATTATTACTATATCCCCCCAGAATCAATTCAAGAGAACTGTGCGGCGCAATATTTTTGAATTGATTAACAACCGCCGCCACGCGATTTCCATTTTCATTTTg                                                                                                                                                                                   <  2:89471/251‑1 (MQ=255)
                                         ccTTCTGCAAAGGTATATAACGCCTTGCGTAAGTCTGTTTTATATCCCAAACCAACATTCACTATTTCGTGTATCACATTGCCGTATTGTTGCCCCCCCTGATTAAACATCAACACATCATGACGCCTTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTTAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaat                                                                                                                                                                                <  1:198985/251‑1 (MQ=255)
                                            tCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaataat                                                                                                                                                                             >  1:112667/1‑251 (MQ=255)
                                            tCTGCAAAGGGATATAACGCCTTGCGTAAGGCTGTTTTATATCCCACACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTGGCCACCCCTGATTAAACATCACCACATCCTGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATACCCCCATAATCAATTCAAGAGAACTTTGCGACACAATATTTCTGAATTGATTCACAACCCCCGCCACGCGATTTCCCTTTTCATTTTGaataat                                                                                                                                                                             >  2:203332/1‑251 (MQ=255)
                                                 aaaGGAATATAACGACTTGCGTAAGGCTGTTTTATATCCCCAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGACACACCTGATTAAACATCAACACATAATGACGCCGTTCAAGATTGGATGCTTGCTTATTACGATATCCACCCATAATCCATTCAAGAGGACTTTGCGGCGCACTATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATATTCTTAt                                                                                                                                                                        >  2:88632/1‑251 (MQ=255)
                                                           aaCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACaa                                                                                                                                                               <  1:194438/250‑1 (MQ=255)
                                                           aaCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAAc                                                                                                                                                              >  1:131494/1‑251 (MQ=255)
                                                                                 tataTCCCAACCAAACAGTCTCTATTTCGTTAATCCAATTGCCTTATGTTTGCCACCCCTGATTAAACATCAACAAATCAGGACTCCTTTAAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTTAATGGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCatga                                                                                                                                         <  2:117890/250‑1 (MQ=255)
                                                                                         aaaccaaaagaaaccaTTTCGTTTATCACATTGCAGTATTGTTGACACCCCTGATTAAACATCAACACATCATGACGCCGTTCCAGAGTAGATGACTGATTATTACGATCTCCCCCCATAATCAATTCACGAGAACTTTGCGGCGCACTAATTCTGAATTGATTCACAACCGCCGCCACGCGATTTCAATTTTCATTTTGAATAATATTATCTACCACAAATAGGCTGATATTTGTCTATGATGGTGCgg                                                                                                                                 >  2:142746/16‑250 (MQ=255)
                                                                                                                         cccGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCt                                                                                                <  1:80479/250‑1 (MQ=255)
                                                                                                                                           cTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGtt                                                                               >  1:373445/1‑250 (MQ=255)
                                                                                                                                                                         ttcaaaaaaataTGCTGGATAATTACGATCTCCCCCCAAAATCAATTCAATAGAACTTTGCGGCGCAATATTCCTGAATTGAGTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCAt                                                 <  2:373445/239‑1 (MQ=255)
                                                                                                                                                                                                       tCCCCCCATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCt                  >  1:81025/1‑251 (MQ=255)
                                                                                                                                                                                                        ccccccATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCtt                 >  1:6453/1‑251 (MQ=255)
                                                                                                                                                                                                            ccATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAAttt                                                                                                           <  1:338254/157‑1 (MQ=255)
                                                                                                                                                                                                            ccATAATCAATTCAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTAATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCAATGATGGTGCGGCGGTGGTGAAAGCAAGAAAttt                                                                                                           >  2:338254/1‑157 (MQ=255)
                                                                                                                                                                                                                        cAAGAGAACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCTTAAGAAAACCGGtata  >  1:141204/1‑250 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
GAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCGGAGGTGGTGAAAGCAAGAAATTTTTCCTCAACCTCCCCTTTCTTGACTGTTACCAGACAAGGGGAATGCGGATGGGGCCATAGCCCATTAACTGTAAAGGTAAATTTTTCCTTAAGAAAACCGGTATA  >  NC_000913/2882893‑2883358

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: