Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 8.5% noncoding (227/245 nt) csrC → small regulatory RNA CsrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,2620AC 8.5% 101.2 / 2.2 42noncoding (227/245 nt)csrCsmall regulatory RNA CsrC
Reads supporting (aligned to +/- strand):  ref base A (16/22);  new base C (2/2);  total (18/24)
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.00e-01

GTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTC  >  NC_000913/4051126‑4051398
                                                                                                                                        |                                                                                                                                        
gTCCGGAGACGTAATGAAACGGATATGGAATCTACACTGATTGTAAAAAAAAAAAAAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCCACAGATTACTCTGTc                                                                                                                                    <  2:487785/143‑1 (MQ=255)
 tCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCg                                                                                                                                   <  2:143742/143‑1 (MQ=255)
    ggAGCCGGAATTAACCGGAACGGGGTGCACCACGGCTTGGTCCGGCTAAAGGAAAAGCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                                                                <  2:366348/143‑1 (MQ=255)
    ggAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                                                                >  2:357032/1‑143 (MQ=255)
              ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                      >  1:523051/1‑143 (MQ=255)
              ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                      >  1:220179/1‑143 (MQ=255)
              ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                      >  1:411792/1‑143 (MQ=255)
                 aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGCTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCCACAGAGTAATCTGTCGCCTTGTTTCTttgctt                                                                                                                   >  1:565219/1‑143 (MQ=255)
                          ggAATCAATACGGATTGTTCCGGCTAAAGGAAAAGCAGTGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                          <  1:448699/143‑1 (MQ=255)
                           gAATCGACACGGATTGGTACGGCTAAAGGAAAGACAGGGTGGGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                         <  2:460250/143‑1 (MQ=255)
                            aaTCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGAAAGAGTAATCTGTCGCCTTTTTTCTTTGCTTCCTTTCTGttg                                                                                                        >  1:125158/1‑142 (MQ=255)
                                    cGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                                                                <  1:357032/143‑1 (MQ=255)
                                       atcgtgccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                             <  1:289544/137‑1 (MQ=255)
                                              cGGCTAAAGGAACAGCGGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                                                      <  2:138379/143‑1 (MQ=255)
                                               ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCGAATCtat                                                                                     >  1:201629/1‑143 (MQ=255)
                                               ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAACTCtat                                                                                     >  1:160373/1‑143 (MQ=255)
                                                  gaaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                                                                  <  2:214257/142‑1 (MQ=255)
                                                          aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                                                          <  1:166848/143‑1 (MQ=255)
                                                            aCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                                                        <  2:27462/143‑1 (MQ=255)
                                                              aGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAATTCTATACTGAATAAAACGGc                                                                      <  2:301973/143‑1 (MQ=255)
                                                                   gtgtTGGCGGCCTGCCACGACTGTAACACCCGTTAACGGTTATGAGTCAGCAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                                                                 <  1:114339/143‑1 (MQ=255)
                                                                   gtggtgGCGGCCTGCACGGAGTGTAAACCGCGTTGAGGATTATGTGTCCGGAATAAAGGCGACATCGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                                                                 <  2:277782/143‑1 (MQ=255)
                                                                   ctgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                                                                 <  1:355077/142‑1 (MQ=255)
                                                                       tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGGAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAGTCCGCCGCGAATCTGTCCTGAATAAACCGGATAAAAGGCg                                                             >  2:30768/1‑143 (MQ=255)
                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                            <  2:523051/143‑1 (MQ=255)
                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                            <  2:220179/143‑1 (MQ=255)
                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                            <  2:411792/143‑1 (MQ=255)
                                                                         gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGACAAAAGGCGACAGAGTCATCTGTCGCCTCTTTTCTTTGCTTGCTTTCTGTTCGATTCCGCCGCAAATCTATACTGTATAAAACGGCTAAAATACGaa                                                           >  1:96970/1‑143 (MQ=255)
                                                                              cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCGAAAAGAGGAACGAtt                                                      >  1:239329/1‑143 (MQ=255)
                                                                                           aaGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat                                         <  2:101953/143‑1 (MQ=255)
                                                                                           aaGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAGTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat                                         <  2:465569/143‑1 (MQ=255)
                                                                                           aaGACCCGGTAAGGGTTATGAGTCAGGAAGAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat                                         <  1:255162/143‑1 (MQ=255)
                                                                                                 cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCAGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCAACTTCaa                                   >  2:86712/1‑143 (MQ=255)
                                                                                                  gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaaa                                  <  1:140499/143‑1 (MQ=255)
                                                                                                          ttATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGCAACCCTCATCTTCCAACTCACgca                          >  1:551199/1‑143 (MQ=255)
                                                                                                                 tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAgg                   >  2:523371/1‑143 (MQ=255)
                                                                                                                          aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa          <  2:126654/143‑1 (MQ=255)
                                                                                                                            ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCACTAGTCATGCAAAAg        >  1:76205/1‑143 (MQ=255)
                                                                                                                             gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCCAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACCCACGCAGCCAAGGTCATGCAAACgc       >  2:124789/1‑143 (MQ=255)
                                                                                                                               gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAGCCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgcgc     >  1:436523/1‑143 (MQ=255)
                                                                                                                               gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAACCCGTCATCTTCAAACTCACGCAGCGAAGGTCATGCAAAAgcgc     >  1:138559/1‑143 (MQ=255)
                                                                                                                                  agagTAATCTGTCGCCTTTTTGCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAGAACGGCTAAAGGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc  <  2:526508/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
GTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTC  >  NC_000913/4051126‑4051398

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: