Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 346,983 G→A 6.7% P210L (CCG→CTG)  ybbO ← predicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE346,9830GA 6.7% 171.8 / 2.1 59P210L (CCG→CTG) ybbOpredicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base A (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGACC  >  minE/346919‑347051
                                                                |                                                                    
taaAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATTCCGGGAttt                                                                <  1:5110085/71‑1 (MQ=255)
  aaaGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCAGGATTTTc                                                              <  1:3253168/71‑1 (MQ=255)
  aaaGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCAGGATTTTc                                                              <  1:2764610/71‑1 (MQ=255)
       aTGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTg                                                         >  1:856812/1‑71 (MQ=255)
       aTGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTg                                                         >  1:1547835/1‑71 (MQ=255)
           cGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGttt                                                     <  1:2211133/71‑1 (MQ=255)
              aCTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTAt                                                   >  1:825170/1‑70 (MQ=255)
               cTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCa                                                 >  1:5245117/1‑71 (MQ=255)
                 ttGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACt                                               >  1:5165763/1‑71 (MQ=255)
                    tCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACttt                                             <  1:5533512/70‑1 (MQ=255)
                      caccacCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGttt                                                     <  1:4957382/60‑1 (MQ=255)
                      caccacCGCTTCCGGTCCCAACGTAAAGCGG‑‑GGCGATGCCGGGATTTTCGACTGGttt                                                     <  1:3085279/58‑1 (MQ=255)
                       accaccGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGt                                         <  1:5480999/71‑1 (MQ=255)
                         caccGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTa                                                    >  1:953029/1‑58 (MQ=255)
                         caccGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTc                                        >  1:1247378/1‑70 (MQ=255)
                          accGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGGTTTTCGACTGGTTTATCACTTTGCGTCt                                       >  1:107586/1‑70 (MQ=255)
                          accGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTg                                      >  1:4491382/1‑71 (MQ=255)
                          accGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTg                                      >  1:4345130/1‑71 (MQ=255)
                             gCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGtt                                   <  1:702171/71‑1 (MQ=255)
                              cTTCCGGTCCCAACGTAAAGCGGGCGGCGACGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGtt                                   >  1:1614627/1‑70 (MQ=255)
                               ttCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCg                                                             <  1:1676218/43‑1 (MQ=255)
                                tCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTa                                                    <  1:238480/51‑1 (MQ=255)
                                  cGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACg                               <  1:2670684/70‑1 (MQ=255)
                                     tCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGt                           >  1:3689610/1‑71 (MQ=255)
                                       ccAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCa                                                 <  1:2556918/47‑1 (MQ=255)
                                       ccAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCa                                                 <  1:463969/47‑1 (MQ=255)
                                       ccAACGTAAAGCGGGGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCgg                        >  1:1657934/1‑71 (MQ=255)
                                        cAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTc                          <  1:1354723/69‑1 (MQ=255)
                                         aaCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGt                                         >  1:3970366/1‑53 (MQ=255)
                                         aaCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTCCGt                                         >  1:1182406/1‑53 (MQ=255)
                                          aCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTa                                                    <  1:4056097/41‑1 (MQ=255)
                                          aCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTa                                                    <  1:2715665/41‑1 (MQ=255)
                                          aCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATc                                                  <  1:4501424/43‑1 (MQ=255)
                                          aCGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCgg                        >  1:2013838/1‑69 (MQ=255)
                                           cttAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTa                                                    <  1:4319471/38‑1 (MQ=255)
                                            gTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGa                     <  1:4523240/70‑1 (MQ=255)
                                            gTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGa                     <  1:5166002/70‑1 (MQ=255)
                                              aaaGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGc                  >  1:3217670/1‑71 (MQ=255)
                                                  cgggcggCGATGCCAGGATTTTCGACTGGTTTATCACTTTGCGTCTGATTGACGTTGTCGGTGAAGCGa                >  1:1345336/1‑69 (MQ=255)
                                                  cgggcggCGATGCCAGGATTTTCGACTGGTTTATCACTTTGCGTCTGATTGACGTTGTCGGTGAAGCGAg               >  1:284547/1‑70 (MQ=255)
                                                   gggcggCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTa             >  1:5082243/1‑71 (MQ=255)
                                                      cggcgATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTg                                  <  1:228188/47‑1 (MQ=255)
                                                      cggcgATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGa          <  1:1515228/71‑1 (MQ=255)
                                                      cggcgATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGa          <  1:3795054/71‑1 (MQ=255)
                                                          gATGCCGGGATTTTCGACTGGTTTATCACTTTGCTTCTGGTTGACGTTGTCGGTGaa                    <  1:4722551/57‑1 (MQ=255)
                                                          gATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTc                          >  1:5659257/1‑51 (MQ=255)
                                                           aTGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATgg      >  1:3866249/1‑70 (MQ=255)
                                                           aTGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATgg      >  1:888177/1‑70 (MQ=255)
                                                            tGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGa    >  1:5301489/1‑71 (MQ=255)
                                                            tGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGa    >  1:5452601/1‑71 (MQ=255)
                                                            tGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGa    >  1:2201149/1‑71 (MQ=255)
                                                             gCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGtt                                   <  1:2292242/39‑1 (MQ=255)
                                                             gCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAg                   >  1:2425501/1‑55 (MQ=255)
                                                             gCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGAc   >  1:711306/1‑71 (MQ=255)
                                                             gCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGAc   >  1:897082/1‑71 (MQ=255)
                                                             gCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGAc   >  1:1126796/1‑71 (MQ=255)
                                                             gCCGGGATTTTCGACTGTTTATCACTTTGCGTCTGGtt                                   >  1:93808/1‑38 (MQ=38)
                                                              ccGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGAcc  >  1:5488389/1‑71 (MQ=255)
                                                              ccGGGATTTTCGACTGGTTTATCACTTAGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGa          >  1:4065448/1‑63 (MQ=255)
                                                                |                                                                    
TAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATCACTTTGCGTCTGGTTGACGTTGTCGGTGAAGCGAGTACGAATGGGACC  >  minE/346919‑347051

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: