Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 454,695 C→G 7.5% E130D (GAG→GAC)  ybgH ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE454,6950CG 7.5% 279.3 / 8.1 93E130D (GAG→GAC) ybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (44/42);  new base G (5/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCC  >  minE/454625‑454764
                                                                      |                                                                     
cGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                                                       >  1:1172243/1‑71 (MQ=255)
cGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                                                       >  1:1569829/1‑71 (MQ=255)
   taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                                    >  1:2161032/1‑71 (MQ=255)
        gACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                                    <  1:184754/66‑1 (MQ=255)
           cccACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                                                                  <  1:859548/65‑1 (MQ=255)
           cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTc                                                            >  1:3847756/1‑71 (MQ=255)
             cACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGc                                                          >  1:608391/1‑71 (MQ=255)
               cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                                                        <  1:1312282/71‑1 (MQ=255)
                gTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGa                                                       <  1:2009645/71‑1 (MQ=255)
                 ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                      <  1:276893/71‑1 (MQ=255)
                   aCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                      >  1:373191/1‑69 (MQ=255)
                   aCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                      >  1:1799032/1‑69 (MQ=255)
                    cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                      >  1:2578201/1‑68 (MQ=255)
                    cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                                    <  1:2559078/70‑1 (MQ=255)
                     ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                  >  1:2965736/1‑71 (MQ=255)
                     ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                  >  1:2391635/1‑71 (MQ=255)
                      cgccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                  >  1:3608211/1‑70 (MQ=255)
                       gccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGcc                                                         <  1:2241698/62‑1 (MQ=255)
                       gccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                  >  1:3541892/1‑69 (MQ=255)
                       cccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCTGTTGGCTCATACAGCTCGcc                                                         <  1:3619297/61‑1 (MQ=255)
                        cagcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAg                                               >  1:2326625/3‑71 (MQ=255)
                         cgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGc                                              >  1:2477066/1‑71 (MQ=255)
                          gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             <  1:829819/71‑1 (MQ=255)
                          gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             >  1:405055/1‑71 (MQ=255)
                          gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             <  1:3281296/71‑1 (MQ=255)
                          gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             >  1:1069030/1‑71 (MQ=255)
                          gcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCt                                             >  1:3766968/1‑71 (MQ=255)
                           cATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCtt                                            <  1:3802070/71‑1 (MQ=255)
                            aTACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtga                                           <  1:1096568/71‑3 (MQ=255)
                              aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                                              <  1:2158513/50‑1 (MQ=255)
                              aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACg                                         >  1:3177783/1‑71 (MQ=255)
                              aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTAACg                                         >  1:165886/1‑71 (MQ=255)
                               cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        >  1:3877817/1‑71 (MQ=255)
                               cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        <  1:501747/71‑1 (MQ=255)
                               cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        <  1:1233663/71‑1 (MQ=255)
                                aTCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATAc                                                                 <  1:1081706/45‑1 (MQ=255)
                                  cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        >  1:2540867/1‑68 (MQ=255)
                                  cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        >  1:1202482/1‑68 (MQ=255)
                                   aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                        >  1:2064997/1‑67 (MQ=255)
                                    gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                                    <  1:3110265/38‑1 (MQ=255)
                                    gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             <  1:3103711/61‑1 (MQ=255)
                                     cGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                                                        >  1:2622821/1‑49 (MQ=255)
                                     cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                                  >  1:3321801/1‑71 (MQ=255)
                                     cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                                  >  1:402654/1‑71 (MQ=255)
                                      gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                      <  1:3643367/50‑1 (MQ=255)
                                      gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                                  >  1:1724687/1‑70 (MQ=255)
                                       agaAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                                    >  1:78781/1‑35 (MQ=255)
                                       agaAACCGCCATCACGACGCGGATCGGTTGGCTCATAc                                                                 >  1:2698670/1‑38 (MQ=255)
                                       agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                                              <  1:3723749/41‑1 (MQ=255)
                                       agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                                   >  1:1792505/1‑68 (MQ=255)
                                       aaaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttt                                <  1:407424/69‑2 (MQ=255)
                                        gaAACCGCCATCACGACGCGGATCGGTTGGCTCATTCAGCTCGCCGAGCAGACAGCTTACGTTAGAt                                   >  1:2123033/1‑67 (MQ=255)
                                        gaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaa                               <  1:559351/71‑1 (MQ=255)
                                           aCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGtt                                       >  1:2765125/1‑60 (MQ=255)
                                           aCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCg                                                           >  1:1895018/1‑40 (MQ=255)
                                            ccGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                           <  1:1102261/71‑1 (MQ=255)
                                            ccGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                           <  1:683214/71‑1 (MQ=255)
                                             cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             <  1:2614407/52‑1 (MQ=255)
                                             cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAgg                          >  1:2730951/1‑71 (MQ=255)
                                             cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                           <  1:3693585/70‑1 (MQ=255)
                                             cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCAGAGCAGACAGCTTACGTTAGATTTAAACAgg                          >  1:2724593/1‑71 (MQ=255)
                                              gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                                   >  1:3615566/1‑61 (MQ=255)
                                              gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                                   >  1:2097698/1‑61 (MQ=255)
                                               ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                             <  1:484344/50‑1 (MQ=255)
                                               ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                           >  1:783409/1‑52 (MQ=255)
                                               ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                           >  1:50847/1‑52 (MQ=255)
                                               ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                           >  1:968076/1‑52 (MQ=255)
                                                cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                                    <  1:3866321/42‑1 (MQ=255)
                                                cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa                       <  1:1801350/71‑1 (MQ=255)
                                                 aTCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGc                                                     >  1:3345456/1‑40 (MQ=255)
                                                 aTCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCAt                      <  1:813998/71‑1 (MQ=255)
                                                  tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa                       >  1:246003/1‑69 (MQ=255)
                                                  tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa                     <  1:3780608/71‑1 (MQ=255)
                                                  tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa                     <  1:2751043/71‑1 (MQ=255)
                                                   cACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACa                            >  1:968035/1‑63 (MQ=255)
                                                    acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                              >  1:1682499/1‑60 (MQ=255)
                                                    acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGc                   >  1:3451288/1‑71 (MQ=255)
                                                      gacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg                 <  1:5780/71‑1 (MQ=255)
                                                       acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                               >  1:2425160/1‑40 (MQ=255)
                                                       acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACATCTTACGTTAGATTTa                                <  1:2108123/55‑1 (MQ=255)
                                                        cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                               >  1:3792344/1‑39 (MQ=255)
                                                        cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                               >  1:870356/1‑39 (MQ=255)
                                                        cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg                 >  1:1228929/1‑69 (MQ=255)
                                                          cgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGa             <  1:2779507/71‑1 (MQ=255)
                                                             aTCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGa          <  1:2758569/71‑1 (MQ=255)
                                                              tCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGa          >  1:709944/1‑70 (MQ=255)
                                                              tCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGAt         <  1:476305/71‑1 (MQ=255)
                                                                 gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                                       <  1:2172007/38‑1 (MQ=255)
                                                                 gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                                 <  1:645129/44‑1 (MQ=255)
                                                                 gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                                 <  1:2049359/44‑1 (MQ=255)
                                                                     gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                              >  1:2483421/1‑43 (MQ=255)
                                                                     gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc  <  1:816773/71‑1 (MQ=255)
                                                                     gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc  <  1:3879563/71‑1 (MQ=255)
                                                                     gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc  <  1:386479/71‑1 (MQ=255)
                                                                     gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc   <  1:1096232/70‑1 (MQ=255)
                                                                     ccTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc   <  1:3247900/69‑1 (MQ=255)
                                                                      cTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc   >  1:281619/1‑69 (MQ=255)
                                                                      |                                                                     
CGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCC  >  minE/454625‑454764

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: ‑