Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 228,508 T→G 8.0% T65T (ACA→ACC phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE228,5080TG 8.0% 222.5 / 13.6 88T65T (ACA→ACCphoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base G (4/3);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

ACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATA  >  minE/228439‑228574
                                                                     |                                                                  
aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgtg                                                                   >  1:1044959/1‑71 (MQ=255)
aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgt                                                                    >  1:1334701/1‑70 (MQ=255)
 ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                                  <  1:1845/71‑1 (MQ=255)
 ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                                  <  1:1843239/71‑1 (MQ=255)
   ctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                                >  1:3898814/1‑71 (MQ=255)
   ctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTAAATTTGTGttt                                                                >  1:3251103/1‑71 (MQ=255)
   ctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTAAATTTGTGttt                                                                >  1:3246444/1‑71 (MQ=255)
       gCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGc                                                             <  1:1266696/70‑1 (MQ=255)
        cTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCt                                                           <  1:1702383/71‑1 (MQ=255)
        cTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCt                                                           <  1:3022681/71‑1 (MQ=255)
         tttATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCtt                                                          >  1:1626583/1‑71 (MQ=255)
          ttattaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCttt                                                         >  1:2477113/1‑71 (MQ=255)
           tattaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCttt                                                         <  1:3693190/70‑1 (MQ=255)
             ttaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGGTGATCGTTAATTtgtg                                                                   >  1:3904983/1‑58 (MQ=255)
               aCCGGCAAACTCTGCGTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:2876531/68‑1 (MQ=255)
                ccGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGaaa                                                     >  1:373982/1‑69 (MQ=255)
                ccGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGaaa                                                     >  1:766122/1‑69 (MQ=255)
                ccGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAAcc                                                   <  1:3232582/71‑1 (MQ=255)
                ccGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAAc                                                    <  1:112809/70‑1 (MQ=255)
                 cGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGc                                                             <  1:1266698/60‑1 (MQ=255)
                  ggCAAACTCTGCTTCCCAACGACCATAACTAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                                 <  1:1667525/71‑1 (MQ=255)
                  ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                                  <  1:545640/54‑1 (MQ=255)
                  ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                                 <  1:2956999/71‑1 (MQ=255)
                  ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCa                                                  >  1:1931688/1‑70 (MQ=255)
                   gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCg                                                              >  1:3803766/1‑57 (MQ=255)
                   gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCg                                                              >  1:309208/1‑57 (MQ=255)
                     aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAAc                                              >  1:3156295/1‑71 (MQ=255)
                      aaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:2880943/61‑1 (MQ=255)
                      aaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:3210592/61‑1 (MQ=255)
                      aaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:1977697/61‑1 (MQ=255)
                      aaCTCTGCTTCCCAACGACCATAACCAGACAGTTGATCGTTAATTTGTGTTTCGcc                                                            <  1:2793746/56‑1 (MQ=255)
                       aCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                                             >  1:3924915/1‑70 (MQ=255)
                       aCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                                             >  1:1077088/1‑70 (MQ=255)
                         tctGCTTCCCAACGACCATAACCAGTCCGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:3609622/58‑1 (MQ=255)
                         tctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                       <  1:3020037/58‑1 (MQ=255)
                            gCTTCCCAACGACCATATCCAGTCAGTTGATCGTTAATTTGTGtt                                                                 <  1:2165696/45‑1 (MQ=255)
                            gCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGtt                                                                 <  1:780597/45‑1 (MQ=255)
                            gCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGtt                                                                 <  1:3150124/45‑1 (MQ=255)
                            gCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGtt                                                                 <  1:1561232/45‑1 (MQ=255)
                            gCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCa                                                  <  1:1836206/60‑1 (MQ=255)
                             cTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGatat                                        <  1:3692783/69‑1 (MQ=255)
                             cTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATaa                                      >  1:3150205/1‑71 (MQ=255)
                              ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                     >  1:3209535/1‑71 (MQ=255)
                              ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                     >  1:425616/1‑71 (MQ=255)
                                cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCt                                                           <  1:1874380/47‑1 (MQ=255)
                                cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGaaa                                                     >  1:2679860/1‑53 (MQ=255)
                                cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAActc                                   >  1:3380070/1‑71 (MQ=255)
                                cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                     <  1:551892/69‑1 (MQ=255)
                                 ccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAAcc                                                   <  1:1667624/54‑1 (MQ=255)
                                  cAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCttt                                                         >  1:301102/1‑47 (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:2812419/1‑70 (MQ=255)
                                      gACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTc                              >  1:3404868/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)
                                          aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                                 <  1:3191900/47‑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)
                                            aaCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGatat                                        <  1:2028902/54‑1 (MQ=255)
                                            aaCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGat                                          >  1:89456/1‑52 (MQ=255)
                                             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:3127936/1‑71 (MQ=255)
                                                      ttGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGt             >  1:292341/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:43192/1‑71 (MQ=255)
                                                             ttAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACt      <  1:1135909/71‑1 (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:1862248/54‑1 (MQ=255)
                                                                 tttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTAc                   <  1:1970715/54‑1 (MQ=255)
                                                                  ttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCata  >  1:787080/1‑70 (MQ=255)
                                                                  ttGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCata  >  1:1125748/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)
                                                                     |                                                                  
ACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATA  >  minE/228439‑228574

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: