Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 367,830 G→A Y197Y (TAC→TAT)  mhpR ← mhp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913367,8300GA98.5% 215.3 / ‑4.2 65Y197Y (TAC→TAT) mhpRmhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (40/24);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCG  >  NC_000913/367733‑367926
                                                                                                 |                                                                                                
ccAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATa                                                                                                >  1:1079550/1‑100 (MQ=255)
  aGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCTATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAAc                                                                                              <  2:388301/100‑1 (MQ=255)
    aTTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAAccc                                                                                            <  1:1565399/100‑1 (MQ=255)
     ttCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGCTCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCct                                                                                           >  1:1258959/1‑100 (MQ=255)
                aaTCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTgcgcgcg                                                                                <  2:671268/100‑1 (MQ=255)
                 aTCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGAGGTAGTTCTGTCCATAACCCTCTTTgcgcgcg                                                                                >  1:482681/1‑99 (MQ=255)
                    aCCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTgagcgcgcgc                                                                             >  2:745209/1‑99 (MQ=255)
                      ccGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATACCCCTCTTTtcgcccgcgcgc                                                                           >  2:198431/1‑99 (MQ=255)
                      ccGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGcc                                                                          <  1:1172620/100‑1 (MQ=255)
                           gttCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAAt                                                                     >  1:1245762/1‑100 (MQ=255)
                             tCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATgg                                                                   >  1:1487551/1‑100 (MQ=255)
                              cACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGc                                                                  >  1:1351581/1‑100 (MQ=255)
                                cTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCtt                                                                <  1:850826/100‑1 (MQ=255)
                                cTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCtt                                                                >  2:1519968/1‑100 (MQ=255)
                                  gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTcc                                                              <  2:819589/100‑1 (MQ=255)
                                   cgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCa                                                             <  1:1190933/100‑1 (MQ=255)
                                   cgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCa                                                             >  1:1554988/1‑100 (MQ=255)
                                    gcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAg                                                            <  2:847468/100‑1 (MQ=255)
                                       gCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCtt                                                         >  1:1581205/1‑100 (MQ=255)
                                         ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                                       >  1:212777/1‑100 (MQ=255)
                                         ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                                       <  1:1519977/100‑1 (MQ=255)
                                         ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                                       >  1:1183485/1‑100 (MQ=255)
                                           tACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAAc                                                     >  1:372683/1‑100 (MQ=255)
                                           tACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAAc                                                     >  2:1460271/1‑100 (MQ=255)
                                              ggCGATAGAGGCGATCTTCTCCTCCTGATCTCAGCCGCGGTAGTGCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACgg                                                   <  2:14763/99‑1 (MQ=255)
                                               gCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGtt                                                 >  1:492805/1‑100 (MQ=255)
                                                 gATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTcc                                               <  2:218807/100‑1 (MQ=255)
                                                 gATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTc                                                >  2:183467/1‑99 (MQ=255)
                                                  aTAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTccc                                              >  2:247997/1‑100 (MQ=255)
                                                   tAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCg                                             <  2:1581195/100‑1 (MQ=255)
                                                     gagGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                                           <  1:1128275/100‑1 (MQ=255)
                                                     gagGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                                           <  2:485787/100‑1 (MQ=255)
                                                     gagGCGATCTTCTCCACCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCCCGCGCCCCACCATGGCCTCCCGCTTTAATGCTTcc                                               >  2:1331416/1‑96 (MQ=255)
                                                      agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                          >  1:403884/1‑100 (MQ=255)
                                                      agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                          >  2:505138/1‑100 (MQ=255)
                                                      agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                                          >  1:1024392/1‑100 (MQ=255)
                                                       ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGcc                                         >  1:1227419/1‑100 (MQ=255)
                                                         cGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTAAACGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAg                                       <  2:372683/100‑1 (MQ=255)
                                                             cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGa                                   >  2:1488485/1‑100 (MQ=255)
                                                              ttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGAt                                  <  1:1362690/100‑1 (MQ=255)
                                                                ctcctcCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAg                                >  2:1449660/1‑100 (MQ=255)
                                                                ctcctcCTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAg                                >  1:696504/1‑100 (MQ=255)
                                                                 tcctccTGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGt                               <  1:1482732/100‑1 (MQ=255)
                                                                   ctcctGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtca                             <  2:133111/100‑1 (MQ=255)
                                                                     cctGATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtcatc                           >  2:381356/1‑100 (MQ=255)
                                                                        gATCCCAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCAcc                        >  2:1072214/1‑100 (MQ=255)
                                                                          tCCCAGCCGCGGTAGTTCTGTCCATACCCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCgg                      <  2:1229173/100‑1 (MQ=255)
                                                                           cccAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCggg                     >  2:1253971/1‑100 (MQ=255)
                                                                             cggCCGTGGTAGTTTTTTCCATAACCCTCTTTGCGTGCGCGCGTCATAATGACTTCCAGCTTTATCGGTTCCCTTGCCACGTGATAATCCTCACTggg                     >  2:1341073/3‑98 (MQ=255)
                                                                             cAGCCGCGGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGc                   >  1:1049993/1‑100 (MQ=255)
                                                                                 cgcgGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCggg                >  2:1436157/1‑99 (MQ=255)
                                                                                 cgcgGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGAGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCgg                 >  2:622274/1‑98 (MQ=255)
                                                                                  gcgGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAg              >  2:1171063/1‑100 (MQ=255)
                                                                                   cgGTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAgg             >  1:1158827/1‑100 (MQ=255)
                                                                                    ggTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGc            >  2:1567707/1‑100 (MQ=255)
                                                                                    ggTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGc            <  2:1183480/100‑1 (MQ=255)
                                                                                    ggTAGTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGc            <  1:81851/100‑1 (MQ=255)
                                                                                        gTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAAc        >  1:1344007/1‑100 (MQ=255)
                                                                                        gTTCTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAAc        <  1:1119459/100‑1 (MQ=255)
                                                                                           cTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACAtt     >  2:584638/1‑100 (MQ=255)
                                                                                           cTGTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACAt      >  2:409245/1‑99 (MQ=255)
                                                                                             gTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTc   >  2:138702/1‑100 (MQ=255)
                                                                                             gTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTc   <  1:1460279/100‑1 (MQ=255)
                                                                                             gTCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTc   >  2:807625/1‑100 (MQ=255)
                                                                                              tCCATAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCg  <  1:233539/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CCAGATTCAGACAGCCAATCACCCGTTGTTCACTGCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCG  >  NC_000913/367733‑367926

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