Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,196,331 C→G intergenic (‑232/+1) thiC ← / ← rsd phosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,3310CG96.9% 236.1 / ‑0.6 65intergenic (‑232/+1)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (34/29);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTC  >  NC_000913/4196235‑4196424
                                                                                                |                                                                                             
cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCa                                                                                            <  2:617621/100‑1 (MQ=255)
cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCACCGCAAATGGGGAATTGAATGTAAATTACGAGTTAACAGCGCAGAAGTCa                                                                                            >  2:482629/1‑100 (MQ=255)
  cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAg                                                                                          >  2:798762/1‑100 (MQ=255)
   gAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGc                                                                                         >  2:1304132/1‑100 (MQ=255)
   gAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGc                                                                                         <  2:1391738/100‑1 (MQ=255)
   gAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGc                                                                                         <  1:1491841/100‑1 (MQ=255)
   gAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGc                                                                                         >  1:1733547/1‑100 (MQ=255)
    aaTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCa                                                                                        >  2:635060/1‑100 (MQ=255)
       aaaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGa                                                                                     >  2:462563/1‑100 (MQ=255)
           ggTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGtt                                                                                 <  1:806951/100‑1 (MQ=255)
            gTCTCANCCAGTTAAGGCACTCCGACAAGAAAGTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGNGTTANCAGCGCAGAAGTCAAGCAGGATGttt                                                                                <  1:1576608/100‑1 (MQ=255)
            gTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGttt                                                                                >  2:995840/1‑100 (MQ=255)
             tctcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTg                                                                               <  2:490700/100‑1 (MQ=255)
                cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAcg                                                                            <  2:1574832/100‑1 (MQ=255)
                  gCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAcgcg                                                                          >  2:1882082/1‑100 (MQ=255)
                  gCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAcgcg                                                                          >  2:225869/1‑100 (MQ=255)
                   ccAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgg                                                                         >  2:1829329/1‑100 (MQ=255)
                      gTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgggcg                                                                      <  2:1805031/100‑1 (MQ=255)
                         aaGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                   >  1:317489/1‑100 (MQ=255)
                         aaGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                   <  2:1855398/100‑1 (MQ=255)
                         aaGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                   <  1:1286127/100‑1 (MQ=255)
                         aaGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                   <  1:984138/100‑1 (MQ=255)
                            gCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                <  1:1599263/100‑1 (MQ=255)
                              aCTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAg                                                              <  1:1304138/100‑1 (MQ=255)
                              aCTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAg                                                              <  2:1296794/100‑1 (MQ=255)
                              aCTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAg                                                              <  1:995843/100‑1 (MQ=255)
                               cTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGc                                                             <  1:462562/100‑1 (MQ=255)
                                tCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAAGTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                                                            <  1:1884958/100‑1 (MQ=255)
                                 ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc                                                           >  2:1971793/1‑100 (MQ=255)
                                 ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc                                                           <  1:798764/100‑1 (MQ=255)
                                  cGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAcc                                                          >  1:250472/1‑100 (MQ=255)
                                   gACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca                                                         >  1:1986278/1‑100 (MQ=255)
                                    aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                        <  1:1238321/100‑1 (MQ=255)
                                    aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                        <  1:194296/100‑1 (MQ=255)
                                      aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                                      >  2:1049710/1‑100 (MQ=255)
                                      aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                                      <  2:98170/100‑1 (MQ=255)
                                      aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                                      >  2:1981089/1‑100 (MQ=255)
                                        gAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGa                                                    >  2:1046606/1‑100 (MQ=255)
                                                ccccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCtt                                            >  2:302414/1‑100 (MQ=255)
                                                  ccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTAt                                          <  2:1733535/100‑1 (MQ=255)
                                                          ggggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGc                                  <  2:953443/100‑1 (MQ=255)
                                                           gggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGca                                 >  2:1585765/1‑100 (MQ=255)
                                                            ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcac                                >  1:498112/1‑100 (MQ=255)
                                                            ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcac                                >  2:409805/1‑100 (MQ=255)
                                                             gCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcaca                               >  2:263281/1‑100 (MQ=255)
                                                                ttGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAg                            >  1:1230190/1‑100 (MQ=255)
                                                                  gAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                          <  2:1532116/100‑1 (MQ=255)
                                                                  gAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                          >  2:1932153/1‑100 (MQ=255)
                                                                  gAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                          >  1:1633713/1‑100 (MQ=255)
                                                                  gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                          >  1:812738/1‑100 (MQ=255)
                                                                    aTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                        >  1:1122641/1‑100 (MQ=255)
                                                                      gTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                      >  2:99944/1‑100 (MQ=255)
                                                                      gTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                      <  1:350957/100‑1 (MQ=255)
                                                                       tAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCtt                     >  1:1377070/1‑100 (MQ=255)
                                                                         aaTTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                   >  2:344585/1‑100 (MQ=255)
                                                                            tACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAg                 >  2:1274219/1‑99 (MQ=255)
                                                                               gcgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc             <  1:819308/100‑1 (MQ=255)
                                                                               gcgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc             <  1:765441/100‑1 (MQ=255)
                                                                               gcgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc             <  2:323706/100‑1 (MQ=255)
                                                                                cgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCt            >  2:727547/1‑100 (MQ=255)
                                                                                    aaCAGCGCAGAAGTCTAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGcc        <  1:1981101/100‑1 (MQ=255)
                                                                                    aaCAGCGCAGAAGTCAAGCAGGATGTTTGACGCTGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGcc        >  2:1313409/1‑100 (MQ=255)
                                                                                     aCAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCg       <  1:471237/100‑1 (MQ=255)
                                                                                      cAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGa      >  2:1487875/1‑100 (MQ=255)
                                                                                          gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc  >  1:1536512/1‑100 (MQ=255)
                                                                                                |                                                                                             
CCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTC  >  NC_000913/4196235‑4196424

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