Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 30.7% noncoding (88/109 nt) spf → Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,9860GT30.7% 21.0 / 6.8 13noncoding (88/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (2/2);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAA  >  NC_000913/4049881‑4050131
                                                                                                         |                                                                                                                                                 
aaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                                                                           >  1:186421/1‑114 (MQ=255)
aaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                                                                           <  2:186421/114‑1 (MQ=255)
                    aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                 <  2:451784/88‑1 (MQ=255)
                    aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                 >  1:451784/1‑88 (MQ=255)
                                tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                          <  4:574273/83‑1 (MQ=255)
                                tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                          >  3:574273/1‑83 (MQ=255)
                                             cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                           >  1:624632/1‑117 (MQ=255)
                                             cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                           <  2:624632/117‑1 (MQ=255)
                                                                 cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                          >  4:1009921/1‑50 (MQ=37)
                                                                 cGTAATCGGATTTGGCTGAATATGTTAGCCGCCCCAGTCATTactgactg                                                                                                                                          <  3:1009921/50‑1 (MQ=25)
                                                                                             ccgacCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAGATGATCCGTAATGCCTGATGCGCTATGTTTCGCAAGCCAACCGTAGAATTGTAATCTATTGAATTTACGGGCCGGGTACGCCACATCCGGCACAAGCATTa            <  2:299024/148‑1 (MQ=255)
                                                                                                  ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                                                                >  1:666846/1‑91 (MQ=255)
                                                                                                  ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                                                                <  2:666846/91‑1 (MQ=255)
                                                                                                   cagtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                    <  2:118890/102‑1 (MQ=255)
                                                                                                   cagtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                    >  1:118890/1‑102 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                                                                        <  2:1052519/50‑1 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                                                                        >  1:1052519/1‑50 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                 >  1:629191/1‑121 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                 <  2:629191/121‑1 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                 <  2:577459/121‑1 (MQ=255)
                                                                                                   cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                 >  1:577459/1‑121 (MQ=255)
                                                                                                         gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaa  >  1:339710/1‑146 (MQ=255)
                                                                                                         gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaa  <  2:339710/146‑1 (MQ=255)
                                                                                                         |                                                                                                                                                 
AAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAA  >  NC_000913/4049881‑4050131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: