Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,974 A→G 13.1% R365R (CGT→CGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,9740AG13.1% 157.0 / 15.4 61R365R (CGT→CGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (33/20);  new base G (5/3);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

GTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTG  >  minE/1014906‑1015044
                                                                    |                                                                      
gTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACg                                                                      >  1:3123220/1‑71 (MQ=255)
   ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGtt                                                                   <  1:1281701/71‑1 (MQ=255)
   ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGt                                                                    <  1:83058/70‑1 (MQ=255)
   ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGt                                                                    <  1:2504190/70‑1 (MQ=255)
    cGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                                                                  <  1:610351/71‑1 (MQ=255)
      gtgaGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAg                                                                <  1:2741368/68‑1 (MQ=255)
       tAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAgg                                                               >  1:3111127/1‑71 (MQ=255)
        aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGt                                                              >  1:603516/1‑71 (MQ=255)
        aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGAAGGTTAAgg                                                               >  1:2035429/1‑70 (MQ=255)
          gCGATAAGCGTGAAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCa                                                            <  1:3038885/71‑1 (MQ=255)
            gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                                                >  1:1494838/1‑65 (MQ=255)
            gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                                                >  1:1300614/1‑65 (MQ=255)
            gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                                                >  1:3914917/1‑65 (MQ=255)
             aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGtt                                                                   <  1:2684268/61‑1 (MQ=255)
              tAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcgc                                                        <  1:161114/71‑1 (MQ=255)
               aaGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcgcg                                                       >  1:219754/1‑71 (MQ=255)
                aGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACgcgc                                                      <  1:2803762/71‑1 (MQ=255)
                   gtgtAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCAcc                                                   >  1:780386/1‑71 (MQ=255)
                      tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCt                                                >  1:3112831/1‑71 (MQ=255)
                        gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTaa                                                                 >  1:3276116/1‑52 (MQ=255)
                        gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                                                                  <  1:1870742/51‑1 (MQ=255)
                            cGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCt                                          <  1:2180474/71‑1 (MQ=255)
                             gggATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCaa                                                           >  1:3678980/1‑53 (MQ=255)
                               gATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTtct                                       <  1:2101188/71‑1 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                                                >  1:3759899/1‑42 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                                                >  1:1618018/1‑42 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                                                             <  1:548376/45‑1 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                                                          >  1:188160/1‑48 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                                                          >  1:3275765/1‑48 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCg                                                  >  1:1738502/1‑56 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg                                   >  1:1089664/1‑71 (MQ=255)
                                   cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg                                   >  1:1534982/1‑71 (MQ=255)
                                    tGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGc                                  >  1:1858761/1‑71 (MQ=255)
                                     gCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAAc                                                          >  1:1511350/1‑46 (MQ=37)
                                     gCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGcc                                 >  1:3034778/1‑70 (MQ=255)
                                      cGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGc                                  <  1:1757436/69‑1 (MQ=255)
                                       gACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCttt                               >  1:704309/1‑71 (MQ=255)
                                          aaCTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                            >  1:591585/1‑71 (MQ=255)
                                            cTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAGGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                            >  1:2612827/1‑69 (MQ=255)
                                             tgtgCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCtt                                >  1:193981/1‑64 (MQ=255)
                                                  gAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcgc                           <  1:3066474/64‑1 (MQ=255)
                                                     gATCTCTTTCATCGGACGGTTAAGGTGAACGCGCACCGCTTCGCCt                                          <  1:2470173/46‑1 (MQ=255)
                                                      aTCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTccc                                    <  1:1930734/51‑1 (MQ=255)
                                                       tctcTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCt                                          >  1:3132604/1‑44 (MQ=255)
                                                        ctctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAg                         <  1:2703865/60‑1 (MQ=255)
                                                        ctctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTc                   >  1:818005/1‑66 (MQ=255)
                                                         tctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAg                         <  1:1604314/59‑1 (MQ=255)
                                                          ctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                       >  1:1701736/1‑58 (MQ=255)
                                                          ctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                       >  1:4069418/1‑58 (MQ=255)
                                                          ctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTcc                                     >  1:546026/1‑46 (MQ=255)
                                                          ctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGat             >  1:2833588/1‑70 (MQ=255)
                                                           tttCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGctc                      >  1:3367253/1‑60 (MQ=255)
                                                           tttCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCgg                <  1:2316345/66‑1 (MQ=255)
                                                             tCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                       >  1:1688759/1‑55 (MQ=255)
                                                             tCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata          <  1:2174847/70‑1 (MQ=255)
                                                             tCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGata            <  1:4284/68‑1 (MQ=255)
                                                              cATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATAtt        >  1:1361354/1‑71 (MQ=255)
                                                               aTCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATAttt       <  1:4089249/71‑1 (MQ=255)
                                                                 cGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGc     >  1:1349687/1‑71 (MQ=255)
                                                                  ggACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGcc    >  1:3936515/1‑71 (MQ=255)
                                                                   gACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCt   >  1:328254/1‑71 (MQ=255)
                                                                    aCGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTg  >  1:612937/1‑71 (MQ=255)
                                                                    |                                                                      
GTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTG  >  minE/1014906‑1015044

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: