Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 16,966 A→G noncoding (15/59 nt) sokC → small regulatory RNA antitoxin SokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091316,9660AG87.5% 29.8 / ‑3.2 17noncoding (15/59 nt)sokCsmall regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (5/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATT  >  NC_000913/16720‑17147
                                                                                                                                                                                                                                                        |                                                                                                                                                                                        
gcgttttcggcctccccccccGCCGTTGCACTTCCTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGGCTGGCCAGTTCGGATGGGAACCTCACAGAGGTCTTTTCTCGTAACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGcc                                                                                                                                                                                      <  1:302872/234‑1 (MQ=255)
gcgcTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGGTCAAGATTGCGATCAGAGAGGCAAACACCCCCCTATGCTTCGTCATTGCTCTCTTCTCCTTTGACTTCCGGTCAGTAAGAGGCACTCTACATGTGTTTAGCGTATAGGGGGcc                                                                                                                                                                                      >  2:381748/1‑251 (MQ=255)
 gtttgttggcctccccccccGCCGTGGCATTTACTCGGATTCTTAGGCCGTTAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCt                                                                                                                                                                                     <  1:462563/235‑1 (MQ=255)
   cTTATGGGATTACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGTTAAAGATGGAGGTCAGGGGGGCCATCATCTCCTTTTGCTGCTTCATTTCTCTCTTTTCCTTGACCTTTCGGTCAGTTAAAGGCACTTTACATGTGTTCCACATATAAGGGGCCTCg                                                                                                                                                                                   >  2:100885/1‑251 (MQ=255)
      ggggcacccccccccGCCGTTCCACTTACTCGGATTCGTAAGCCGGAAAAACAGCAACCCCCGTCTGGCCGGTTCGGATGTGACCCTCACAGAGGTCTTTTCCCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt                                                                                                                                                                               <  1:381748/240‑1 (MQ=255)
       gggcccccccccccGCCGTTGCACTTACTCGGTTTCGTAAGCCGTAAAAACAGCCACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGtt                                                                                                                                                                              <  2:34132/241‑1 (MQ=255)
                                  cTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGCTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGAtttt                                                                                                                                                   >  2:262001/1‑251 (MQ=255)
                                          cGTAAGCAGTAAAAACAGCAACCTCCGTCTCGCCATTTCGGATGTGATCCTCACGGAGGTCTTTCCCGGTTCCCACCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCt                                                                                                                                           <  2:105912/251‑1 (MQ=255)
                                           gTAAGCCGTGAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGACCCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTg                                                                                                                                          <  1:99659/251‑1 (MQ=255)
                                                                         gccatttgggaattaacccccacaGAAGTCTTTTCTCGCTCCACGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCAGTTCTCTCTTCTCCTTGTCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGTGTTGATAGTAAAATATCACGCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGcc                                                                                                           <  1:434992/231‑1 (MQ=255)
                                                                                   aTGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGCGATAAAAAACACGCTCAGGGCGGCAAACATCGCCTTATGCTGCTTTATTGCTCTCTTCTCCTTGAACTTACGGGCCGTAAGAGGCACTCTACATGTGTTCCGCATATAGGGGGCCTAGGGTTGATGGTAAAATATCACTCGTGGCTTTTTTCTTTCTGCCGTTCAGCCAATTCCTGAGACAGACAGCCTCAAGCAccc                                                                                                 >  2:190323/1‑251 (MQ=255)
                                                                                                               ttcccAGCCCCGCCCCTGCGGCGGTAATACAGAGGCCGACCAGGGCGACAATCACCTCCTTATGCGGTTTCATGCCTCTCTTCTCCTTGCCTTTACGGTCGGTAAGAGGCACTCCAAATGGGTTCACAATTTAGGGGCCCCGGGGTTGATAAAAAAATATCCCTGGGGTTTTTTCTCAACCTCCCGTTCAGCTAATTCCGGAAACAGATACCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAAcc                                                                     <  1:23256/247‑1 (MQ=255)
                                                                                                               ttACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGAAAATCCTCGCCTTATGCAGCTTCATTGCTCTCTTCTCCTTGACCTTTCGGTCCGTAAGAGGCACTCTACATGTGTTCTGCATATAGGGGGCCTCGGGGTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGGGCCAACTAGTCTCAAGCACCCCCCCCTTTATTATCCCTCTCTTTTAccc                                                                    >  2:203224/1‑251 (MQ=255)
                                                                                                                          cgccACTACGGCGGTGATACAGATGACGATCAGGGGGACAATCATCGCCTTTTGCTGCTTCCTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATAACTCGGGGCTTTTCTCTTTCTGCCGTTTAGCTAATGCCTGCGACAGATAGCCTCAAGCACCCGCCGCTATTATATCCCTCTCTTTAACCCAttctttttt                                                          >  2:459319/1‑244 (MQ=255)
                                                                                                                                     cggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTa                                                >  2:190485/1‑250 (MQ=255)
                                                                                                                                      ggTGATACAGAGAAGGATCGGGGCAAAAATCATCGCCTTTTGCTGCTTCATTGCTCTCTTCCCCTTGCCCTGAGGGTCAGAAAGGGCCCCTCTACATGTGTTCAGCATATGGGGGGCCGGGGGTTAATAGAAAAATATCACTCGGGTTTTTTCTTTCTCTGCCGTTCAGCAAAGGCCTGAGACAGATACCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAAt                                              <  1:459319/251‑1 (MQ=255)
                                                                                                                                                                                  cctTCATGGCTCCTTTCTCCTGGACCTATCGGTCAGTAAGAGGCCCTCTACATGTGTTCAGCATATGGGGGGTCTCGGGTTGATATTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTAAAGCATCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAtt  <  2:236544/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                        
GCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATT  >  NC_000913/16720‑17147

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