Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 33.3% 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,9860GT33.3% 31.2 / 18.5 21noncoding (88/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (4/3);  total (11/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 = 1.00e+00

AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGC  >  NC_000913/4049843‑4050125
                                                                                                                                               |                                                                                                                                           
aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATg                                                                                                                                        >  4:802807/1‑149 (MQ=255)
    aaaaTACAAAAAGGTCTTTCCTAACTTAAACAAAAAGAATAAAATTAATCGGGGAAGGTTCAGGGGTTAGATGTTCTATCTTTTCAACCCTTTAATTCCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                                                                    <  3:1000669/149‑1 (MQ=255)
     aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                   <  3:802807/149‑1 (MQ=255)
                                                   gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                                                             <  1:1693765/93‑1 (MQ=255)
                                                   gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                                                             >  2:1693765/1‑93 (MQ=255)
                                                                      tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                    <  4:62790/99‑1 (MQ=255)
                                                                      tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                    >  3:62790/1‑99 (MQ=255)
                                                                             ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                           >  1:62571/1‑149 (MQ=255)
                                                                                  tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                                                                <  1:1001936/75‑1 (MQ=255)
                                                                                  tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                                                                >  2:1001936/1‑75 (MQ=255)
                                                                                  tCTTTCAGACCTTTTACTTCACGTAATCGCATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTCTTTTATTGGGCGAAAGAAAAGATACGTAATGCCTGCTGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                      >  1:1130146/1‑149 (MQ=255)
                                                                                             ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                          <  3:100965/102‑1 (MQ=255)
                                                                                             ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                          >  4:100965/1‑102 (MQ=255)
                                                                                                                aTTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACt                                                                                                                                     >  3:1213656/1‑40 (MQ=25)
                                                                                                                aTTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACt                                                                                                                                     <  4:1213656/40‑1 (MQ=21)
                                                                                                                    ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATc                    <  2:62571/149‑1 (MQ=255)
                                                                                                                      cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                          >  2:515620/1‑77 (MQ=255)
                                                                                                                      cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                          <  1:515620/77‑1 (MQ=255)
                                                                                                                                      ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAGTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGc  >  2:892013/1‑149 (MQ=255)
                                                                                                                                       cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCgg                                <  2:504234/118‑1 (MQ=255)
                                                                                                                                       cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCgg                                >  1:504234/1‑118 (MQ=255)
                                                                                                                                       cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg    >  2:50366/1‑146 (MQ=255)
                                                                                                                                       cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg    <  1:50366/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGC  >  NC_000913/4049843‑4050125

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