Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,709 C→A 9.3% L456L (CTG→CTT lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,7090CA 9.3% 206.1 / 4.9 64L456L (CTG→CTTlldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base A (4/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTG  >  minE/2562640‑2562776
                                                                     |                                                                   
gCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCa                                                                    <  1:3304159/71‑1 (MQ=255)
  cAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAg                                                                   >  1:1955246/1‑70 (MQ=255)
    aTTTGTTGTGCTGCGGTGTCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                                <  1:3506186/71‑1 (MQ=255)
    aTTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                                >  1:961878/1‑71 (MQ=255)
     tttGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                                >  1:2686338/1‑70 (MQ=255)
     tttGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                                >  1:684665/1‑70 (MQ=255)
      ttgttgTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAAt                                                               <  1:514334/70‑1 (MQ=255)
        gttgtGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAgg                                                                  <  1:3162486/65‑1 (MQ=255)
        gttgtGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAgg                                                                  <  1:1949096/65‑1 (MQ=255)
              ctgcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAg                                                                   >  1:1939526/1‑58 (MQ=255)
                 cGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcca                                                    <  1:2200930/70‑1 (MQ=255)
                       cTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCAAGCCAGCCGAgg                                             <  1:1715006/71‑1 (MQ=255)
                        ttGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                                       <  1:1569462/60‑1 (MQ=255)
                        ttGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAgg                                             >  1:1927720/1‑70 (MQ=255)
                          gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATAc                                                             >  1:1903400/1‑52 (MQ=255)
                          gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                                       >  1:619563/1‑58 (MQ=255)
                          gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                                       >  1:1131862/1‑58 (MQ=255)
                             gcgcGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGc                                       <  1:1087162/71‑1 (MQ=255)
                              cgcgGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagccag                                                   >  1:3417283/1‑58 (MQ=255)
                                 ggcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATAccc                                                           >  1:3138321/1‑47 (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)
                                       cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcc                                                     >  1:866392/1‑47 (MQ=255)
                                       cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcc                                                     >  1:988312/1‑47 (MQ=255)
                                       cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                  >  1:3869225/1‑66 (MQ=255)
                                       cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                  >  1:3873898/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)
                                                aGATGAGGTATCCGACCCGGTCAGGAATACCCccagcca                                                    >  1:606549/1‑39 (MQ=255)
                                                aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCa                     <  1:850830/70‑1 (MQ=255)
                                                aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                    >  1:3828736/1‑71 (MQ=255)
                                                 gATGAGGTATCCGACCCGGTCAGGAATACCCccagccag                                                   <  1:3538920/39‑1 (MQ=255)
                                                 gATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                    <  1:995191/70‑1 (MQ=255)
                                                  aggAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACg                                         <  1:1772796/46‑1 (MQ=19)
                                                  aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa                  >  1:1929560/1‑71 (MQ=255)
                                                  aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa                  >  1:3147750/1‑71 (MQ=255)
                                                   tGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGa                                            >  1:2189709/1‑44 (MQ=255)
                                                    gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                                    <  1:3308563/51‑1 (MQ=255)
                                                    gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                                    <  1:1882487/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:2318651/71‑1 (MQ=255)
                                                               cccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa     <  1:2371450/71‑1 (MQ=255)
                                                                  ggTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt   >  1:974396/1‑70 (MQ=255)
                                                                  ggTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTg  >  1:1907396/1‑71 (MQ=255)
                                                                     |                                                                   
GCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTG  >  minE/2562640‑2562776

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: