Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,348 (A)5→6 intergenic (+338/‑141) 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,3481.A92.9% 33.8 / ‑2.1 14intergenic (+338/‑141)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (10/3);  total (11/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCT  >  NC_000913/17105‑17554
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                               
cTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCACTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGTGCTTACAGCCGTCa‑aaaaactcatc                                                                                                                                                                                                           >  1:88609/1‑244 (MQ=255)
                                  ttAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATCATTGACGACAAGCTGAATTATTTTTGAAATATTGTCATTACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCTGAAATCCTCACCTTAAGCTAAGATTATCTGGGATTCGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGTTGGCGAAAATTCTTCattggatcg                                                                                                                                                                         >  1:248674/1‑243 (MQ=255)
                                                                tctcTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATc                                                                                                                                           >  1:193165/1‑251 (MQ=255)
                                                                           cAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACg                                                                                                                                <  1:182297/251‑1 (MQ=255)
                                                                                     cgacAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTcc                                                                                                                                                       >  1:15264/1‑218 (MQ=255)
                                                                                             tGGATTATTTTTGGAATATTGTCCTAACAAATAACGACGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCACAAAA‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTACTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTa                                                                                                               >  2:96204/1‑250 (MQ=255)
                                                                                                                    ctaacaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCAc                                                                                      >  2:226260/12‑251 (MQ=255)
                                                                                                                                 ccGACTGACAAGACATTAATTATTAATTTTACTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATATGGACTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGG‑TACAAATAGTACAAACGATCTATTCACCTGAAAGCGAAAt                                                                         >  1:190029/1‑251 (MQ=255)
                                                                                                                                                                         cccTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGc                                 >  1:9626/1‑251 (MQ=255)
                                                                                                                                                                         cccTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGACGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGc                                 >  2:113017/1‑251 (MQ=255)
                                                                                                                                                                             cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTATGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCACCGATTCTTTACCAGTGATGCCTc                              <  2:24472/246‑1 (MQ=255)
                                                                                                                                                                              agaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCgg                            >  2:186550/4‑251 (MQ=255)
                                                                                                                                                                                                cTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc           >  2:74961/1‑250 (MQ=255)
                                                                                                                                                                                                        aTCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCt  >  2:75436/1‑251 (MQ=255)
                                                                                                                                                                                                            gggCTTAGGGGCACTCGTTAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGTTCTATTCACCTGAAAGAGAAATAAGAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGGCTCGGGAGGCATTATTCTTATCATTGCCGCt  <  1:20650/246‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                               
CTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCT  >  NC_000913/17105‑17554

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: ‑