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,9660AG93.3% 41.9 / NA 15noncoding (15/59 nt)sokCsmall regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (8/6);  minor base T (0/1);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTAT  >  NC_000913/16736‑17186
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              
cccccGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCCCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑aaa                                                                                                                                                                                                          <  1:180545/250‑1 (MQ=255)
           gccctcatacgcgtgcttAAGCCGTGAAAACAGCAGCCCCCTGGC‑GGCCAGTCCGTATGTCATCCGTCA‑AGAGGT‑TTTTCCTCTCTACCAGCTCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGCCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCCCCTTGCCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATAGAGGTGGCCTCGGGTTGA‑TTGT‑AAAATTTCACTCgg                                                                                                                                                                                               <  2:353164/234‑1 (MQ=255)
                     ggATTCGTAAGCCGTGAAAACATCACCCTCC‑GTCTGGCCAGGTCGGATGTGAACC‑TCACAGAGTT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTtctc                                                                                                                                                                                     <  1:213644/251‑1 (MQ=255)
                                  gTGAAAACAGAAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAACCTCGCCTTATGCTGCTTCATTTCTCTCTTCTCCTTGACCTTACCGTCAGTAAGAGGCACTCTAACTGTGTTTAGCATATCGGGGGCCTCGGGTTGA‑‑GATGAAAAAATACCTCGCGGCTTTTCTCTATCTGCCGTTcc                                                                                                                                                                        >  2:115968/1‑250 (MQ=255)
                                        aCAGCAACCTCC‑GTCGGGCCAGTTCGGATGTGACCC‑TCACAGAGGT‑CTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAt                                                                                                                                                                  <  1:200015/251‑1 (MQ=255)
                                                  cccttcgggCCAGTTCGGATGGGAACC‑CCA‑AGGGGTTCTTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagac                                                                                                                                                         <  2:550978/244‑1 (MQ=255)
                                                   ccctCTGTCCAGGTCGGATGTGAACC‑CCACAGAGTT‑CTTTGCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagaca                                                                                                                                                        <  1:352979/248‑1 (MQ=255)
                                                                                                    aCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGGTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAACACCCGCCGCTATTATATCGCTCTCTTTAAcccct                                                                                                         >  1:6746/1‑249 (MQ=255)
                                                                                                                  cgcgggggtaaatcggATGACCATCAGGGCTACATTCACCTCCTTCTGTTGCTTCATTCCCCTCTTCCCCTTCATCTCCCGGTCGGAAAGGGCCCCTCTACATGTTTTCAGTAGATGGGGGGCCTCGG‑TTGAATAGA‑AAAATATCCCCCGGGTTTTTTTTCTATCTCCCGTTCAGATAATCCCTAAGACATATAGCCAAAACCTCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGAt                                                                                           <  1:522452/236‑1 (MQ=255)
                                                                                                                     ggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCt                                                                                        >  1:204489/1‑251 (MQ=255)
                                                                                                                                                 aCAATCATCGCCTTATGCTGCTTCATTGATATCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGACAACGCCTCGCAtttttttg                                                            >  2:395651/1‑248 (MQ=255)
                                                                                                                                                                        cATTGCTCTCTTCTACGTGAACTTACGGTCAGTAAGAGGCACTCTACCTGTGTTCAGCATATAGAGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTACTCCTGAAGACGCATCGCATATCTGAGGCGTAATTTTTTGATGATTTAat                                     >  2:470459/1‑251 (MQ=255)
                                                                                                                                                                               ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTa                               >  1:300693/1‑250 (MQ=255)
                                                                                                                                                                                       cTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGctctct                                                                                                                  >  1:62723/1‑159 (MQ=255)
                                                                                                                                                                                                           gagGCACTCTTCATGTGGTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCATGAGCCAGATAGCCTCAAGCACCCGCCGCTATTATCTCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGCCGCCGCATCGAATTCTTGTGGCGAATTTTTTTACTGCTTTAATTATTTAACTTTATTTTATATCtttatagtaatttt  >  2:150418/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              
TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTAT  >  NC_000913/16736‑17186

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