Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,390,896 A→C 41.1% intergenic (‑40/+24) rhtC ← / ← recQ threonine efflux system/ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,390,8960AC41.1% 28.1 / 91.0 90intergenic (‑40/+24)rhtC/recQthreonine efflux system/ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (35/18);  new base C (19/18);  total (54/36)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATGCGCACGAATCAGCGC  >  minE/2390849‑2390961
                                               |                                                                 
aTCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATgctgct                                              <  1:2732886/69‑1 (MQ=255)
aTCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATgctgct                                              <  1:3324295/69‑1 (MQ=255)
aTCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATgctgct                                              <  1:618180/69‑1 (MQ=255)
aTCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATgctgct                                              <  1:2904625/69‑1 (MQ=255)
 tCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGctgctg                                             >  1:484476/1‑69 (MQ=255)
 tCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGctgctg                                             >  1:361256/1‑69 (MQ=255)
 tCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGctgctg                                             >  1:1243244/1‑69 (MQ=255)
        acatTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActc                                      <  1:77795/69‑1 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACTGAGTCTTCATACTGGCACtt                                                          <  1:1409689/48‑1 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACtt                                                          <  1:652293/48‑1 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACtt                                                          <  1:132438/48‑1 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACtt                                                          <  1:2851151/48‑1 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1951914/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1926335/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1861147/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1827340/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:3351205/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:170873/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:2465081/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1048406/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:1407206/1‑69 (MQ=255)
         catTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActct                                     >  1:859364/1‑69 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1511002/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1124138/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1149018/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1511580/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:233302/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1524841/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:2676411/61‑1 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2768579/1‑69 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2839371/1‑69 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:1538002/1‑69 (MQ=255)
           tttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:422132/1‑69 (MQ=255)
            ttGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGac                                           <  1:1856912/60‑1 (MQ=255)
               aCTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTActc                                      <  1:85265/62‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACTGAGTCTTCATCCTGGCACttttt                                                       <  1:1497198/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:999840/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:308804/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:530360/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:3255711/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:3375435/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:3425846/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:49143/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:375274/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:1178287/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:1874672/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:2161465/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:2312201/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:1121636/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:2334265/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:2367108/43‑1 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2590349/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2488469/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2848778/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2393053/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:2946862/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:1955334/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:3098376/1‑63 (MQ=255)
                 tCGCGGGGGAAACGTTTACGGAGTCTTCATCCTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:1025290/1‑63 (MQ=255)
                  cgcgGGGGAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:977008/42‑1 (MQ=255)
                      ggggAAACGTTTACGGAGTCTTCATCCTGGCACttttt                                                       <  1:2986153/38‑1 (MQ=255)
                          aaaCGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGact                                          >  1:1331710/1‑47 (MQ=255)
                          aaaCGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAAc                    >  1:2441281/1‑69 (MQ=255)
                          aaaCGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAAc                    >  1:1544351/1‑69 (MQ=255)
                          aaaCGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAAc                    >  1:170863/1‑69 (MQ=255)
                          aaaCGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAAc                    >  1:3019763/1‑69 (MQ=255)
                               tttACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcg               >  1:2432355/1‑69 (MQ=255)
                               tttACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcg               >  1:1662139/1‑69 (MQ=255)
                               tttACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcg               >  1:943628/1‑69 (MQ=255)
                               tttACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcg               >  1:1530573/1‑69 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATgctgct                                              >  1:976311/1‑36 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTc                                   >  1:1790314/1‑47 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:2116677/1‑68 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:2704554/1‑68 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:1100287/1‑68 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:579437/1‑68 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:3024404/1‑68 (MQ=255)
                                 tACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATgcgc              >  1:1805251/1‑68 (MQ=255)
                                       gTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATGCGCACGAATc       >  1:1652462/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:608335/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:1852169/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:667049/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:2734851/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:834595/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:872772/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcgc  >  1:1776226/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcg   >  1:2216928/1‑68 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACATGCGCACGAATCAgcg   >  1:1354551/1‑68 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATCGACAAGCGCACGAATCAgcgc  >  1:1054163/1‑69 (MQ=255)
                                            cATCCTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCGTCACCATAGACATGCGCACGAATCAgcgc  >  1:857258/1‑69 (MQ=255)
                                               |                                                                 
ATCAACATACATTTGACTCGCGGGGGAAACGTTTACGGAGTCTTCATACTGGCACTTTTTTATGCTGCTGACTACTCTTCGTCATCGCCATCAACATGCGCACGAATCAGCGC  >  minE/2390849‑2390961

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