Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,036,339 A→G 41.1% K157K (AAA→AAG)  ygiC → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,036,3390AG41.1% 32.3 / 64.8 90K157K (AAA→AAG) ygiChypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/32);  new base G (19/18);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTCA  >  minE/2036279‑2036397
                                                            |                                                          
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:3152481/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:676055/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:2017441/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:2063577/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:1720780/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:1565984/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:2320050/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:1229322/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                    >  1:254998/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcg                                                     >  1:103215/1‑68 (MQ=255)
     ttAACTCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                               >  1:2343121/1‑69 (MQ=255)
     ttAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                               >  1:1367525/1‑69 (MQ=255)
     ttAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                               >  1:491124/1‑69 (MQ=255)
       aaCGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTc                                             >  1:2610355/1‑69 (MQ=255)
       aaCGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTc                                             >  1:2931198/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                    >  1:3048031/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                    >  1:1521321/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                    >  1:2919771/1‑69 (MQ=255)
                    tGCCGTAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:2644300/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:1466657/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:150295/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:2171793/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:3048000/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                                <  1:3296028/69‑1 (MQ=255)
                         gAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAg                           <  1:3198696/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:2252682/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:3342313/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:2158342/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:863594/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:324277/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:1863169/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:1324150/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:3416464/69‑1 (MQ=255)
                           gggCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:1960240/68‑1 (MQ=255)
                             gCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATg                       <  1:1192694/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCTAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:2498135/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATGGCTTCGTTGAGCTGCGTGAACAGTTTGGCTTcc                 <  1:598027/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1186608/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:321276/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1242769/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:125600/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:3253249/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:325340/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:2420013/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:472924/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1529857/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:2236149/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1182799/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1605443/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1166083/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:338565/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:3190337/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1978607/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                          <  1:1300169/60‑1 (MQ=255)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:881993/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:336815/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:47179/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:3263440/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:3411040/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:332788/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:3372823/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:1227271/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:267378/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:2598489/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:127538/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:2430376/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:2175833/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:102259/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:2010844/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:1925768/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:1410256/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:1116988/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                 >  1:1652926/1‑46 (MQ=39)
                                          tAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTtgc           >  1:3112033/1‑68 (MQ=255)
                                            aCAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTtgctg         <  1:561235/68‑1 (MQ=255)
                                             cAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTtgctgc        <  1:680807/68‑1 (MQ=255)
                                             cAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTtgctgc        <  1:2465681/68‑1 (MQ=255)
                                             cAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCATCCAGTtgctgc        <  1:1880922/68‑1 (MQ=255)
                                              aGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:3347876/58‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1904285/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:1917411/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:377939/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                 <  1:503295/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc   <  1:614470/69‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc   <  1:100415/69‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc   <  1:1555188/69‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc   <  1:2112381/69‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc   <  1:2881489/69‑1 (MQ=255)
                                                  tGCAAGAAAAACTGATCGATCGCTTCGTTGATCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTCa  >  1:1410442/1‑69 (MQ=255)
                                                  tGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTCa  >  1:3368756/1‑69 (MQ=255)
                                                            |                                                          
TCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTCA  >  minE/2036279‑2036397

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