Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,294,050 A→C 18.6% intergenic (+77/‑198) gldA → / → yijF glycerol dehydrogenase, NAD/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,294,0500AC18.6% 121.6 / 38.0 70intergenic (+77/‑198)gldA/yijFglycerol dehydrogenase, NAD/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base C (5/8);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCAC  >  minE/2293988‑2294118
                                                              |                                                                    
ccGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGc                                                              <  1:295183/71‑1 (MQ=255)
 cGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCa                                                             <  1:1818989/71‑1 (MQ=255)
    cTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAata                                                          <  1:2046117/71‑1 (MQ=255)
          cGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGa                                                    <  1:3960317/71‑1 (MQ=255)
           gggAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCaa                                                            <  1:4085014/62‑1 (MQ=255)
              aGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                                <  1:601405/71‑1 (MQ=255)
              aGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                                <  1:4037149/71‑1 (MQ=255)
                  tGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                                                        >  1:1004926/1‑59 (MQ=255)
                  tGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGt                                            >  1:2806764/1‑71 (MQ=255)
                  tGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGt                                            >  1:2664707/1‑71 (MQ=255)
                   gAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtc                                                                   <  1:2129346/47‑1 (MQ=255)
                    accGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCaa                                                            <  1:1663005/51‑1 (MQ=14)
                    aaCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                                                        <  1:763694/57‑1 (MQ=255)
                    aaCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                                                        <  1:649066/57‑1 (MQ=255)
                      cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCtt                                        <  1:3655105/71‑1 (MQ=255)
                       gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCttt                                       >  1:2404164/1‑71 (MQ=255)
                        cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCtttt                                      >  1:944550/1‑71 (MQ=255)
                          aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                     >  1:2957537/1‑70 (MQ=255)
                            aTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                     >  1:3233528/1‑68 (MQ=255)
                            aTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAg                                  >  1:1762764/1‑71 (MQ=255)
                             ttGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAg                                  <  1:1090479/70‑1 (MQ=255)
                             atGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCt                                         <  1:1500061/62‑1 (MQ=255)
                                 cTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCttt                                       <  1:3973634/61‑1 (MQ=255)
                                  tGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCt                                         <  1:2019713/58‑1 (MQ=255)
                                  tGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                                <  1:3700046/51‑1 (MQ=255)
                                  tGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                                <  1:2636353/51‑1 (MQ=255)
                                   gATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                                                        <  1:3025520/42‑1 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCttt                                       >  1:2503167/1‑58 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                           <  1:714379/70‑1 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                           <  1:3630694/70‑1 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                           <  1:1186774/70‑1 (MQ=255)
                                    aTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                           <  1:2068245/70‑1 (MQ=255)
                                     tGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                   >  1:3923444/1‑61 (MQ=255)
                                     tGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTCTTGTa                                   >  1:3955572/1‑61 (MQ=255)
                                      gcgcTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAg                        >  1:1506268/1‑71 (MQ=255)
                                       cgcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCg                                             >  1:3083047/1‑49 (MQ=255)
                                       cgcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                   >  1:1513172/1‑59 (MQ=255)
                                       cgcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAg                        <  1:2085056/70‑1 (MQ=255)
                                       cgcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAg                        <  1:3893042/70‑1 (MQ=255)
                                        gcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGc                      <  1:3695798/71‑1 (MQ=255)
                                         cTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCg                     >  1:1128869/1‑71 (MQ=255)
                                          tACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGc                                          <  1:3171633/49‑1 (MQ=255)
                                           aCGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGcag                              <  1:3758293/60‑3 (MQ=255)
                                           aCGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                   >  1:2967751/1‑71 (MQ=255)
                                           aCGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                   >  1:634413/1‑71 (MQ=255)
                                           aCGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                   >  1:2012027/1‑71 (MQ=255)
                                            cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                   <  1:2181348/70‑1 (MQ=255)
                                            cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAACGCGtt                   <  1:3806914/70‑1 (MQ=255)
                                              cTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                     <  1:3015579/50‑1 (MQ=255)
                                              cTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCAc                >  1:3332386/1‑71 (MQ=255)
                                                tATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                   >  1:2796250/1‑50 (MQ=255)
                                                 aTCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTcc             >  1:1556165/1‑71 (MQ=255)
                                                  tCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCg            <  1:3480846/71‑1 (MQ=14)
                                                   cAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCa                 >  1:2325473/1‑65 (MQ=255)
                                                   cAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                   >  1:657776/1‑47 (MQ=255)
                                                     ggCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCAt         <  1:2421985/71‑1 (MQ=255)
                                                        cTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATa                          >  1:2168502/1‑51 (MQ=255)
                                                        cTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCg      >  1:3070385/1‑71 (MQ=255)
                                                         tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                           <  1:192892/49‑1 (MQ=255)
                                                         tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAgg                       <  1:3873079/53‑1 (MQ=255)
                                                         tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAgg                       <  1:2592740/53‑1 (MQ=255)
                                                           cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                     >  1:567249/1‑37 (MQ=255)
                                                           cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                         <  1:1927264/49‑1 (MQ=14)
                                                           cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                           <  1:191330/47‑1 (MQ=255)
                                                           cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGc                      >  1:1731883/1‑52 (MQ=255)
                                                           cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCgg     <  1:2956246/69‑1 (MQ=255)
                                                           cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCa   >  1:3182553/1‑71 (MQ=255)
                                                           cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCa   >  1:1297534/1‑71 (MQ=255)
                                                             cAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                   >  1:655763/1‑37 (MQ=255)
                                                             cAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCAc  >  1:1030857/1‑70 (MQ=255)
                                                              |                                                                    
CCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCAC  >  minE/2293988‑2294118

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: