Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 876,300 G→A 27.9% G506S (GGT→AGT)  acnA → aconitate hydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE876,3000GA27.9% 69.2 / 28.0 43G506S (GGT→AGT) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base A (9/3);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 3.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGT  >  minE/876230‑876370
                                                                      |                                                                      
aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGTATACGGTTGTACCACCTGTATTg                                                                        <  1:2173299/71‑1 (MQ=255)
  cGATACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGt                                                                      <  1:2549353/71‑1 (MQ=255)
  cGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGt                                                                      <  1:191944/71‑1 (MQ=255)
    aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTaa                                                                    >  1:2188396/1‑71 (MQ=255)
      aCTGACACCGTATCTCGAGGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAAct                                                                  >  1:434586/1‑71 (MQ=255)
      aCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAAc                                                                   >  1:2020226/1‑70 (MQ=255)
          acacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAAc                                                                   <  1:894271/66‑1 (MQ=255)
          acacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTgg                                                              <  1:871491/71‑1 (MQ=255)
            acCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTaa                                                                    <  1:572833/63‑1 (MQ=255)
            acCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTaa                                                                    <  1:2409073/63‑1 (MQ=255)
             ccGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGcc                                                           <  1:2378549/71‑1 (MQ=255)
               gTATCTCGACGAACTGGGGTTTAACCTTGTGGGATATGGTTGTACCACCTGTATTGGTAACTCTGGGCCGc                                                         >  1:3444381/1‑71 (MQ=255)
               gTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGc                                                         >  1:726100/1‑71 (MQ=255)
               gTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGc                                                         >  1:2441125/1‑71 (MQ=255)
                tATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGt                                                                      <  1:1533838/57‑1 (MQ=255)
                 aTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGCTg                                                       >  1:2947308/1‑71 (MQ=255)
                       acgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGAt                                                 <  1:3260409/71‑1 (MQ=255)
                         gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGCTGCCCGATcc                                               >  1:1676332/1‑71 (MQ=255)
                          aaCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt                                              <  1:900285/71‑1 (MQ=255)
                            cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAt                                            >  1:3344077/1‑71 (MQ=255)
                            cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAt                                            >  1:2915583/1‑71 (MQ=255)
                            cTGAGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAt                                            >  1:3441833/1‑71 (MQ=255)
                                  tttAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTaa                                                                    <  1:2351133/41‑1 (MQ=255)
                                    tAACCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGc                                                         >  1:1232921/1‑50 (MQ=255)
                                     aaCCTTGTGGGATACTGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGccc                                                    >  1:3697175/1‑54 (MQ=255)
                                     aaCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGccc                                                    >  1:3686278/1‑54 (MQ=255)
                                      aCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAActc                                                                 <  1:2366772/40‑1 (MQ=255)
                                      aCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAActc                                                                 <  1:2090246/40‑1 (MQ=255)
                                      aCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAActc                                                                 <  1:1507639/40‑1 (MQ=255)
                                      aCCTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCg                                                          >  1:4038017/1‑47 (MQ=255)
                                        cTTGTGGGATACGGTTGTACCACCTGTATTAGTAACTCTGGGCCGCTGCCCGATCCTATCGa                                         <  1:796982/62‑1 (MQ=255)
                                          tgtgGGATACGGTTGTACCACCTGTATTGGTAActc                                                                 <  1:1512539/36‑1 (MQ=255)
                                           gtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc                                                      <  1:929234/46‑1 (MQ=255)
                                           gtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc                                                      <  1:239217/46‑1 (MQ=255)
                                           gtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATc                                                >  1:2214050/1‑52 (MQ=255)
                                              ggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCaaaa                           <  1:3144306/70‑1 (MQ=255)
                                               gATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCaaaaaa                         <  1:1461358/71‑1 (MQ=255)
                                                   cGGTTGTACCACCTGTATTAGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                     <  1:754867/71‑1 (MQ=255)
                                                       tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAtt                  >  1:2486539/1‑70 (MQ=255)
                                                       tGTACCACCTGTATTAGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAtt                  >  1:3725698/1‑70 (MQ=255)
                                                         tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTaa               <  1:3586867/71‑1 (MQ=255)
                                                         tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                <  1:1748563/70‑1 (MQ=255)
                                                         tACCACCTGTATTAGTAACTCTGGGCCGCTGCCCGATCCTATCGaaa                                       >  1:813560/1‑47 (MQ=255)
                                                                    ttGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                >  1:2538541/1‑59 (MQ=255)
                                                                      ggTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcggt  >  1:3570779/1‑71 (MQ=255)
                                                                      |                                                                      
AGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGT  >  minE/876230‑876370

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: