Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 333,763 A→G L224P (CTG→CCG)  yahB ← putative DNA‑bindng transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913333,7630AG96.7% 201.1 / ‑2.7 61L224P (CTG→CCG) yahBputative DNA‑bindng transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (27/32);  total (28/33)
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.84e-01

CCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACAAAGGCTAACGGGT  >  NC_000913/333629‑333895
                                                                                                                                      |                                                                                                                                     
ccTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAAt                                                                                                                                   <  1:478897/139‑1 (MQ=255)
ccTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAAt                                                                                                                                   <  1:556783/139‑1 (MQ=255)
           ttATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGcc                                                                                                                        <  1:195828/139‑1 (MQ=255)
                    tGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTgcagca                                                                                                               >  1:420978/1‑139 (MQ=255)
                    tGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTgcagca                                                                                                               >  2:916294/1‑139 (MQ=255)
                      tCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                                             >  2:1119456/1‑139 (MQ=255)
                       cGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCATCGGAATTGACTCCTGCCCGTGCAGCAAcc                                                                                                            <  2:1065026/139‑1 (MQ=255)
                          gTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGc                                                                                                         <  2:432050/139‑1 (MQ=255)
                          gTAACCAGCAGCGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGc                                                                                                         <  1:993634/139‑1 (MQ=255)
                             aCCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCa                                                                                                      >  2:622580/1‑139 (MQ=255)
                                agcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACtt                                                                                                   >  1:88440/1‑139 (MQ=255)
                                         tGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTaa                                                                                         >  2:1096358/1‑139 (MQ=255)
                                               aTCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTCCAGCAACCAGCCCACTTTTTTATTAATCGTAt                                                                                   >  2:1514355/1‑139 (MQ=255)
                                                tCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATg                                                                                  >  2:795940/1‑139 (MQ=255)
                                                    cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                              <  1:622580/139‑1 (MQ=255)
                                                    cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                              <  1:1119456/139‑1 (MQ=255)
                                                    cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                              <  1:1096358/139‑1 (MQ=255)
                                                      cACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGc                                                                                                                         <  2:1433751/95‑1 (MQ=255)
                                                      cACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGc                                                                                                                         >  1:1433751/1‑95 (MQ=255)
                                                       aCGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                                           <  2:420978/139‑1 (MQ=255)
                                                         ggACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCCGCCCACTTTTTTATTAATCGTATGCGCGGTGtc                                                                         <  1:1514355/139‑1 (MQ=255)
                                                           aCCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcct                                                                       <  1:702351/139‑1 (MQ=255)
                                                            ccATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                                      <  1:795940/139‑1 (MQ=255)
                                                            ccATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                                      <  2:88440/139‑1 (MQ=255)
                                                            ccATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                                      <  1:916294/139‑1 (MQ=255)
                                                             cATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                                >  2:751153/1‑128 (MQ=255)
                                                             cATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                                <  1:751153/128‑1 (MQ=255)
                                                                gTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                                        >  2:1470813/1‑102 (MQ=255)
                                                                gTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                                        <  1:1470813/102‑1 (MQ=255)
                                                                gTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTTAATCGTATGCGCGGTGTCCTccac                                                                   >  1:1295551/1‑139 (MQ=255)
                                                                  aaTCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTTAATCGTATGCGCGGTGTCCTccacca                                                                 <  2:1295551/139‑1 (MQ=255)
                                                                   aTCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACttttt                                                                                                <  2:962598/107‑1 (MQ=255)
                                                                   aTCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACctttt                                                                                                >  1:962598/1‑107 (MQ=255)
                                                                       ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTtta                                                                                              <  1:839345/105‑1 (MQ=255)
                                                                       ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTtta                                                                                              >  2:839345/1‑105 (MQ=255)
                                                                       ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAatat                                                           <  2:579740/139‑1 (MQ=255)
                                                                         cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTcc                                                                                                                            <  1:847972/73‑1 (MQ=255)
                                                                         cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTcc                                                                                                                            >  2:847972/1‑73 (MQ=255)
                                                                         cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                                      <  1:278053/126‑1 (MQ=255)
                                                                         cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                                      >  2:278053/1‑126 (MQ=255)
                                                                         cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCCGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTa                                                         >  1:295070/1‑139 (MQ=255)
                                                                           aaaaaCCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAgg                                                       >  1:574953/1‑139 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                                        <  1:1063385/63‑1 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                                        >  2:1063385/1‑63 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTAt                                                                                   >  2:152620/1‑83 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTAt                                                                                   <  1:152620/83‑1 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATg                           >  1:281274/1‑139 (MQ=255)
                                                                                                       gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATg                           >  1:750892/1‑139 (MQ=255)
                                                                                                             aCATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTccacc                                                                  >  2:386067/1‑94 (MQ=255)
                                                                                                             aCATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTccacc                                                                  <  1:386067/94‑1 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCa                                                                                                      >  1:1391053/1‑49 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCa                                                                                                      <  2:1391053/49‑1 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                                <  2:1439320/70‑1 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                                >  1:1439320/1‑70 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                              >  2:326349/1‑72 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                              <  1:326349/72‑1 (MQ=255)
                                                                                                                       ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACAAAgg           >  2:653850/1‑139 (MQ=255)
                                                                                                                         gAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACAAAGGCt         <  2:574953/139‑1 (MQ=255)
                                                                                                                                ggCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACAAAGGCTAACGGGt  <  2:295070/139‑1 (MQ=255)
                                                                                                                                 gCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAgaga                                                >  1:1454504/1‑92 (MQ=255)
                                                                                                                                 gCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAgaga                                                <  2:1454504/92‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACAAAGGCTAACGGGT  >  NC_000913/333629‑333895

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: