| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 4,049,989 | A→C | 30.5% | noncoding (91/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_000913 | 4,049,989 | 0 | A | C | 30.5% | 26.5 / 21.7 | 23 | noncoding (91/109 nt) | spf | Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent | 
| Reads supporting (aligned to +/- strand): ref base A (8/8); new base C (4/3); total (12/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01 | |||||||||||
AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATT  >  NC_000913/4049843‑4050137                                                                                                                                                  |                                                                                                                                                    aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATg                                                                                                                                                    >  4:802807/1‑149 (MQ=255)    aaaaTACAAAAAGGTCTTTCCTAACTTAAACAAAAAGAATAAAATTAATCGGGGAAGGTTCAGGGGTTAGATGTTCTATCTTTTCAACCCTTTAATTCCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                                                                                <  3:1000669/149‑1 (MQ=255)     aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                               <  3:802807/149‑1 (MQ=255)                                                                      tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                                <  4:62790/99‑1 (MQ=255)                                                                      tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                                >  3:62790/1‑99 (MQ=255)                                                                             ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                       >  1:62571/1‑149 (MQ=255)                                                                                  tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                                                                            >  2:1001936/1‑75 (MQ=255)                                                                                  tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                                                                            <  1:1001936/75‑1 (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                <  1:50366/146‑1 (MQ=255)                                                                                                                                       cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg                >  2:50366/1‑146 (MQ=255)                                                                                                                                                tAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt                                                            >  3:308857/1‑93 (MQ=255)                                                                                                                                                tAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt                                                            <  4:308857/93‑1 (MQ=255)                                                                                                                                                  aTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAattatt  >  3:762337/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATT  >  NC_000913/4049843‑4050137 | 
| 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: ‑ |