Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 228,520 G→A 7.7% F61F (TTC→TTT phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE228,5200GA 7.7% 278.4 / 8.7 90F61F (TTC→TTTphoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base G (47/36);  new base A (4/3);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTA  >  minE/228454‑228589
                                                                  |                                                                     
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)
   ggCAAACTCTGCTTCCCAACGACCATAACTAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                                                <  1:1667525/71‑1 (MQ=255)
   ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                                                <  1:2956999/71‑1 (MQ=255)
   ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCa                                                                 >  1:1931688/1‑70 (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)
        aCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                                                            >  1:1077088/1‑70 (MQ=255)
        aCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                                                            >  1:3924915/1‑70 (MQ=255)
          tctGCTTCCCAACGACCATAACCAGTCCGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                                      <  1:3609622/58‑1 (MQ=255)
          tctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                                                      <  1:3020037/58‑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:425616/1‑71 (MQ=255)
               ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                                                    >  1:3209535/1‑71 (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)
                    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:2628137/71‑1 (MQ=255)
                              aCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCtt                                     <  1:963391/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: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)
                                                        tttCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATg            <  1:1654923/70‑1 (MQ=255)
                                                          tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc                                       >  1:753096/1‑41 (MQ=255)
                                                          tCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc                                       >  1:1778927/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:1628621/1‑49 (MQ=255)
                                                                ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTAt                         >  1:3632668/1‑49 (MQ=255)
                                                                ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt    >  1:2704674/1‑70 (MQ=255)
                                                                ttGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt    >  1:3665960/1‑70 (MQ=255)
                                                                 tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt   >  1:1626692/1‑70 (MQ=255)
                                                                 tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCttt   >  1:3320600/1‑70 (MQ=255)
                                                                 tGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCtt    <  1:2217319/69‑1 (MQ=255)
                                                                  gAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTa  <  1:2628581/70‑1 (MQ=255)
                                                                  |                                                                     
ACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTA  >  minE/228454‑228589

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: