Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,861,362 A→G 22.2% Q150R (CAG→CGG)  miaA → delta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,861,3620AG22.2% 124.7 / 36.5 63Q150R (CAG→CGG) miaAdelta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase
Reads supporting (aligned to +/- strand):  ref base A (19/30);  new base G (3/11);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGG  >  minE/2861294‑2861430
                                                                    |                                                                    
ggAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCagg                                                                    <  1:3579808/71‑1 (MQ=255)
ggAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCagg                                                                    <  1:2083939/71‑1 (MQ=255)
ggAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCag                                                                     >  1:1716106/1‑70 (MQ=255)
    gTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGt                                                                <  1:1277356/71‑1 (MQ=255)
       cGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTGATTGCATCGTCAACTTCgggagg                                                                 >  1:1794981/1‑67 (MQ=255)
       cGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCaggagg                                                                 >  1:648609/1‑67 (MQ=255)
       cGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCaggagg                                                                 >  1:3649504/1‑67 (MQ=255)
       cGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCaggagg                                                                 >  1:3775806/1‑67 (MQ=255)
        gggCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGAt                                                            <  1:1353039/71‑1 (MQ=255)
        gggCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGAt                                                            <  1:2854110/71‑1 (MQ=255)
         ggCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATc                                                           <  1:630735/71‑1 (MQ=255)
           ccAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCg                                                         <  1:2323232/71‑1 (MQ=255)
                aTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGc                                                    >  1:837195/1‑71 (MQ=255)
                aTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGc                                                    <  1:691062/71‑1 (MQ=255)
                aTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGc                                                    <  1:600664/71‑1 (MQ=255)
                aTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGc                                                    <  1:454289/71‑1 (MQ=255)
                aTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGc                                                    <  1:3969568/71‑1 (MQ=255)
                 ttGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCg                                                   >  1:850812/1‑71 (MQ=255)
                  tGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCgg                                                  >  1:1591809/1‑71 (MQ=255)
                      cAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGgcag                                               >  1:914265/1‑70 (MQ=255)
                      cAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGgcag                                               >  1:103596/1‑70 (MQ=255)
                      cAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGgcagc                                              <  1:2722027/71‑1 (MQ=255)
                         cAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAg                                           >  1:3402350/1‑71 (MQ=255)
                           ggcggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGa                                         >  1:932752/1‑71 (MQ=255)
                            gcggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGAt                                        <  1:3415336/71‑1 (MQ=255)
                              ggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTc                                      <  1:3264685/71‑1 (MQ=255)
                              ggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTc                                      <  1:2161522/71‑1 (MQ=255)
                              ggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTc                                      <  1:1915406/71‑1 (MQ=255)
                              ggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTc                                      <  1:2223624/71‑1 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCagg                                                                    >  1:1013446/1‑38 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:3277016/1‑44 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:2390507/1‑44 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:2562661/1‑44 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:2994432/1‑44 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:2896221/1‑44 (MQ=255)
                                 agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAg                                                              >  1:1527162/1‑44 (MQ=255)
                                   agCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCgg                                                  <  1:3017098/54‑1 (MQ=255)
                                   agCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCgg                                                  <  1:251947/54‑1 (MQ=255)
                                   agCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCa                                 <  1:43106/71‑1 (MQ=255)
                                    gCAAGGTTGGGAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:1060735/71‑1 (MQ=255)
                                    gCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:672493/71‑1 (MQ=255)
                                    gCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:484660/71‑1 (MQ=255)
                                    gCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:2629316/71‑1 (MQ=255)
                                    gCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:2643953/71‑1 (MQ=255)
                                     cAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATc                                   >  1:2001740/1‑67 (MQ=255)
                                     cAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaaa                               >  1:1018449/1‑71 (MQ=255)
                                     cAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCaa                                <  1:4112246/70‑1 (MQ=255)
                                       aGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATg                             <  1:485645/71‑1 (MQ=255)
                                           tGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATc                          >  1:2726115/1‑70 (MQ=255)
                                              gAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCaa                                            <  1:2719824/49‑1 (MQ=255)
                                              gAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCaa                                            <  1:526457/49‑1 (MQ=255)
                                              gAGTCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCaa                                            <  1:2861901/49‑1 (MQ=255)
                                               aGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCcaca                      >  1:1725558/1‑70 (MQ=255)
                                                 tCATTGCATCGTCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAg                   >  1:3704579/1‑71 (MQ=255)
                                                 tCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCa                                 <  1:391864/57‑1 (MQ=255)
                                                 tCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACaaa                    >  1:891472/1‑70 (MQ=255)
                                                 tCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAg                   <  1:3538169/71‑1 (MQ=255)
                                                 tCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAg                   <  1:1769537/71‑1 (MQ=255)
                                                 tCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAg                   >  1:3088215/1‑71 (MQ=255)
                                                   aTTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAg                                           <  1:3500415/45‑1 (MQ=255)
                                                    ttGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCcac                       <  1:2262037/64‑1 (MQ=255)
                                                      gCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCttt              >  1:66218/1‑71 (MQ=255)
                                                       cATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTc             <  1:1382317/71‑1 (MQ=255)
                                                        aTCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTcc            <  1:3623262/71‑1 (MQ=255)
                                                         tCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACaa                     <  1:809820/61‑1 (MQ=255)
                                                            tCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCa                                 <  1:3302085/46‑1 (MQ=255)
                                                            tCAACTTCGGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCa                                 <  1:461734/46‑1 (MQ=255)
                                                                   cAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTgg  <  1:818426/70‑1 (MQ=255)
                                                                   cAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTgg  >  1:2560073/1‑70 (MQ=255)
                                                                   cAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTg   <  1:725245/69‑1 (MQ=255)
                                                                   cAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTg   <  1:3355850/69‑1 (MQ=255)
                                                                    |                                                                    
GGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGG  >  minE/2861294‑2861430

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: