Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,175 C→T intergenic (+165/‑314) sokC → / → nhaA small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,1750CT96.7% 91.1 / ‑4.0 30intergenic (+165/‑314)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/10);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

GAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCT  >  NC_000913/16965‑17419
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
ggggCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAat                                                                                                                                                                                                                  >  1:302391/3‑250 (MQ=255)
     cTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAGTTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTgtc                                                                                                                                                                                                            >  1:77699/1‑249 (MQ=255)
     cTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTgtc                                                                                                                                                                                                            >  1:427193/1‑249 (MQ=255)
                                  gCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAaat                                                                                                                                                                               >  1:151532/1‑251 (MQ=255)
                                   cTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAaatt                                                                                                                                                                              >  1:401456/1‑251 (MQ=255)
                                                    ttCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCt                                                                                                                                                              >  2:48024/1‑250 (MQ=255)
                                                                             ccccaaccccccGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCt                                                                                                                                    <  1:100885/242‑1 (MQ=255)
                                                                              cccAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCtt                                                                                                                                   <  1:50099/249‑1 (MQ=255)
                                                                                 aaGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTACTTTTTTAATGATTTCATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGCTTATTTTTGAACTATTGGCCTAACAAGCATCGCCGCCTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCtaccctta                                                                                                                                  >  2:301432/1‑242 (MQ=255)
                                                                                  aGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGc                                                                                                                               >  1:427373/1‑251 (MQ=255)
                                                                                           cgcTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTAt                                                                                                                      <  1:203224/251‑1 (MQ=255)
                                                                                                       gCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAggg                                                                                                          >  1:60793/1‑251 (MQ=255)
                                                                                                              ttAACCCATTCTGTTTTATCGATTCTAATCCTTAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATTATTTAATTATTTAACTTGAATTTATCTCATTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGAATATCTGGGCTTAGGGTCAc                                                                                                      <  2:396938/248‑1 (MQ=255)
                                                                                                                aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGacaaga                                                                                                                                                                                  >  1:389162/1‑168 (MQ=255)
                                                                                                                aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGacaaga                                                                                                                                                                                  <  2:389162/170‑3 (MQ=255)
                                                                                                                aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTg                                                                                                 >  1:417908/1‑251 (MQ=255)
                                                                                                                 aCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGa                                                                                                <  1:79782/251‑1 (MQ=255)
                                                                                                                     attctGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAgcgc                                                                                            >  1:37814/5‑251 (MQ=255)
                                                                                                                            tttATCGATTCTACTCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGc                                                                                     >  2:330748/1‑251 (MQ=255)
                                                                                                                            tttATCGATTCTAATCCTGACGCCGCCTCGCATTTTTGTGGCGTAATTTTTTACTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTATTTAATCCCTCAGAAATCCTAACCTTAAGCTATGATTATCTTGCCTTAGGGTCACTCGTGAGCGCTTAccgc                                                                                     >  2:91218/1‑248 (MQ=255)
                                                                                                                              tATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTGTCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCg                                                                                   >  2:406917/1‑251 (MQ=255)
                                                                                                                                               aaGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCaaaa‑aactcatctc                                                                  >  1:505923/1‑241 (MQ=255)
                                                                                                                                                    gCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAAGCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑AAC‑TCATCTCACCGCt                                                           <  2:417908/251‑1 (MQ=255)
                                                                                                                                                                           tttttAATTATTTAATTATTTAACTTTAATTTAACTCTTTAACGTAATTATTGACGACAATCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGATTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCTCCTTAATCTATGATTATCTGGGCTTAGGTTCACTCGTGAGCGTTTACAGCCGTCAAAA‑AAC‑TCATCCCACCGTTGATGTCGCAAATTCTTCAAtgg                                     <  2:494893/250‑3 (MQ=255)
                                                                                                                                                                                 aTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTATAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑AAC‑TCATCTCACCGATGATGGCGCACATTATTCAATGGATCGTa                               >  2:397388/1‑250 (MQ=255)
                                                                                                                                                                                            attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTTGATTATTTTTGAAATATTGTCCAAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCGTAAGCTATGATTATCGGGGCTTAGGGTCACTCGTGAGCGCTGTCAGCCGTCAAAA‑AAC‑ACATCTCTCCGCGGATGGCGGAAATTCTTCAATGGATCGTAAAAAACGagtt                    <  2:308009/250‑4 (MQ=255)
                                                                                                                                                                                            attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAACCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGagtt                    <  2:151532/250‑4 (MQ=255)
                                                                                                                                                                                                 aCTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTACTTAATCCCTCAGAAATCCTCACCTTAAGCTATAATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAACA‑C‑TAATATCACCGCTGATGGCACAAATTCTTCAATGGATCGTAAAAAACgtgtaattcct              >  2:278499/1‑241 (MQ=255)
                                                                                                                                                                                                         ttATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGTAATCTTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑AAC‑TCATCTCACCGCTGATGACGCAAATTCTTCAATGGCTCGTAAAAAACGAGTTATTCATACACtat       >  2:426814/1‑250 (MQ=255)
                                                                                                                                                                                                             ctctTTATCGCACATATTGACGACAAGCTGGATTATTTTTGAAATATTGTCATAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCCCCTTAAGCTATGATTATCTGGGCTTAGGGTCACTAGTGAGCGCTTACAGCCGTCAAAA‑A‑CCTCATATCACCGCTGATGGCGCAAATTCTTCAATGGTTCGAAAACAAAGAGTTATTCCTACACTATCATCt  >  2:378540/1‑251 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCT  >  NC_000913/16965‑17419

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