Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,739 A→G 40.0% T452T (ACT→ACC pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,7390AG40.0% 28.8 / 51.7 75T452T (ACT→ACCpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (38/7);  new base G (14/16);  total (52/23)
Fisher's exact test for biased strand distribution p-value = 7.95e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

ATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACA  >  minE/223685‑223795
                                                      |                                                        
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGTCCCGATAGa                                            >  1:362569/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGc                                                    >  1:3232653/1‑61 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:1850641/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2605085/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2146257/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2500862/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2507784/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2628711/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2777490/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:2899595/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:3164543/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:3317888/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:3433735/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:53245/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:561076/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:637656/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:710974/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:926195/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:1085409/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:1146301/1‑69 (MQ=255)
aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                            >  1:1454564/1‑69 (MQ=255)
   aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         <  1:771991/69‑1 (MQ=255)
   aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         <  1:530628/69‑1 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:2822349/1‑69 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:1357405/1‑69 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:1377902/1‑69 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:1431329/1‑69 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:1534298/1‑69 (MQ=255)
      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                      >  1:2143720/1‑69 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGa                                                >  1:2671643/1‑53 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:1338379/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:402019/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:456705/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:505456/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:2909261/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:2870050/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:2805459/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:2626795/1‑60 (MQ=255)
            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         >  1:772333/1‑60 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:2115368/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1754786/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:340507/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1411383/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:448726/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:839464/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:2280958/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:2326173/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:2419012/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1664091/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1736025/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1202462/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1841450/67‑1 (MQ=255)
                            tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:580694/67‑1 (MQ=255)
                               tCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc                                         <  1:1876366/41‑1 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:1693694/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:187153/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:3037908/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:4209/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:58482/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:1892480/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:1591842/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:2692502/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:263461/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:2587173/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:227151/1‑68 (MQ=255)
                                ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:1674299/1‑68 (MQ=255)
                                ccGGAGACTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:1162406/1‑68 (MQ=255)
                                acGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg             >  1:2784648/2‑68 (MQ=255)
                                 cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGaa           <  1:2072798/69‑1 (MQ=255)
                                 cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGaa           <  1:1658458/69‑1 (MQ=255)
                                     gagTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:520597/58‑1 (MQ=255)
                                       gtgtGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt                  <  1:1011388/56‑1 (MQ=255)
                                        tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca    <  1:1905813/69‑1 (MQ=255)
                                        tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca    <  1:2951569/69‑1 (MQ=255)
                                          tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcaca  >  1:1655019/1‑69 (MQ=255)
                                                      |                                                        
ATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACA  >  minE/223685‑223795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: