Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,656 G→A 23.0% T34T (ACG→ACA ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,6560GA23.0% 167.4 / 46.0 87T34T (ACG→ACAydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base A (6/14);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

TCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACAC  >  minE/982587‑982725
                                                                     |                                                                     
tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                                                                      >  1:1360865/1‑71 (MQ=255)
 ctGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAtt                                                                     >  1:1404636/1‑71 (MQ=255)
  tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                                                                     >  1:311942/1‑70 (MQ=255)
  tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                                                                     >  1:386150/1‑70 (MQ=255)
    gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                                                                  <  1:2970778/71‑1 (MQ=255)
    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:238927/71‑1 (MQ=255)
    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:3844028/71‑1 (MQ=255)
         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                                                   >  1:3916346/1‑65 (MQ=255)
         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAcc                                                              <  1:3676827/70‑1 (MQ=255)
         cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                                                                  <  1:2735242/66‑1 (MQ=255)
          ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa                                                            <  1:41107/71‑1 (MQ=255)
           gCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGa                                                            >  1:2761011/1‑70 (MQ=255)
           gCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                                                                    <  1:1738321/62‑1 (MQ=255)
            cGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaa                                                                 <  1:23565/64‑1 (MQ=255)
             gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                         <  1:553774/71‑1 (MQ=255)
              tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa                                                                >  1:2926376/1‑63 (MQ=255)
              tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa                                                                >  1:1087032/1‑63 (MQ=255)
                tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                         <  1:1761545/68‑1 (MQ=255)
                tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                      <  1:851923/71‑1 (MQ=255)
                tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                      <  1:3158471/71‑1 (MQ=255)
                  tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTg                                                    >  1:2929104/1‑71 (MQ=255)
                   cccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                   >  1:1490063/4‑71 (MQ=255)
                   ccAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  >  1:2051826/1‑56 (MQ=255)
                    cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc                                                   <  1:3745069/70‑1 (MQ=255)
                     caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGaaa                                                          >  1:3023933/2‑62 (MQ=255)
                         cAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTg                                                       <  1:3667758/61‑1 (MQ=255)
                          aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCg                                             <  1:784787/70‑1 (MQ=255)
                          aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgc                                            <  1:52041/71‑1 (MQ=255)
                           ggATGGAACGTTCACTGGCATGACGTATACCAAAAGCAAAACATTTaaaa                                                                >  1:2531996/1‑50 (MQ=255)
                            gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                                                   <  1:572479/46‑1 (MQ=255)
                             aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGttt                                                                    <  1:671834/44‑1 (MQ=255)
                             aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg                                         <  1:2816653/71‑1 (MQ=255)
                                gCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                                           <  1:1522643/50‑1 (MQ=255)
                                gCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctggctg                                      >  1:1633841/1‑71 (MQ=255)
                                 aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                      <  1:3835257/69‑1 (MQ=255)
                                 aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                     <  1:3140303/70‑1 (MQ=255)
                                 aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                     <  1:2702312/70‑1 (MQ=255)
                                  aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                    <  1:570782/71‑1 (MQ=255)
                                   cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                   >  1:1324791/1‑55 (MQ=255)
                                   cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                   <  1:1362761/71‑2 (MQ=255)
                                   cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                   <  1:2415209/71‑1 (MQ=255)
                                   cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                    >  1:2869985/1‑70 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:707003/39‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:2297997/39‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:3811901/39‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                                                  <  1:3501450/39‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                  <  1:2718327/71‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                  <  1:2181387/71‑1 (MQ=255)
                                    gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                    <  1:1293730/69‑1 (MQ=255)
                                     ttCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                                           <  1:3314056/45‑1 (MQ=255)
                                     ttCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaag                                 <  1:519500/71‑1 (MQ=255)
                                      tCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaag                                 >  1:777994/1‑70 (MQ=255)
                                       cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                         >  1:27844/1‑45 (MQ=255)
                                       cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                >  1:3461625/1‑70 (MQ=255)
                                       cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                >  1:2709370/1‑70 (MQ=255)
                                       cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                <  1:2544152/70‑1 (MQ=255)
                                       cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                >  1:2005320/1‑70 (MQ=255)
                                       cACCGCAATGACGTATACAAAAAGCAAAACGTTTAAAAc                                                               >  1:2901478/1‑39 (MQ=255)
                                         ccGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                    <  1:2051769/64‑1 (MQ=255)
                                          cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg                                                             <  1:2518115/38‑1 (MQ=255)
                                          cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                            >  1:860743/1‑71 (MQ=255)
                                          cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                            >  1:385486/1‑71 (MQ=255)
                                            cAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                              >  1:54342/1‑51 (MQ=255)
                                                gACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                   >  1:2420514/1‑42 (MQ=38)
                                                    tataCCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggtt                                       <  1:2769810/50‑3 (MQ=255)
                                                    tataCCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                   <  1:3588852/54‑2 (MQ=255)
                                                     ataCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                              >  1:2049951/1‑58 (MQ=255)
                                                      taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                 <  1:2102060/70‑1 (MQ=255)
                                                       aCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                                           <  1:2539005/43‑1 (MQ=255)
                                                         cAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATg             <  1:3491763/71‑1 (MQ=255)
                                                          aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtga                     >  1:3530296/1‑62 (MQ=255)
                                                          aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATgg            <  1:815703/71‑1 (MQ=255)
                                                           aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa           >  1:880873/1‑71 (MQ=255)
                                                           aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa           >  1:2436610/1‑71 (MQ=255)
                                                           aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa           >  1:1128136/1‑71 (MQ=255)
                                                            aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                                           >  1:1036833/1‑38 (MQ=255)
                                                            aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt          <  1:543831/71‑1 (MQ=255)
                                                                 aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                              >  1:56866/1‑46 (MQ=255)
                                                                 aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                       <  1:1946963/53‑1 (MQ=255)
                                                                 aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGa     <  1:1398796/71‑1 (MQ=255)
                                                                 aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCACGTGAGTGATATGGATTTCGa     >  1:113236/1‑71 (MQ=255)
                                                                 aaaCATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                   <  1:3593249/41‑2 (MQ=12)
                                                                 aaaCATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                   <  1:2567329/41‑2 (MQ=12)
                                                                  aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGa     >  1:252705/1‑70 (MQ=255)
                                                                   aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                              >  1:3720713/1‑44 (MQ=255)
                                                                   aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                        >  1:977316/1‑50 (MQ=255)
                                                                    cGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac  <  1:1396034/71‑1 (MQ=255)
                                                                     |                                                                     
TCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACAC  >  minE/982587‑982725

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: