Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,402,395 T→C 40.8% Q212Q (CAA→CAG yigA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,402,3950TC40.8% 34.3 / 67.7 93Q212Q (CAA→CAGyigAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (15/23);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGCGTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCATCGCTT  >  minE/2402328‑2402461
                                                                   |                                                                  
cgcgTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTAATGAAGTAACTGCGTTCCtt                                                                   >  1:1344312/1‑69 (MQ=255)
  cgTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGc                                                                 >  1:1139668/1‑69 (MQ=255)
  cgTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGc                                                                 >  1:388979/1‑69 (MQ=255)
  cgTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGc                                                                 >  1:1558144/1‑69 (MQ=255)
      cAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCtt                                                             >  1:2231019/1‑69 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1490004/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2371280/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1990617/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:250249/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2661029/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2784217/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2799197/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1804009/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:3000112/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1509107/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:3214307/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1063523/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:142211/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:930876/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:3425948/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1109014/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:1277398/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:3303437/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:656702/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:676312/63‑1 (MQ=255)
          ccAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:756500/63‑1 (MQ=255)
           cAACGCGCCAGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2256177/62‑1 (MQ=255)
                cTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:534779/69‑1 (MQ=255)
                cTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:2348189/69‑1 (MQ=255)
                cTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:2949501/69‑1 (MQ=255)
                cTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:154638/69‑1 (MQ=255)
                cTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:1471516/69‑1 (MQ=255)
                   cAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTg                                                   <  1:3120304/66‑1 (MQ=255)
                    aGAAGCTCCGGCAACATCAGTGCAATTTCATGAAGTAACTGCGTTCCCTGccc                                                               <  1:2521264/53‑1 (MQ=255)
                                    tCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGGCTGGTa                                >  1:996053/1‑68 (MQ=255)
                                      aGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAAc                             >  1:2901077/1‑69 (MQ=255)
                                      aGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAAc                             >  1:1645323/1‑69 (MQ=255)
                                        cgcAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:795040/1‑69 (MQ=255)
                                        cgcAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:314355/1‑69 (MQ=255)
                                        cgcAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:445965/1‑69 (MQ=255)
                                        cgcAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:450071/1‑69 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCGTCGCGACTGGTAAACAg                           >  1:2251982/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:540843/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:1223965/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:686334/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:254281/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:2611121/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:40606/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:1472575/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:330088/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:2709333/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:2759328/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:2149421/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:1590231/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:1059849/1‑68 (MQ=255)
                                         gcAATTTCATGAAGTAACTGCGTTCCCTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAg                           >  1:3000111/1‑68 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGTTAAACAGCACGa                      <  1:2166656/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:3085014/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:1974331/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:311932/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:401976/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:3154642/69‑1 (MQ=255)
                                             tttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:3196935/69‑1 (MQ=255)
                                              ttCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGa                      <  1:2694219/68‑1 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:1661600/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:2275202/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:994550/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:2301687/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:2119330/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:2649587/1‑69 (MQ=255)
                                               tCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGaca                    >  1:3256666/1‑69 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTTGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:1970506/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTa                                >  1:1485201/1‑47 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:496715/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:293244/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:296444/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:73192/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:1161321/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:2398886/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:2458002/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:2796373/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAg          <  1:1389181/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACCACACCCAAATCAg          <  1:3118693/69‑1 (MQ=255)
                                                         aCTGCGTTCCTTGCCATTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCAAAATCAg          <  1:1930220/69‑1 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGc         >  1:2799934/1‑68 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:1939164/1‑69 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:1273166/1‑69 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:114406/1‑69 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:1086048/1‑69 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:1965809/1‑69 (MQ=255)
                                                           tGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCa        >  1:314754/1‑69 (MQ=255)
                                                             cGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCATc      >  1:1338002/1‑69 (MQ=255)
                                                                  cTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCATCGCtt  >  1:3376611/1‑68 (MQ=255)
                                                                  cTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCATCGCtt  >  1:3372514/1‑68 (MQ=255)
                                                                   |                                                                  
CGCGTTCAATCCAACGCTCCAGAAGCTCCGGCAACATCAGCGCAATTTCATGAAGTAACTGCGTTCCTTGCCCTTGTTGATAGTGACTGGCATCGCGACTGGTAAACAGCACGACACCCAAATCAGCATCGCTT  >  minE/2402328‑2402461

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