Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,258,086 G→A 8.4% D1019D (GAC→GAT rpoB ← RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,258,0860GA 8.4% 221.2 / 8.1 71D1019D (GAC→GATrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CAGATCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTCCAGCC  >  minE/2258017‑2258153
                                                                     |                                                                   
cAGATCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGt                                                                    <  1:3035439/71‑1 (MQ=255)
cAGATCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGt                                                                    <  1:321552/71‑1 (MQ=255)
 aGATCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGt                                                                    >  1:1690847/1‑70 (MQ=255)
   aTCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCa                                                                  <  1:3297650/70‑1 (MQ=255)
   aTCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCAt                                                                 <  1:3065970/71‑1 (MQ=255)
    tcgtcgCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATa                                                                >  1:1152802/1‑71 (MQ=255)
      gtcgCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACt                                                              <  1:2274815/71‑1 (MQ=255)
      gtcgCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACt                                                              <  1:1922313/71‑1 (MQ=255)
        cgCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATAGTGc                                                            >  1:3100279/1‑71 (MQ=255)
         gCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCt                                                           <  1:1359600/71‑1 (MQ=255)
         gCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCt                                                           >  1:3399572/1‑71 (MQ=255)
            cTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTAGAACTCGTGTTTCAGTTCGTCATACTGCTcag                                                        >  1:3031702/1‑71 (MQ=255)
              gggTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCa                                                                  <  1:1737439/59‑1 (MQ=255)
               ggTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTcagc                                                       <  1:354859/69‑1 (MQ=255)
               ggTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTcagc                                                       >  1:3496858/1‑69 (MQ=255)
                 tGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCa                                                                  <  1:2719764/56‑1 (MQ=255)
                 tGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATAc                                                               <  1:911616/59‑1 (MQ=255)
                   aTTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTg                                                 >  1:1816217/1‑71 (MQ=255)
                   aTTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTg                                                 >  1:1750844/1‑71 (MQ=255)
                   aTTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTg                                                 >  1:1449103/1‑71 (MQ=255)
                    ttttGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGt                                                <  1:4096855/71‑1 (MQ=255)
                    ttttGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGt                                                <  1:710741/71‑1 (MQ=255)
                      ttGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTc                                              >  1:2286996/1‑71 (MQ=255)
                           gcgTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCt                                         <  1:1251805/71‑1 (MQ=255)
                            cgTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGtt                                               >  1:2646018/1‑64 (MQ=255)
                                tCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAttt                                    >  1:3786300/1‑71 (MQ=255)
                                tCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAttt                                    >  1:834087/1‑71 (MQ=255)
                                tCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAttt                                    >  1:2730069/1‑71 (MQ=255)
                                tCGCTTCGAGTTTCTTCTCGAACTAGTGTTTCAGTTCGt                                                                    >  1:2628206/1‑39 (MQ=255)
                                 cGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAtttt                                   >  1:838442/1‑71 (MQ=255)
                                 cGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAttt                                    >  1:695702/1‑70 (MQ=255)
                                 cGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAttt                                    >  1:3546127/1‑70 (MQ=255)
                                  gCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTg                                  <  1:3501211/71‑1 (MQ=255)
                                  gCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTg                                  >  1:331770/1‑71 (MQ=255)
                                   cTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGt                                 <  1:1077415/71‑1 (MQ=255)
                                   cTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGt                                 <  1:368351/71‑1 (MQ=255)
                                      cGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGttt                               <  1:3129692/70‑1 (MQ=255)
                                      cGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtc                              >  1:1063812/1‑71 (MQ=255)
                                      cGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtc                              >  1:3035317/1‑71 (MQ=255)
                                        aGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtctc                            <  1:389079/71‑1 (MQ=255)
                                         gTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtctgt                           >  1:2063901/1‑69 (MQ=255)
                                         gTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTctct                           <  1:2710819/71‑1 (MQ=255)
                                           ttcttcTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTc                         >  1:2334901/1‑71 (MQ=255)
                                           ttcttcTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTc                         >  1:2329930/1‑71 (MQ=255)
                                           ttcttcTCGAACTCGTGTTTCAGTTCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTc                         >  1:2797421/1‑71 (MQ=255)
                                            tcttctCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGtt                                               >  1:879234/1‑48 (MQ=255)
                                            tcttctCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGtt                                               >  1:1498163/1‑48 (MQ=255)
                                            tcttctCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGtt                                               >  1:1509749/1‑48 (MQ=255)
                                              ttctCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCt                                         >  1:3734248/1‑52 (MQ=255)
                                                ctcGAACTCGTGTTTAAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGt                       <  1:2414323/68‑1 (MQ=255)
                                                   gAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtct                             >  1:2781414/1‑59 (MQ=255)
                                                   gAACTCGTGTTTCAGTTCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGt                   <  1:2451996/69‑1 (MQ=255)
                                                    aaCTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAg                <  1:3311348/71‑1 (MQ=255)
                                                    aaCTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAg                <  1:1124213/71‑1 (MQ=255)
                                                      cTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTg                                  <  1:3688239/51‑1 (MQ=255)
                                                      cTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtctc                            >  1:1670435/1‑57 (MQ=255)
                                                      cTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTtctc                            >  1:366249/1‑57 (MQ=255)
                                                         gtgtTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGccc            <  1:2918352/70‑1 (MQ=255)
                                                           gtTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAtttt                                   <  1:1934807/45‑1 (MQ=255)
                                                           gtTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAtttt                                   <  1:4054079/45‑1 (MQ=255)
                                                           gtTTCAGTGCGTCATACTGCTCAGCCAGCTGTTCCAGCTGAtttt                                   <  1:32427/45‑1 (MQ=255)
                                                                 gTTCGTCATACTGCTCCGCCAGCTGTTCCAGCTGAt                                      <  1:3880486/36‑1 (MQ=255)
                                                                  ttCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:3850814/71‑1 (MQ=255)
                                                                  ttCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:1338804/71‑1 (MQ=255)
                                                                  ttCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:2712584/71‑1 (MQ=255)
                                                                  ttCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:3738504/71‑1 (MQ=255)
                                                                  ttCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:2057741/71‑1 (MQ=255)
                                                                  ttCATCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagcc  <  1:3338619/71‑1 (MQ=255)
                                                                   tCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTccagc   >  1:343946/1‑69 (MQ=255)
                                                                    cGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGttt                               >  1:1771984/1‑40 (MQ=255)
                                                                    cGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGt                   >  1:3214732/1‑52 (MQ=255)
                                                                     gTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTc                         <  1:714123/45‑1 (MQ=255)
                                                                     |                                                                   
CAGATCGTCGCCCTGGGTGATTTTGCGGCGTTTCGCTTCGAGTTTCTTCTCGAACTCGTGTTTCAGTTCGTCATACTGCTCAGCCAGCTGTTCCAGCTGATTTTGTTTCTCTTCGTCTGTCAGGCCCAGCTCCAGCC  >  minE/2258017‑2258153

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