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,9660AG96.4% 76.1 / ‑5.3 28noncoding (15/59 nt)sokCsmall regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/12);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCAT  >  NC_000913/16726‑17177
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
aTGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATTACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt                                                                                                                                                                                                             >  1:179370/1‑251 (MQ=255)
    gATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATTACGATCCAGGCGACAATCAACGCCTTATGCTTCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAATAAGAGGCACTCTCCATGTGTTCAGACAATAGGGGGCCCCGGGGTGAt                                                                                                                                                                                                        >  1:99925/1‑251 (MQ=255)
        aCTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAAAGATGAAGATCAGGGGGGAAATCAACGCCTTATGCTGCTTCAATGCTTTCTTCTCCTTGACCCTACGGTCCGTAAGAGGCACTCTAAATGTGTTCAGCAAATAGGGGGCCCAGGGTTGAtagt                                                                                                                                                                                                     >  1:323244/1‑247 (MQ=255)
              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)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
ATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCAT  >  NC_000913/16726‑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: