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.4% 117.2 / ‑4.4 39*92* (TAG→TAAfrmRregulator protein that represses frmRAB operon
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/19);  total (19/20)
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.99e-01

TTTAATTAGCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCA  >  NC_000913/379473‑379727
                                                                                                                                     |                                                                                                                         
tttAATTAGCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTAttt                                                                                                                      <  1:898398/139‑1 (MQ=255)
        gCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAgg                                                                                                              >  2:866915/1‑139 (MQ=255)
         cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGa                                                                                                                  <  1:582231/134‑1 (MQ=255)
         cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGa                                                                                                                  >  2:582231/1‑134 (MQ=255)
         cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGAt                                                                                                                 >  2:306655/1‑135 (MQ=255)
         cACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGAt                                                                                                                 <  1:306655/135‑1 (MQ=255)
          aCTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCa                                                                                                            >  2:1002049/1‑139 (MQ=255)
             tCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGa                                                                                                         <  1:680493/139‑1 (MQ=255)
              cACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAATTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGaa                                                                                                        >  1:589223/1‑139 (MQ=255)
               aCCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       <  2:349861/139‑1 (MQ=255)
               aCCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       <  2:576495/139‑1 (MQ=255)
                 cTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAATTGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCa                                                                                                     <  2:589223/139‑1 (MQ=255)
                      ttCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCa                                                                                                >  1:769727/1‑139 (MQ=255)
                            ggtgCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAgtg                                                                                          <  1:1002049/139‑1 (MQ=255)
                                                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       >  1:251689/1‑85 (MQ=255)
                                                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       <  2:251689/85‑1 (MQ=255)
                                                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTg                                                            >  1:51219/1‑128 (MQ=255)
                                                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTg                                                            <  2:51219/128‑1 (MQ=255)
                                                                     aGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCAtt                                                 <  1:866915/139‑1 (MQ=255)
                                                                           tGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATGTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTc                                           <  2:769727/139‑1 (MQ=255)
                                                                               aCGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAg                                                                                                    >  1:97820/1‑78 (MQ=255)
                                                                               aCGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAg                                                                                                    <  2:97820/78‑1 (MQ=255)
                                                                                                   ctctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       <  2:985371/55‑1 (MQ=255)
                                                                                                   ctctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAAc                                                                                                       >  1:985371/1‑55 (MQ=255)
                                                                                                     ctcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTACCCGGATATGGCTTTCaa                 >  1:723217/1‑139 (MQ=255)
                                                                                                      tcGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAg                >  2:579305/1‑139 (MQ=255)
                                                                                                       cGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGata                           <  1:934786/127‑1 (MQ=255)
                                                                                                       cGCTCTTCCTCAATATGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGata                           >  2:934786/1‑127 (MQ=255)
                                                                                                                   atatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATc   <  1:351601/139‑1 (MQ=255)
                                                                                                                    tatGGTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCa  >  1:60946/1‑139 (MQ=255)
                                                                                                                        gTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCa              >  1:869056/1‑123 (MQ=255)
                                                                                                                        gTAATAGATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCa              <  2:869056/123‑1 (MQ=255)
                                                                                                                              gATTCAGTTATTTAAGATAGGCACGAACCAGTTCATTAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCaa                 <  2:659198/114‑1 (MQ=255)
                                                                                                                              gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGt                                            <  2:146654/87‑1 (MQ=255)
                                                                                                                              gATTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGt                                            >  1:146654/1‑87 (MQ=255)
                                                                                                                              caTTCAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCaa                 >  1:659198/2‑114 (MQ=255)
                                                                                                                                  cAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAAc                                       >  1:200511/1‑88 (MQ=255)
                                                                                                                                  cAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGata                           <  2:176684/100‑1 (MQ=255)
                                                                                                                                  cAGTTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGata                           >  1:176684/1‑100 (MQ=255)
                                                                                                                                     |                                                                                                                         
TTTAATTAGCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAATATGGTAATAGATTCAGCTATTTAAGATAGGCACGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCA  >  NC_000913/379473‑379727

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