Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 228,538 A→G 8.2% S55S (AGT→AGC phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE228,5380AG 8.2% 231.7 / 10.0 85S55S (AGT→AGCphoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base A (49/29);  new base G (4/3);  total (53/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCC  >  minE/228468‑228607
                                                                      |                                                                     
cTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATaa                                                                       >  1:3150205/1‑71 (MQ=255)
 ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                                                      >  1:425616/1‑71 (MQ=255)
 ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                                                      >  1:3209535/1‑71 (MQ=255)
   cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAActc                                                                    >  1:3380070/1‑71 (MQ=255)
   cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                                                      <  1:551892/69‑1 (MQ=255)
      aaCGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTgg                                                                 <  1:314234/71‑1 (MQ=255)
         gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGt                                                                >  1:1422945/1‑69 (MQ=255)
         gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                              >  1:157693/1‑71 (MQ=255)
         gACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTc                                                               >  1:3404868/1‑70 (MQ=255)
         gACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTc                                                               >  1:2812419/1‑70 (MQ=255)
         gACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCg                                                              >  1:844957/1‑71 (MQ=255)
          aCCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGc                                                             <  1:3226807/71‑1 (MQ=255)
             aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGACACCAAAACGGATATAACTCTGGTCGCCAt                                                          <  1:492455/71‑1 (MQ=255)
             aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCAt                                                          <  1:1180725/71‑1 (MQ=255)
             aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGACTTTGAAACCAAAACGGATATAACTCTGGTCGCCAt                                                          >  1:1137391/1‑71 (MQ=255)
              tagccAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                                                         <  1:2273121/68‑1 (MQ=12)
              tagccAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                                                         >  1:1835283/4‑71 (MQ=12)
              tagccAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                                                         >  1:289367/4‑71 (MQ=12)
                aCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATaa                                                                       >  1:722913/1‑55 (MQ=255)
                aCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCtt                                                       <  1:963391/71‑1 (MQ=255)
                aCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCtt                                                       <  1:2628137/71‑1 (MQ=255)
                 ccAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCttt                                                      <  1:166367/71‑1 (MQ=255)
                  cagtcagtTGATCGTTAATTTGTGTTTCGCGTTTGAAACCAAAACGGATATAACTCTGGTc                                                               >  1:700112/1‑61 (MQ=255)
                   agtgagtTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCtttcg                                                    <  1:3266790/71‑3 (MQ=19)
                    gtgagtTGATCGTTAATTTGGGTTTCGCCTTTAAAAGCCAGACGCGCATAGCTCTGGTCGCCATCt                                                        <  1:3865825/63‑1 (MQ=255)
                         ttGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTc                                                               >  1:3925720/1‑54 (MQ=255)
                         ttGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGt                                              >  1:292341/1‑71 (MQ=255)
                         ttGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGt                                              >  1:3127936/1‑71 (MQ=255)
                            aTCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCtt                                                       <  1:300585/59‑1 (MQ=255)
                             tCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTAt                                           >  1:3878702/1‑70 (MQ=255)
                             tCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATc                                          <  1:3818685/71‑1 (MQ=255)
                             tCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATc                                          <  1:2265918/71‑1 (MQ=255)
                              cGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATc                                          >  1:155143/1‑70 (MQ=255)
                               gTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCa                                         >  1:1034926/1‑70 (MQ=255)
                               gTTAATTTGTGTTTCGCCTTTGAAAACAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCa                                         <  1:2401921/70‑1 (MQ=255)
                                ttAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCa                                                           >  1:694900/1‑51 (MQ=255)
                                ttAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCa                                                           >  1:2862300/1‑51 (MQ=255)
                                ttAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACt                                       <  1:1135909/71‑1 (MQ=255)
                                ttAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACt                                       >  1:43192/1‑71 (MQ=255)
                                 tAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACt                                       >  1:1290886/1‑70 (MQ=255)
                                  aaTTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTc                                      <  1:1224362/70‑1 (MQ=255)
                                    tttGTGTTTCGCCTTTGAAACCAAAACGGATATAATTCTGTTCg                                                              >  1:3542332/1‑44 (MQ=255)
                                    tttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTAc                                                    <  1:1970715/54‑1 (MQ=255)
                                    tttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTAc                                                    <  1:1862248/54‑1 (MQ=255)
                                     ttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCata                                   >  1:1125748/1‑70 (MQ=255)
                                     ttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCata                                   >  1:787080/1‑70 (MQ=255)
                                       gtgtTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGc                                                >  1:187798/1‑55 (MQ=255)
                                        tgtTTCGCCTTTGAAACCAAAACGGATATAACTCTg                                                                  <  1:2717931/36‑1 (MQ=255)
                                        tgtTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCAc                                        >  1:1557438/1‑62 (MQ=255)
                                          tttCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATg                              <  1:1654923/70‑1 (MQ=255)
                                            tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc                                                         >  1:1778927/1‑41 (MQ=255)
                                            tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc                                                         >  1:753096/1‑41 (MQ=255)
                                            tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTa                                                     >  1:1008030/1‑45 (MQ=255)
                                            tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCa                            <  1:1959772/70‑1 (MQ=255)
                                            tCGCATTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTc                                      >  1:2035740/1‑60 (MQ=255)
                                             cGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTAc                                                    >  1:3409579/1‑45 (MQ=255)
                                             cGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATg                          >  1:372582/1‑71 (MQ=255)
                                             cGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATATTTACTGGCGTTATCACTCATATAATGCATg                          >  1:3378649/1‑71 (MQ=255)
                                                cTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTAACTGGCGt                                              >  1:60679/1‑48 (MQ=255)
                                                 tttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTa                                            <  1:3373877/49‑1 (MQ=255)
                                                  ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTAt                                           >  1:3632668/1‑49 (MQ=255)
                                                  ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTAt                                           >  1:1628621/1‑49 (MQ=255)
                                                  ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt                      >  1:2704674/1‑70 (MQ=255)
                                                  ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt                      >  1:3665960/1‑70 (MQ=255)
                                                   tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt                     >  1:3320600/1‑70 (MQ=255)
                                                   tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt                     >  1:1626692/1‑70 (MQ=255)
                                                   tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt                      <  1:2217319/69‑1 (MQ=255)
                                                    gAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTa                    <  1:2628581/70‑1 (MQ=255)
                                                     aaaCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAAc                  >  1:993791/1‑71 (MQ=255)
                                                     aaaCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAAc                  <  1:3341268/71‑1 (MQ=255)
                                                      aaCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCata                                   <  1:3870080/53‑1 (MQ=255)
                                                      aaCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACt                 <  1:1959813/71‑1 (MQ=255)
                                                      aaCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACt                 <  1:3910593/71‑1 (MQ=255)
                                                       aCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACtt                >  1:2276480/1‑71 (MQ=255)
                                                         cAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTc                                      >  1:2790916/1‑47 (MQ=255)
                                                         cAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTg              >  1:1979329/1‑71 (MQ=255)
                                                            aaCGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCAtata                                 >  1:549648/1‑49 (MQ=255)
                                                            aaCGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGcc            >  1:1347099/1‑70 (MQ=255)
                                                            aaCGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGc             >  1:1118761/1‑69 (MQ=255)
                                                             aCGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCa           >  1:2272089/1‑70 (MQ=255)
                                                               ggATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCa           >  1:90862/1‑68 (MQ=255)
                                                                gATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGa       <  1:1826869/71‑1 (MQ=255)
                                                                  tataACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGa       <  1:1811353/69‑1 (MQ=255)
                                                                  tataACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACa     <  1:3520924/71‑1 (MQ=255)
                                                                   ataACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt                     >  1:2493854/1‑54 (MQ=255)
                                                                   ataACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt                     >  1:2824907/1‑54 (MQ=255)
                                                                   ataACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt                     >  1:3611570/1‑54 (MQ=255)
                                                                     aaCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATc   >  1:1997662/1‑70 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATcc  <  1:2778880/70‑1 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATc   >  1:2298119/1‑69 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATc   >  1:2916033/1‑69 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATc   >  1:3676129/1‑69 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATc   >  1:1612215/1‑69 (MQ=255)
                                                                      aCTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAAAtt                >  1:1309816/1‑56 (MQ=255)
                                                                      |                                                                     
CTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCC  >  minE/228468‑228607

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: