Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 379,606 C→T *92* (TAG→TAA frmR ← regulator protein that represses frmRAB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913379,6060CT97.8% 140.2 / ‑4.1 46*92* (TAG→TAAfrmRregulator protein that represses frmRAB operon
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/21);  total (25/21)
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.20e-01

CACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCCGCGC  >  NC_000913/379482‑379742
                                                                                                                            |                                                                                                                                        
cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGc                                                                                                                            >  2:731827/1‑139 (MQ=255)
    tCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGa                                                                                                                        <  2:1070534/139‑1 (MQ=255)
         ttttttCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCGACGGCAGCACGTGATTTCATCACTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAg                                                                                                                   <  2:276437/139‑1 (MQ=255)
                            tCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACg                                                                                                >  2:280182/1‑139 (MQ=255)
                              aaTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGa                                                                                              <  1:731827/139‑1 (MQ=255)
                               aTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTAtt                                                                                                                                      >  2:875002/1‑98 (MQ=255)
                               aTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTAtt                                                                                                                                      <  1:875002/98‑1 (MQ=255)
                                      cGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGa                                                                                      >  1:834789/1‑139 (MQ=255)
                                                                 aTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGt                                                           <  2:939537/139‑1 (MQ=255)
                                                                          cagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                                                                         >  1:136815/1‑84 (MQ=255)
                                                                          cagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                                                                         <  2:136815/84‑1 (MQ=255)
                                                                           agcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTc                                                 <  1:471401/139‑1 (MQ=255)
                                                                              aCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCg                                              >  2:905355/1‑139 (MQ=255)
                                                                                        aTCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCtt                                    >  2:332676/1‑139 (MQ=255)
                                                                                          ctctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTc                                  >  1:1392576/1‑139 (MQ=255)
                                                                                            ctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                                                               >  2:811979/1‑92 (MQ=255)
                                                                                            ctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTcgcg                                                                               <  1:811979/92‑1 (MQ=255)
                                                                                              cGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGc                              >  1:1178119/1‑139 (MQ=255)
                                                                                                ctctTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                      <  1:1091400/49‑1 (MQ=255)
                                                                                                ctctTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                      >  2:494747/1‑49 (MQ=255)
                                                                                                ctctTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                      >  2:1091400/1‑49 (MQ=255)
                                                                                                ctctTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                      <  1:494747/49‑1 (MQ=255)
                                                                                                       tCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGcc                     <  2:834789/139‑1 (MQ=255)
                                                                                                         aaTATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCAt                   <  1:905355/139‑1 (MQ=255)
                                                                                                          atatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATc                  <  2:291573/139‑1 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCt                                                                            >  1:815684/1‑78 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCt                                                                            <  2:815684/78‑1 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                                                                          >  2:926204/1‑80 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                                                                          <  1:926204/80‑1 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCa                                                         <  2:177217/97‑1 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCa                                                         >  1:177217/1‑97 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATAGCTTTCAAGCACTTCTGCCATAAGc               >  1:973479/1‑139 (MQ=255)
                                                                                                               gTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                                                                    <  1:1403/52‑1 (MQ=255)
                                                                                                               gTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTgtcgtc                                                                                                    >  2:1403/1‑52 (MQ=255)
                                                                                                                tAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCgg                                                                              >  2:341490/1‑73 (MQ=255)
                                                                                                                tAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCgg                                                                              <  1:341490/73‑1 (MQ=255)
                                                                                                                tAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTc                         <  2:896471/126‑1 (MQ=255)
                                                                                                                tAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTc                         >  1:896471/1‑126 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                                                                          <  2:434720/72‑1 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                                                                          >  1:434720/1‑72 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAg                               <  2:531914/115‑1 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAg                               >  1:531914/1‑115 (MQ=255)
                                                                                                                        tCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCcgc    >  2:211629/1‑139 (MQ=255)
                                                                                                                         cAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCcgcg   >  1:511722/1‑139 (MQ=255)
                                                                                                                         cAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCcgcg   >  1:260799/1‑139 (MQ=255)
                                                                                                                          aGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCCgcgc  >  1:1166612/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                        
CACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGCCGCGC  >  NC_000913/379482‑379742

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