Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,450,298 A→G W78R (TGG→CGG)  yifE ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,450,2980AG99.6% 884.3 / ‑4.9 228W78R (TGG→CGG) yifEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (110/117);  minor base T (1/0);  total (111/117)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 251 total aligned reads displayed.

CTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTC  >  minE/2450230‑2450364
                                                                    |                                                                  
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:2151364/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:2430492/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:4291582/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:4199564/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:3079328/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:3955183/1‑71 (MQ=255)
cTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcgca                                                                  >  1:328879/1‑71 (MQ=255)
 tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                 >  1:933082/1‑71 (MQ=255)
 tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                 >  1:2603854/1‑71 (MQ=255)
   cGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcg                                                                <  1:19403/70‑1 (MQ=255)
   cGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcg                                                                <  1:772259/70‑1 (MQ=255)
   cGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcg                                                                <  1:874553/70‑1 (MQ=255)
   cGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGCACCGCAcg                                                                <  1:2581289/70‑1 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             >  1:3727607/1‑71 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             <  1:4033884/71‑1 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             >  1:1001411/1‑71 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             >  1:685579/1‑71 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             >  1:2017060/1‑71 (MQ=255)
     gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                             <  1:2974384/71‑1 (MQ=255)
                gAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGt                                                  <  1:2495455/71‑1 (MQ=255)
                gAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGt                                                  <  1:2339550/71‑1 (MQ=255)
                  aaGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCa                                                >  1:1714021/1‑71 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:2996790/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:1765232/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:2010597/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:2501936/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:3097402/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:3225157/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:2781747/71‑1 (MQ=255)
                     ggTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:1422028/71‑1 (MQ=255)
                      gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTgg                                            >  1:1097202/1‑71 (MQ=255)
                      gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTgg                                            >  1:3624260/1‑71 (MQ=255)
                      gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTgg                                            >  1:4307584/1‑71 (MQ=255)
                      gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             >  1:208078/1‑70 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCt                                                       >  1:3548957/1‑59 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                           >  1:2363859/1‑71 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                           >  1:1200830/1‑71 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                           >  1:2686990/1‑71 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                           >  1:2678202/1‑71 (MQ=255)
                       tgtgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGc                                           >  1:819774/1‑71 (MQ=255)
                         tgAAAACGTTTTGGACGCTTAATTCGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTg                                             <  1:1566325/67‑1 (MQ=255)
                           aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCAc                                       >  1:2461266/1‑71 (MQ=255)
                           aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCAc                                       >  1:703993/1‑71 (MQ=255)
                             aaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGt                                     <  1:2280128/71‑1 (MQ=255)
                                 ttttGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 <  1:3669814/71‑1 (MQ=255)
                                 ttttGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 <  1:4271415/71‑1 (MQ=255)
                                 ttttGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 <  1:2862687/71‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:3322970/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:1801589/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:3804050/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:3283371/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:1036360/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:2544005/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:1604015/70‑1 (MQ=255)
                                     ggACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                              <  1:2313498/70‑1 (MQ=255)
                                        cGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                          <  1:1579005/71‑1 (MQ=255)
                                        cGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                          <  1:3655408/71‑1 (MQ=255)
                                        cGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                          <  1:2229099/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:3026347/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:3904819/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:4176663/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:3370793/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:2390805/71‑1 (MQ=255)
                                                  gcgTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACg                <  1:1168366/71‑1 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACt                                              >  1:3682075/1‑36 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:1342004/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:2292546/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:3754666/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:1195741/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:3059962/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:2631125/1‑58 (MQ=255)
                                                       atatACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATa                        >  1:3422933/1‑58 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                   <  1:2816387/44‑1 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                                   <  1:409239/44‑1 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt        >  1:1497118/1‑71 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt        >  1:2651355/1‑71 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt         >  1:3935818/1‑70 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt         >  1:4022947/1‑70 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt         >  1:2721568/1‑70 (MQ=255)
                                                          taCTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt         >  1:3492520/1‑70 (MQ=255)
                                                          taCTTGGACCGCAAGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt        >  1:1744438/1‑71 (MQ=255)
                                                            cTTGGACCTCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:1695849/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGtt                                    >  1:1466697/1‑41 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCg                                  >  1:854968/1‑43 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:354046/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:3216240/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:523819/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:4349745/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:1511125/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                                 >  1:3827752/1‑44 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtct      >  1:2527000/1‑71 (MQ=255)
                                                            cTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtct      >  1:4142643/1‑71 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc     <  1:3160654/69‑1 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc     <  1:16246/69‑1 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc     <  1:2059404/69‑1 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc     <  1:386739/69‑1 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc     <  1:4084459/69‑1 (MQ=255)
                                                               ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTctct    <  1:1903656/70‑1 (MQ=255)
                                                                gACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTc  <  1:4200813/71‑1 (MQ=255)
                                                                gACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTc  <  1:568620/71‑1 (MQ=255)
                                                                    |                                                                  
CTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTC  >  minE/2450230‑2450364

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: