Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,036,328 C→T 46.8% L154L (CTG→TTG)  ygiC → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,036,3280CT46.8% 8.8 / 66.2 79L154L (CTG→TTG) ygiChypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (19/18);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTC  >  minE/2036279‑2036396
                                                 |                                                                    
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:2320050/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:254998/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:1720780/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:2017441/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:1565984/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:2063577/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:1229322/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:676055/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcga                                                   >  1:3152481/1‑69 (MQ=255)
tCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTgatcg                                                    >  1:103215/1‑68 (MQ=255)
     ttAACTCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                              >  1:2343121/1‑69 (MQ=255)
     ttAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                              >  1:491124/1‑69 (MQ=255)
     ttAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCt                                              >  1:1367525/1‑69 (MQ=255)
       aaCGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTc                                            >  1:2931198/1‑69 (MQ=255)
       aaCGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTc                                            >  1:2610355/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                   >  1:3048031/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                   >  1:2919771/1‑69 (MQ=255)
                aaCTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTg                                   >  1:1521321/1‑69 (MQ=255)
                    tGCCGTAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                               <  1:2644300/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                               <  1:2171793/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                               <  1:150295/69‑1 (MQ=255)
                    tGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTg                               <  1:1466657/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:2158342/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:3342313/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:2252682/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:1863169/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:3416464/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:324277/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:1324150/69‑1 (MQ=255)
                          aGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:863594/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:3253249/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:325340/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:2420013/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:321276/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:3190337/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1978607/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1166083/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1182799/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1186608/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:338565/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1242769/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1605443/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:2236149/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1529857/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:472924/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:125600/69‑1 (MQ=255)
                                   aCCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGt                         <  1:1300169/60‑1 (MQ=255)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:47179/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:881993/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:3411040/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:3372823/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:336815/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:3263440/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:332788/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:1116988/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:267378/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:2598489/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:1227271/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:127538/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:1925768/1‑46 (MQ=39)
                                          tAACAGTTTGCAAGAAAAGCTGATCGATCGCTTCGTTGAGCTGCGt                                >  1:1410256/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:377939/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:1917411/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTcc                <  1:503295/55‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc  <  1:100415/69‑1 (MQ=255)
                                                 cTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCAtctc  <  1:614470/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)
                                                 |                                                                    
TCAGCTTAACGCCGGTAACTTGCCGGAGGGCAGCGACCAGTTTAACAGTCTGCAAGAAAAACTGATCGATCGCTTCGTTGAGCTGCGTGAACAGTATGGCTTCCAGTTGCTGCATCTC  >  minE/2036279‑2036396

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