Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 550,934 T→C 9.1% K4K (AAA→AAG ybiV ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE550,9340TC 9.1% 218.7 / 10.4 77K4K (AAA→AAGybiVpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/5);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

ATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGA  >  minE/550866‑550997
                                                                    |                                                               
aTAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACttt                                                               <  1:4236989/71‑1 (MQ=255)
aTAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACttt                                                               <  1:4415545/71‑1 (MQ=255)
aTAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACttt                                                               <  1:4845396/71‑1 (MQ=255)
 tAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACttt                                                               >  1:2973472/1‑70 (MQ=255)
   aaaCGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTa                                                             >  1:2230703/1‑70 (MQ=255)
   aaaCGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTAc                                                            >  1:2761889/1‑71 (MQ=255)
         ggttggttGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACttt                                                               <  1:1543378/62‑1 (MQ=255)
                   aCGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTc                                            >  1:128666/1‑71 (MQ=255)
                     gTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCt                                          >  1:3742891/1‑71 (MQ=255)
                     gTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCt                                          >  1:5190347/1‑71 (MQ=255)
                        ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTa                                                             >  1:4984247/1‑49 (MQ=255)
                        ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTa                                                             >  1:4632488/1‑49 (MQ=255)
                        ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTa                                                             >  1:4946738/1‑49 (MQ=255)
                        ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTa                                                             >  1:1879992/1‑49 (MQ=255)
                        ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGc                                       >  1:2404106/1‑71 (MQ=255)
                           gCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCt                                                         <  1:70340/50‑1 (MQ=255)
                           gCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCt                                                         <  1:3761266/50‑1 (MQ=255)
                           gCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCt                                                         <  1:2873868/50‑1 (MQ=255)
                           gCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCt                                                         <  1:765400/50‑1 (MQ=255)
                           gCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCt                                                         <  1:846474/50‑1 (MQ=255)
                             gtcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCTCa                                                       >  1:4047421/1‑50 (MQ=255)
                             gtcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACCTTTACGCTCa                                                       >  1:427499/1‑50 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:1514052/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:913198/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:1299078/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:2774294/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:2088142/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:2425023/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:2982195/53‑1 (MQ=255)
                              tcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:2977197/71‑1 (MQ=255)
                               cgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTTTGCTCCTGGCTCGaaaa                                >  1:107857/1‑71 (MQ=255)
                                gtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:4665260/51‑1 (MQ=255)
                                gtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAAt                               >  1:1747415/1‑71 (MQ=255)
                                 ttAATAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCAt                                                      <  1:683418/47‑1 (MQ=255)
                                 ttAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCAt                                                      <  1:4158328/47‑1 (MQ=255)
                                 ttAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCAt                                                      <  1:3923280/47‑1 (MQ=255)
                                 ttAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCAt                                                      <  1:956968/47‑1 (MQ=255)
                                 ttAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCAt                                                      <  1:3674537/47‑1 (MQ=255)
                                   aaGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCa                                                       >  1:636195/1‑44 (MQ=255)
                                   aaGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCa                                                       >  1:395323/1‑44 (MQ=255)
                                           tACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAAc                    <  1:3820587/71‑1 (MQ=255)
                                             ccGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAg                  >  1:3973948/1‑71 (MQ=255)
                                             ccGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAg                  >  1:5163351/1‑71 (MQ=255)
                                             ccGTCCATGCCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAg                  >  1:3389272/1‑71 (MQ=255)
                                                tCCATGTCTGTGACGATAACTTTTACGCTCATAgg                                                   <  1:5172810/35‑1 (MQ=255)
                                                 ccATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTc                                            <  1:2632372/41‑1 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:3097496/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:5095454/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:2945563/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:2119032/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:2068874/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGaaa                           >  1:4586046/1‑57 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGtaa             >  1:81579/1‑71 (MQ=255)
                                                  cATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGtaa             >  1:46324/1‑71 (MQ=255)
                                                   aTGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAAc                    <  1:1171226/63‑1 (MQ=255)
                                                      tCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGtaataa          >  1:762443/1‑70 (MQ=255)
                                                      tCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGtaataa          >  1:3642527/1‑70 (MQ=255)
                                                      tCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAAc         <  1:2686864/71‑1 (MQ=255)
                                                      tCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAAc         <  1:4854070/71‑1 (MQ=255)
                                                      tCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAAc         <  1:2164194/71‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGa                             >  1:1508038/1‑47 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGa                             >  1:792570/1‑47 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:3584167/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:103121/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:5017026/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:231315/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:4723932/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:3167098/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:4619111/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAgtg                <  1:2252070/60‑1 (MQ=255)
                                                          tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAAtg     <  1:3288939/71‑1 (MQ=255)
                                                            aCGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAAtgt    <  1:2505390/70‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:2840765/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:1382643/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:300992/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:2891208/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:1827974/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGaaa                                 <  1:525170/38‑1 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGa  >  1:4179569/1‑69 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGa  >  1:2599813/1‑69 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGa  >  1:4693659/1‑69 (MQ=255)
                                                               aTAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGa  >  1:3599391/1‑69 (MQ=255)
                                                                    |                                                               
ATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGA  >  minE/550866‑550997

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: