Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 16,981 G→A noncoding (30/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,9810GA91.7% 54.3 / ‑3.2 25noncoding (30/59 nt)sokCsmall regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (10/12);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCAT  >  NC_000913/16740‑17177
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                     
cgccgATGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTTCAAGTGAAGATCAGGGCGGAAATCATCGCCTTATGCCGCTTCCTTGCTCTCTTCTCCTTTACCTTTCGTGCCGTAAAAGGCACTCCACATTTTTTCCGCATTTAGGGGGCCTCGGGGTGATAATAAAata                                                                                                                                                                                              >  2:449098/1‑250 (MQ=255)
   cgTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACCAATGACGATCAAGGCGACAATCATCGCCCTATGCTTCTTCATTGCCCTCTTCTCCCTGACCTTTCCGTCAGTAAGCGGCACTCCACATGTGTTCAGCCTATAGGGGGCCTCGGGTTGATGAGTAAATATCAc                                                                                                                                                                                          >  1:459687/1‑251 (MQ=255)
    gTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGGCGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCCGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGGTGATAGTAAAATATCAcc                                                                                                                                                                                         >  2:11312/1‑250 (MQ=255)
     ttgcaattACTCGGATTCGTAAGCCGTGAAATCAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAACGCCGCCACTACGGCGGTGATGAAGATGAAGGTCAGGGCGACACTCATCGCCTCATGCTGCTTAATTGCTCTCTTCTCCTTGCCCTTACGGTCCCTAAAAGGCACTCTACCTGTGTTCGACATATAGGGGGCCTCGGGTTGTTAGTAAAAAATCACTc                                                                                                                                                                                        >  1:471609/7‑251 (MQ=255)
     ttcagctTACTTGGACTCTTAAGCCGTAAAAACAGCAACCTCCGGCGGGCCAGTTCGGAGGTGAACCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCAc                                                                                                                                                                                          <  1:466798/244‑1 (MQ=255)
      tGCTCTTATCCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGGACCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTc                                                                                                                                                                                        <  1:124858/250‑1 (MQ=255)
         ccgacgctctccagattcgCCGAGAAAACAGCAACCCCCGTCGGCCCAGTGCGGATGTGAACCCCACAGAGGTTTTTTCTCTTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCgggg                                                                                                                                                                                    <  2:39072/233‑1 (MQ=255)
               ccGGATTCGTAAGCCGTAAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTtc                                                                                                                                                                              <  1:65341/250‑1 (MQ=255)
                cGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTACAGATGACGATCAGGGCGACAATCATCTCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCCTACGGTCAGTAAAAGGCCCTCTACCTGTGTTCAGCATTTAGGGGGCCCCGGGTTGGTCGTAAAATATCACTCGGGGCtttttt                                                                                                                                                                             >  2:239557/1‑249 (MQ=255)
                cGGAGTCGTAAGCCGGGAAACCAGCACCCGCCGTCTGGCCAGTTCGGATGTGATCTTCACAGGGGTCTTGTCCCGTGACCAGCGCCGCCACTACGGCTGTGATACAGAGGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCGTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTtct                                                                                                                                                                             <  1:442177/251‑1 (MQ=255)
                          aGCCGTAAAAACAGCACCCTCGGTCTGGCCATTGCGGGTTGGAACCTCACAGAGTTCTTTCCTCGTACCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCGACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCg                                                                                                                                                                   <  1:448148/251‑1 (MQ=255)
                           gCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCCTATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCggggggatcggaagagcacacc                                                                                                                                                                  >  1:39072/1‑233 (MQ=255)
                                   accagcagCCTCCGTCGGGCCAGTTCGTATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTaa                                                                                                                                                          <  1:58093/249‑1 (MQ=255)
                                                                 gtgAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGCTCAGGGCGACAATCATCGCCTTTTGCTGCTTCCTTGCTCTCTTCTCCTTGACCTTACGGTCCGTAAGAGGCACTCTACCTGTGTTCAGCATTTAGGGGGCCTCGGGGTGGTAGTAAAATATCACTCGGGGCTTTTTTCTATCTGCCGTTCAGCTAATGCCTTCGACAGATAGACTCAAGCACccgc                                                                                                                            >  2:381066/1‑251 (MQ=255)
                                                                 gggAACCCCACAGAGGTCTTTTCTCGTTACCGTCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACccg                                                                                                                             <  1:82526/248‑1 (MQ=255)
                                                                                                            aCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTTCCGTTCAGCTAAAGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGAtt                                                                                 >  1:369855/1‑251 (MQ=255)
                                                                                                                                    ggggCGACAACCTCCGCTTTATGCTTCTTCTTTCCTCTCTTCCCCTTGCCTTTACGGTCAGAAAGAGCCCCTCTACATGTGTACAGAATAGAGGGGGCCGGGGTTTGATAAAAAAATACCCCCGGGGGTTTTTCTTTATCTCCCGTTCAGCTAATGCCTGAGACAAATACCCACAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtt                                                         <  1:381066/250‑1 (MQ=255)
                                                                                                                                      ggCGACAATCATCGCCTTATGCTGCTTCCTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtttt                                                       >  1:91166/1‑251 (MQ=255)
                                                                                                                                                    ccTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAtttt                                         >  1:411026/1‑251 (MQ=255)
                                                                                                                                                         tgctgcTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAAt                                    >  1:152647/1‑251 (MQ=255)
                                                                                                                                                                     ccccttttccccTTGACTTTACGGTAAGTAAGGGCCCTTTTACATGTTTTCAGCATATGGGGGGCCTGGGGTTGATAAAAAAATACCACTCGGGTTTTTTCTTTATCTGCCGTTCAGCTAATGCCTGAGACAGATACCCTAAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAattat                         <  1:239557/240‑1 (MQ=255)
                                                                                                                                                                               ccttgccctgccgggcagaaaGAGCCACTCTAAATGTGTTCAGCATATGGGGGGCCTCGGTTTGATAGAAAAATCTCCCTGGGGGTTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACATATCGCCACAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAAt              <  1:11312/232‑1 (MQ=255)
                                                                                                                                                                                   gtcctgacgggaaGTAAGAGGCACTCTACTTGTTTTCAGATTATGGGGGGCCTCGGTTTGTTAGTAAAATATCCCTCGGGGTTTTTCTCTATCTCCCGTTCATCTAATCCCGGAAACAGATTGCCCCAAGCTCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTAt          <  1:229718/239‑1 (MQ=255)
                                                                                                                                                                                         gtggtcataaaGAGCCACTCTACATGTTTTCAGCGTATGGGGGCCCTCGGTTGTAGAGTAAAATATCCCTGGGGTTTTTTTTCTATCTGCCTGTCAGCTAATGTCGGAGACAGATCGCCTCAAACACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCttt    <  1:109444/242‑2 (MQ=255)
                                                                                                                                                                                           ggTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCCTTCTGTTTTATCGCTTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCtttat  >  2:201977/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                     
CGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCAT  >  NC_000913/16740‑17177

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