Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,172,453 C→A 56.6% L82L (CTC→CTA yhbJ → hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,172,4530CA56.6% 20.8 / 144.4 129L82L (CTC→CTAyhbJhypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (19/37);  new base A (28/45);  total (47/82)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 131 total aligned reads displayed.

TGCCGGAGTCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCA  >  minE/2172386‑2172515
                                                                   |                                                              
tGCCGGAGTCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTcc                                                               >  1:3391521/1‑69 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:1032879/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:950259/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:930540/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:2946107/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:2938666/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:2890947/69‑1 (MQ=255)
        tCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:2672732/69‑1 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTACCCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:1200709/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCTACTCCTGTTCCTGGa                                                     >  1:503345/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       >  1:2254918/1‑67 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:316028/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:167312/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:2036322/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:3362294/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:2550860/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:2730975/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:1288611/1‑69 (MQ=255)
          aCCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGa                                                     >  1:3291044/1‑69 (MQ=255)
           ccAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTg                                                       <  1:653442/66‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:2769795/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:2425724/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:3327958/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:3021811/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:1140551/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:584787/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:1224063/69‑1 (MQ=255)
                    ttCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:820777/69‑1 (MQ=255)
                      cGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:2252823/69‑1 (MQ=255)
                      cGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:1440626/69‑1 (MQ=255)
                      cGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:1252118/69‑1 (MQ=255)
                      cGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:2916099/69‑1 (MQ=255)
                       gAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:678214/68‑1 (MQ=255)
                             gCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTa                                  <  1:2450935/69‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTCATACCt                                    <  1:436355/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2675920/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:1107651/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2989759/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:3348376/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:826671/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:1070123/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:1309547/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2824986/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2521780/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:3215476/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:295008/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:1594336/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2243956/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2090794/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:2074056/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:1637444/65‑1 (MQ=255)
                               gATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:3287623/65‑1 (MQ=255)
                                 tGAGTAACCTGCCTGACGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCt                                    <  1:933427/63‑1 (MQ=255)
                                     tAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATtcgtc                           >  1:3242469/1‑68 (MQ=255)
                                        cctgcctgACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGt                                           <  1:2660238/49‑1 (MQ=255)
                                        cctgcctgACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTaa                                         <  1:3037956/51‑1 (MQ=255)
                                              tgACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGa                 <  1:2274100/69‑1 (MQ=255)
                                                aCGTTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:2527101/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:3172903/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:669670/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1189093/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:3273865/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:711517/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1038673/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:352602/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:509962/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:572165/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:3031269/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1277776/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1510406/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1742080/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:175107/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1874416/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1887305/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1932190/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:20098/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:1070298/1‑51 (MQ=255)
                                                aCGCTTTCTCACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTa                                  >  1:2282108/1‑50 (MQ=255)
                                                aCGCTTTCTCACCGCAACTAATGTTCCTGGATGCCGACCGTAATACCTTaa                                 >  1:371506/1‑51 (MQ=255)
                                                 cGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCAACCGTAATACCTTAATTCGTCGTTACAGTGacac              >  1:2526688/1‑69 (MQ=255)
                                                       ctcAGCGCAACTACTGTTCCTGGATGCCGAGCGTAATACCTTAATtc                              <  1:2699240/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:2314895/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:102904/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:3146972/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:1036795/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:1368826/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:3219345/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:70602/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:2989629/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:2974375/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:643128/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:2953955/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:1051709/47‑1 (MQ=255)
                                                       ctcACCGCAACTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:2301200/47‑1 (MQ=255)
                                                        tcACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGcc         <  1:732997/67‑1 (MQ=255)
                                                        tcACCGCAACTCCTGTTCCTGGATGACGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGcc         <  1:1331218/67‑1 (MQ=255)
                                                             gCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGc   >  1:873062/1‑68 (MQ=255)
                                                             gCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGc   >  1:393768/1‑68 (MQ=255)
                                                             gCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCa  >  1:562225/1‑69 (MQ=255)
                                                               aaCTACTGTTCCTGGATGCCGACCGTAATACCTTAATtc                              <  1:1787723/39‑1 (MQ=255)
                                                                   |                                                              
TGCCGGAGTCACCAGAAATATTCGAACAGGCGATGAGTAACCTGCCTGACGCTTTCTCACCGCAACTCCTGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCA  >  minE/2172386‑2172515

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