Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,696 G→A 18.8% A49A (GCG→GCA yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,6960GA18.8% 123.1 / 16.4 79A49A (GCG→GCAyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (25/39);  new base A (7/8);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

AGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCT  >  minE/1781636‑1781761
                                                            |                                                                 
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:3330570/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:1233766/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:2527170/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:73585/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:2515442/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:1926195/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:1447130/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:3415152/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:790044/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:62715/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                           <  1:1915364/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                          <  1:1425186/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                          <  1:1568876/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                          <  1:141617/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                          <  1:720225/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                          <  1:1150781/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGTCACTAttgtt                                                    <  1:2981319/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:1908303/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:2577819/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:2844175/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:3250214/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:142868/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:1426249/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                                    <  1:795239/69‑1 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGc                                                             >  1:2427475/1‑56 (MQ=37)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCggg                                                 >  1:2574209/1‑68 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                                >  1:358299/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                                >  1:2081009/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                                >  1:1910621/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                                >  1:562729/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                                >  1:1077078/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                               >  1:2650515/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                               >  1:1487992/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                               >  1:206445/1‑69 (MQ=255)
                      ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                                     >  1:2387500/1‑69 (MQ=255)
                      ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                                     >  1:1378778/1‑69 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:964196/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:675052/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:635200/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:1895089/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:3033260/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:194329/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:2835426/68‑1 (MQ=255)
                           cGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                  <  1:3029205/67‑1 (MQ=255)
                            ggCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGAt                               >  1:2041821/1‑69 (MQ=255)
                                  aaCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGc                         >  1:1583131/1‑69 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          >  1:12834/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:2505300/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:3083514/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:2233189/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:187804/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          >  1:2963043/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:2883012/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:2846702/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:552721/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:1066511/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          <  1:419064/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          >  1:2692847/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                          >  1:505960/1‑67 (MQ=255)
                                       ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc                    <  1:2593498/69‑1 (MQ=255)
                                       ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc                    <  1:2381755/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                 <  1:2935050/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                 <  1:763862/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                 <  1:1239142/69‑1 (MQ=255)
                                            atatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                 <  1:2522312/67‑1 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAtt                >  1:2963841/1‑67 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTg               >  1:2494517/1‑68 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATCGc              >  1:539002/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCg                        >  1:2816003/1‑58 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc              >  1:2399167/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc              >  1:2087856/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc              >  1:2820023/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg             >  1:477606/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg             >  1:187557/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg             >  1:3341696/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg             >  1:2841937/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGGTGGGATGCGttt                                         >  1:2824633/1‑41 (MQ=255)
                                               ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa            >  1:3246407/1‑69 (MQ=255)
                                               ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa            >  1:1280030/1‑69 (MQ=255)
                                                           cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt  <  1:39859/67‑1 (MQ=255)
                                                           cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt  <  1:556798/67‑1 (MQ=255)
                                                           cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt  <  1:1527102/67‑1 (MQ=255)
                                                           cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt  <  1:568684/67‑1 (MQ=255)
                                                           cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt  <  1:2904288/67‑1 (MQ=255)
                                                            |                                                                 
AGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCT  >  minE/1781636‑1781761

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: