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

GCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGG‑TTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGAC‑AGCCTCAAGCACCCGCCGCTATTAT‑ATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATT  >  NC_000913/16720‑17147
                                                                                                                                                                                                                                                        |                                                                                                                                                                                        
gcgttttcggcctccccccccGCCGTTGCACTTCCTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGGCTGGCCAGTTCGGATGGGAACCTCACAGAGGTCTTTTCTCGTAACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGcc                                                                                                                                                                                      <  1:302872/234‑1 (MQ=255)
gcgcTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGGGGTCAAGATTGCGATCAGAGAGGCAAACACCCCCCTATGCTTCGTCATTGCTCTCTTCTCCTTTGAC‑TTCCGGTCAGTAAGAGGCACTCTACATGTGTTTAGCGTATAGGGGGcc                                                                                                                                                                                      >  2:381748/1‑251 (MQ=255)
 gtttgttggcctccccccccGCCGTGGCATTTACTCGGATTCTTAGGCCGTTAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCt                                                                                                                                                                                     <  1:462563/235‑1 (MQ=255)
   cTTATGGGATTACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGTGTTAAAGATGGAGGTCAGGGGGGCCATCATCTCCTTTTGCTGCTTCATTTCTCTCTTTTCCTT‑GACCTTTCGGTCAGTTAAAGGCACTTTACATGTGTTCCACATATAAGGGGCCTCg                                                                                                                                                                                   >  2:100885/1‑251 (MQ=255)
      ggggcacccccccccGCCGTTCCACTTACTCGGATTCGTAAGCCGGAAAAACAGCAACCCCCGTCTGGCCGGTTCGGATGTGACCCTCACAGAGGTCTTTTCCCGTTACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGG‑t                                                                                                                                                                               <  1:381748/240‑1 (MQ=255)
       gggcccccccccccGCCGTTGCACTTACTCGGTTTCGTAAGCCGTAAAAACAGCCACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGG‑tt                                                                                                                                                                              <  2:34132/241‑1 (MQ=255)
                                  cTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGCTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGG‑TTGATGGTAAAATATCACTCGGGGAtttt                                                                                                                                                   >  2:262001/1‑251 (MQ=255)
                                          cGTAAGCAGTAAAAACAGCAACCTCCGTCTCGCCATTTCGGATGTGATCCTCACGGAGGTCTTTCCCGGTTCCCA‑CCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGG‑TTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCt                                                                                                                                           <  2:105912/251‑1 (MQ=255)
                                           gTAAGCCGTGAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGACCCTCACAGAGGTCTTTTCTCGTTCCCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGG‑TTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTg                                                                                                                                          <  1:99659/251‑1 (MQ=255)
                                                                         gccatttgggaattaacccccacaGAAGTCTTTTCTCGCT‑CCACGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCAGTTCTCTCTTCTCCTT‑GTCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGTG‑TTGATAGTAAAATATCACGCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGAT‑AGcc                                                                                                           <  1:434992/231‑1 (MQ=255)
                                                                                   aTGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGCGATAAAAAACACGCTCAGGGCGGCAAACATCGCCTTATGCTGCTTTATTGCTCTCTTCTCCTT‑GAACTTACGGGCCGTAAGAGGCACTCTACATGTGTTCCGCATATAGGGGGCCTAGGG‑TTGATGGTAAAATATCACTCGTGGCTTTTTTCTTTCTGCCGTTCAGCCAATTCCTGAGACAGAC‑AGCCTCAAGCAccc                                                                                                 >  2:190323/1‑251 (MQ=255)
                                                                                                               ttcccA‑GCCCCGCCCCTGCGGCGGTAATACAGAGGCCGACCAGGGCGACAATCACCTCCTTATGCGGTTTCATGCCTCTCTTCTCCTT‑GCCTTTACGGTCGGTAAGAGGCACTCCAAATGGGTTCACAATTTAGG‑GGCCCCGGGGTTGATAAAAAAATATCCCTGGGGTTTTTTCTCAACCTCCCGTTCAGCTAATTCCGGAAACAGAT‑ACCCTCAAGCACCCGCCGCTATTAT‑ATCGCTCTCTTTAAcc                                                                     <  1:23256/247‑1 (MQ=255)
                                                                                                               ttACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGAAAATCCTCGCCTTATGCAGCTTCATTGCTCTCTTCTCCTT‑GACCTTTCGGTCCGTAAGAGGCACTCTACATGTGTTCTGCATATAGGGGGCCTCGGG‑GTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGGGCCA‑ACTAGTCTCAAGCACCCCCCCCT‑TTATTATCCCTCTCTTTTAccc                                                                    >  2:203224/1‑251 (MQ=255)
                                                                                                                          cgccACTACGGCGGTGATACAGATGACGATCAGGGGGACAATCATCGCCTTTTGCTGCTTCCTTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGG‑TTGATAGTAAAATATAACTCGGGGCTTTTCTCTTTCTGCCGTTTAGCTAATGCCTGCGACAGAT‑AGCCTCAAGCACCCGCCGCTATTAT‑ATCCCTCTCTTTAACCCAttctttttt                                                          >  2:459319/1‑244 (MQ=255)
                                                                                                                                     cggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGG‑TTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGAC‑AGCCTCAAGCACCCGCCGCTATTAT‑ATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTa                                                >  2:190485/1‑250 (MQ=255)
                                                                                                                                      ggTGATACAGAGAAGGATCGGGGCAAAAATCATCGCCTTTTGCTGCTTCATTGCTCTCTTCCCCTT‑GCCCTGAGGGTCAGAAAGGGCCCCTCTACATGTGTTCAGCATATGGGGGGCCGGGGG‑TTAATAGAAAAATATCACTCGGGTTTTTTCTTTCTCTGCCGTTCAGCAAAGGCCTGAGACAGAT‑ACCCTCAACCACCCGCCGCTATTAT‑ATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAAt                                              <  1:459319/251‑1 (MQ=255)
                                                                                                                                                                                  cctTCATGGCTCCTTTCTCCTG‑GACCTATCGGTCAGTAAGAGGCCCTCTACATGTGTTCAGCATATGGGGGGTCTCGGG‑TTGATATTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGAT‑AGCCTAAAGCATCCGCCGCTATTAT‑ATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAtt  <  2:236544/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                        
GCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCA‑GCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTT‑GACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGG‑TTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGAC‑AGCCTCAAGCACCCGCCGCTATTAT‑ATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATT  >  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: ‑