Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,733 G→A 8.5% F448F (TTC→TTT lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,7330GA 8.5% 243.7 / 8.7 82F448F (TTC→TTTlldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base G (38/37);  new base A (5/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCG  >  minE/2562663‑2562802
                                                                      |                                                                     
cTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCAAGCCAGCCGAgg                                                                       <  1:1715006/71‑1 (MQ=255)
 ttGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAgg                                                                       >  1:1927720/1‑70 (MQ=255)
      gcgcGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGc                                                                 <  1:1087162/71‑1 (MQ=255)
           gcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                                            <  1:1758913/71‑1 (MQ=255)
              aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACggt                                                                 >  1:1096630/1‑62 (MQ=255)
              aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACgg                                                                  >  1:2521122/1‑62 (MQ=255)
              aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGagaag                                                           <  1:83398/69‑1 (MQ=255)
                cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                                            >  1:3873898/1‑66 (MQ=255)
                cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                                            >  1:3869225/1‑66 (MQ=255)
                  gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGaa                                                                     >  1:3295031/1‑55 (MQ=255)
                   ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                                                    <  1:1327871/71‑1 (MQ=255)
                   ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                                                    <  1:124496/71‑1 (MQ=255)
                     cGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAATGTGAAt                                                  <  1:505031/71‑1 (MQ=255)
                     cGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaa                                                   >  1:3290278/1‑70 (MQ=255)
                       ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                                                    <  1:1370431/67‑1 (MQ=255)
                       ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAAt                                                  <  1:2929700/69‑1 (MQ=255)
                       ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGc                                                <  1:2770260/71‑1 (MQ=255)
                        tAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaacat                                                >  1:1457708/1‑67 (MQ=255)
                         aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCa                                               <  1:850830/70‑1 (MQ=255)
                         aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                                              >  1:3828736/1‑71 (MQ=255)
                          gATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                                              <  1:995191/70‑1 (MQ=255)
                           aggAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACg                                                                   <  1:1772796/46‑1 (MQ=19)
                           aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa                                            >  1:3147750/1‑71 (MQ=255)
                           aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa                                            >  1:1929560/1‑71 (MQ=255)
                            tGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGa                                                                      >  1:2189709/1‑44 (MQ=255)
                             gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                                                              <  1:1882487/51‑1 (MQ=255)
                             gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                                                              <  1:3308563/51‑1 (MQ=255)
                             gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGAcc                                          >  1:618523/1‑71 (MQ=255)
                              aGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCg                                         <  1:2615893/71‑1 (MQ=255)
                               ggTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAgg                                                       >  1:1576350/1‑56 (MQ=255)
                               ggTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCgg                                        <  1:1222223/71‑1 (MQ=255)
                                 tATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtg                                      <  1:721815/71‑1 (MQ=255)
                                  aTCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATg                                             <  1:2799648/63‑1 (MQ=255)
                                  aTCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtgt                                     >  1:3640938/1‑71 (MQ=255)
                                   tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaaga                                                          >  1:1271012/1‑49 (MQ=255)
                                   tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATg                                             <  1:1101573/62‑1 (MQ=255)
                                    ccGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                                              >  1:3584682/1‑60 (MQ=255)
                                      gACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc                                 <  1:1735141/71‑1 (MQ=255)
                                       aCCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc                                 >  1:154673/1‑70 (MQ=255)
                                        cccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa                               <  1:2371450/71‑1 (MQ=255)
                                        cccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa                               <  1:2318651/71‑1 (MQ=255)
                                           ggTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt                             >  1:974396/1‑70 (MQ=255)
                                           ggTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTg                            >  1:1907396/1‑71 (MQ=255)
                                                 gAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTg                                                     <  1:1736152/40‑1 (MQ=255)
                                                 gAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTg                                                     <  1:3160904/40‑1 (MQ=255)
                                                 gAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCTGTGTGCGCCAGTGCCAgcg                      >  1:2761252/1‑71 (MQ=255)
                                                 gAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcg                      >  1:1748500/1‑71 (MQ=255)
                                                  aaTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc                     >  1:3154040/1‑71 (MQ=255)
                                                   aTACTCCCAGCCAGCCGAGAAACGGTGAgaagaac                                                        >  1:332702/1‑34 (MQ=25)
                                                   aTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc                     >  1:1022546/1‑70 (MQ=255)
                                                   aTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc                     >  1:1247580/1‑70 (MQ=255)
                                                    tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc                     >  1:2790105/1‑69 (MQ=255)
                                                    tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCa                   <  1:925556/71‑1 (MQ=255)
                                                      cccccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt                             <  1:2733568/59‑1 (MQ=255)
                                                      cccccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAgt                 <  1:2901893/71‑1 (MQ=255)
                                                       ccccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAAt                                                  <  1:1279357/37‑1 (MQ=255)
                                                        cccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAgtgt               >  1:661893/1‑71 (MQ=255)
                                                        cccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAgtgt               <  1:679790/71‑1 (MQ=255)
                                                         ccagccagCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtg                                      >  1:1106044/1‑47 (MQ=255)
                                                         ccagccagCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:1472290/1‑71 (MQ=255)
                                                          cagccagcCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:3657369/1‑70 (MQ=255)
                                                          cagccagcCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:2741785/1‑70 (MQ=255)
                                                          cagccagcCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:2600238/1‑70 (MQ=255)
                                                          cagccagcCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:2404804/1‑70 (MQ=255)
                                                          cagccagcCGAGGAAAGGCGAGAAGAAGGTGAATGCATGACCg                                         >  1:2072616/1‑43 (MQ=255)
                                                           agccagccGAGGAACGGCGAGAAGAAGGTGAATGCa                                               >  1:1916934/1‑36 (MQ=255)
                                                           agccagccGAGGAACGGCGAGAAGAAGGTGAATGCa                                               >  1:2186394/1‑36 (MQ=255)
                                                           agccagccGAGGAACGGCGAGAAGAAGGTGAATGCa                                               >  1:3944/1‑36 (MQ=255)
                                                             ccagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtgt                                     <  1:3798058/44‑1 (MQ=255)
                                                             ccagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAg                              >  1:2814234/1‑51 (MQ=255)
                                                              cagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtgatga         >  1:869274/1‑71 (MQ=255)
                                                              cagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtgatga         <  1:3822908/71‑1 (MQ=255)
                                                              cagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtgatga         <  1:530700/71‑1 (MQ=255)
                                                              cagccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtgat           <  1:3447356/69‑1 (MQ=255)
                                                               agccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtga            <  1:718904/67‑1 (MQ=255)
                                                               agccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtga            <  1:1631394/67‑1 (MQ=255)
                                                                gccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACa       <  1:385335/71‑1 (MQ=255)
                                                                 ccGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAgtg                <  1:615884/61‑1 (MQ=255)
                                                                  cGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTtga            >  1:2449187/1‑64 (MQ=255)
                                                                  cGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGt     >  1:2545349/1‑71 (MQ=255)
                                                                  cGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGt     >  1:2901704/1‑71 (MQ=255)
                                                                   gAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGtt              >  1:1059032/1‑61 (MQ=255)
                                                                   gAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTc    <  1:2389987/71‑1 (MQ=255)
                                                                   gAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTc    <  1:362665/71‑1 (MQ=255)
                                                                     ggAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAg                  >  1:419162/1‑55 (MQ=255)
                                                                     ggAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCg  <  1:2302443/71‑1 (MQ=255)
                                                                     ggAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCCCCAGTGTTGATGAGAGTCCg  <  1:505821/71‑1 (MQ=255)
                                                                      gAACGGCGAGAAGAAGGTGAATGCATGACCGgtgt                                     >  1:1414417/1‑35 (MQ=255)
                                                                      |                                                                     
CTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCG  >  minE/2562663‑2562802

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: