Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,561,049 C→T 36.3% I352I (ATC→ATT evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,561,0490CT36.3% 35.1 / 29.9 55I352I (ATC→ATTevgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base T (6/14);  total (19/36)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGC  >  minE/1560987‑1561099
                                                              |                                                  
cTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGc                                               >  1:1139391/1‑68 (MQ=255)
cTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGc                                               >  1:1550895/1‑68 (MQ=255)
cTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGc                                               >  1:2090237/1‑68 (MQ=255)
cTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGc                                               >  1:1757019/1‑68 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:1079813/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:413381/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:850665/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:2992425/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:2116561/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:1376444/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:1716603/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:1964776/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCa                                                  <  1:1890835/58‑1 (MQ=255)
       ttattaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTCCACACAATATCCa                                                  <  1:591888/58‑1 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATc                                                    >  1:1835321/1‑54 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATc                                                    >  1:717032/1‑54 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATc                                                    >  1:3335182/1‑54 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATc                                                    >  1:3101336/1‑54 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAg                                     >  1:2826853/1‑69 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAg                                     >  1:1975683/1‑69 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAg                                     >  1:1707230/1‑69 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAg                                     >  1:930588/1‑69 (MQ=255)
         attaCCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAg                                     >  1:608100/1‑69 (MQ=255)
                     aCAGGTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:2366251/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:2702613/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:3070791/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:2223725/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:1883466/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:791409/69‑1 (MQ=255)
                      cAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCggag                        <  1:133349/69‑1 (MQ=255)
                          tttAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGAt                     >  1:3374128/1‑68 (MQ=255)
                          tttAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGAt                     >  1:254356/1‑68 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:565050/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:681734/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:698674/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:739864/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:32501/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:3107849/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:2393304/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:2940034/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:1825543/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:2914537/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:1915221/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:1229613/42‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATgg                   <  1:1334282/69‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATgg                   <  1:2225227/69‑1 (MQ=255)
                           ttAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATgg                   <  1:2908700/69‑1 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGc                                               >  1:200649/1‑40 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:2968925/41‑1 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCTGGGACACAGCTTAACCCCGGTGGATgg                   >  1:226175/1‑68 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCTGGGACACAGCTTAACCCCGGTGGATgg                   >  1:516805/1‑68 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCTGGGACACAGCTTAACCCCGGTGGATgg                   >  1:1002280/1‑68 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCTGGGACACAGCTTAACCCCGGTGGATgg                   >  1:2092329/1‑68 (MQ=255)
                            tAAATTTTTCTCCGATCACCGTTTCACATAATATTCATGCTGGGACACAGCTTAACCCCGGTGGATgg                   >  1:3242945/1‑68 (MQ=255)
                               aTTTTTCTCCGATCACCGTTTCACATAATATTCATGCt                                              <  1:392103/38‑1 (MQ=38)
                                tttttCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGatat              <  1:2241453/69‑1 (MQ=255)
                                tttttCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGatat              <  1:2245446/69‑1 (MQ=255)
                                            cACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGgcgc  >  1:971639/1‑69 (MQ=255)
                                            cACCGTTTCACACAATATCCAAGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGgcgc  >  1:1226410/1‑69 (MQ=255)
                                                              |                                                  
CTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGC  >  minE/1560987‑1561099

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