Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,680,559 A→G 43.0% V1472A (GTC→GCC)  yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,680,5590AG43.0% 22.7 / 60.5 79V1472A (GTC→GCC) yfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/26);  new base G (14/20);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

AACCGCTGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAAGAG  >  minE/1680497‑1680627
                                                              |                                                                    
aaCCGCTGGCATCCACACGCAGCCATAACGGCTGAT‑‑CCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:984209/69‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2672436/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1854918/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1807321/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2975892/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1771258/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:3276389/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2884621/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1638464/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1637584/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2678465/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1559184/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2493122/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:940804/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2315802/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1073083/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:120479/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:115245/68‑1 (MQ=255)
   cGCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:2438868/68‑1 (MQ=255)
    gCTGGCATCCACACGCAGCCATAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGa                                                              <  1:1391552/67‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:3020711/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:2546047/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:230676/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:2222294/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:1965859/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:1084908/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:1290848/69‑1 (MQ=255)
      tGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGc                                                          <  1:58859/69‑1 (MQ=255)
            ccATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCa                                                    <  1:2031353/69‑1 (MQ=255)
            ccATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCa                                                    <  1:490225/69‑1 (MQ=255)
            ccATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCa                                                    <  1:3418691/69‑1 (MQ=255)
             cATACGCAACCATATCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAg                                                   >  1:2936162/1‑69 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCACGGTGAc                                                                     <  1:441434/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:1572722/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:2682271/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:2147845/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:2924345/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:13408/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:315440/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGAc                                                                     <  1:24394/51‑1 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:2878141/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:795794/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:804600/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:2337065/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:231525/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:2233902/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:1534627/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:1856276/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTg                                                               >  1:1549046/1‑57 (MQ=255)
             cATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGCTGACAAGTTg                                                               >  1:1562701/1‑57 (MQ=255)
                       cATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGttt                                         <  1:2812084/69‑1 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:3431623/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:375963/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:168742/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:3203811/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:1695328/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:3062129/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:1771819/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:535648/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:782745/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:2858118/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:2251698/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:226992/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:2634178/1‑69 (MQ=255)
                        aTAACGGCTGATCGCCACTGTTGGTCACCTGCAAGGTGGCAAGTTGATCGCTATTCAGATTGCTGgttt                                        >  1:1449309/1‑69 (MQ=255)
                           aCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTgcg                                     >  1:1426422/1‑69 (MQ=255)
                           aCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTgcg                                     >  1:991023/1‑69 (MQ=255)
                           aCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTgcg                                     >  1:3375013/1‑69 (MQ=255)
                           aCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTgcg                                     >  1:3083090/1‑69 (MQ=255)
                           aCGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGATATTCAGATTGCTGTTTTgcg                                     >  1:2210371/1‑69 (MQ=255)
                                                         aGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgag        <  1:1949601/68‑1 (MQ=255)
                                                         aGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgag        <  1:2043046/68‑1 (MQ=255)
                                                         aGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgag        <  1:650446/68‑1 (MQ=255)
                                                         aGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgag        <  1:2304288/68‑1 (MQ=255)
                                                         aGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgag        <  1:1160334/68‑1 (MQ=255)
                                                           gTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgaga       >  1:2462477/1‑67 (MQ=255)
                                                           gTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTgaga       >  1:2393599/1‑67 (MQ=255)
                                                            tGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAag    >  1:2152357/1‑69 (MQ=255)
                                                            tGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTC‑‑CTGCTCAGCTGAGAAagag  <  1:1460597/69‑1 (MQ=255)
                                                              aCAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAagag  >  1:1658336/1‑69 (MQ=255)
                                                              aCAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAagag  >  1:638516/1‑69 (MQ=255)
                                                              aCAAGTTGATCGCTATTCAGATTGCTGTTTTGAGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAagag  >  1:474337/1‑69 (MQ=255)
                                                              |                                                                    
AACCGCTGGCATCCATACGCAACCATAACGGCTGATCGCCACTGTTGCTCACCTGCAAGGTGACAAGTTGATCGCTATTCAGATTGCTGTTTTGCGCTTTCTCGCCTGTCAGCTGCTCAGCTGAGAAAGAG  >  minE/1680497‑1680627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: