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,6060CT96.0% 71.2 / ‑6.4 25*92* (TAG→TAAfrmRregulator protein that represses frmRAB operon
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (14/10);  total (15/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

CACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGC  >  NC_000913/379482‑379737
                                                                                                                            |                                                                                                                                   
cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGc                                                                                                                       >  1:721577/1‑139 (MQ=255)
      aCCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                 <  1:669362/139‑1 (MQ=255)
      aCCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                                 <  1:358569/139‑1 (MQ=255)
                              aaTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGa                                                                                         <  2:659289/139‑1 (MQ=255)
                              aaTTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGa                                                                                         >  1:633268/1‑139 (MQ=255)
                                     aCGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTg                                                                                  >  2:331283/1‑139 (MQ=255)
                                                  gTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGt                                                                     <  1:239921/139‑1 (MQ=255)
                                                     tACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTAtt                                                                                                                                 >  1:281757/1‑76 (MQ=255)
                                                     tACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTAtt                                                                                                                                 <  2:281757/76‑1 (MQ=255)
                                                                      aCGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAAc                                                 <  1:18481/139‑1 (MQ=255)
                                                                         gcagcaCGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGtt                                              >  2:39011/1‑139 (MQ=255)
                                                                                             tcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAg                          <  2:406094/139‑1 (MQ=255)
                                                                                                           tatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCt                                                                                   >  1:179494/1‑68 (MQ=255)
                                                                                                           tatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCt                                                                                   <  2:179494/68‑1 (MQ=255)
                                                                                                           tatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCa            >  1:733889/1‑139 (MQ=255)
                                                                                                             tGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGc          >  1:105464/1‑139 (MQ=255)
                                                                                                                   tAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATgg                                  <  2:705878/109‑1 (MQ=255)
                                                                                                                   tAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATgg                                  >  1:705878/1‑109 (MQ=255)
                                                                                                                   tAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTa    >  2:767322/1‑139 (MQ=255)
                                                                                                                   tAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTa    >  1:688321/1‑139 (MQ=255)
                                                                                                                   tAGATTCAGCTATTTAAGATAGGCACGAACCATTTCAATAGTGTCGTCAACGGATTGGATGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTGTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTa    >  2:46664/1‑139 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGatat                                    >  2:169362/1‑105 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGatat                                    <  1:169362/105‑1 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGc  >  2:308783/1‑139 (MQ=255)
                                                                                                                     gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGc  >  1:424451/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                   
CACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCATTAGC  >  NC_000913/379482‑379737

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: