Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,949 A→G 7.7% L374L (TTG→CTG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,9490AG 7.7% 180.7 / 6.4 65L374L (TTG→CTG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base G (3/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCG  >  minE/1014881‑1015018
                                                                    |                                                                     
aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                                     >  1:1672560/1‑71 (MQ=255)
  cAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                                   <  1:3369983/71‑1 (MQ=255)
     tatCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGc                                                                 >  1:2048384/1‑70 (MQ=255)
     tatCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCg                                                                <  1:3427878/71‑1 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGa                                                               <  1:3551243/70‑1 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:1233175/1‑71 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:307299/1‑71 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:3102014/1‑71 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:144169/1‑71 (MQ=255)
       tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:3304851/1‑71 (MQ=255)
         aCGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGa                                                            <  1:105015/71‑1 (MQ=255)
           gaccgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                                          <  1:1713631/71‑1 (MQ=255)
            accgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtct                                                         <  1:4000449/71‑1 (MQ=255)
               gacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCtt                                                      >  1:2447373/1‑71 (MQ=255)
               gacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCtt                                                      <  1:1596212/71‑1 (MQ=255)
                  gaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCa                                                   >  1:2998540/1‑71 (MQ=255)
                    tAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCg                                                                <  1:2524493/56‑1 (MQ=255)
                     aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                                      >  1:2194138/1‑49 (MQ=255)
                     aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAAAGGGATACTGCGACAACTGTGCGAGGATCTCTTTc                                                    >  1:3635067/1‑67 (MQ=255)
                      atggTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCg                                                <  1:3178752/67‑1 (MQ=255)
                         gTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACg                                            >  1:3123220/1‑71 (MQ=255)
                           gCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATc                                                 >  1:3545611/1‑64 (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)
                               ttAAGCGATAAGCGTGTAGAAACGGGATACTGCGACaa                                                                       <  1:2348954/38‑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)
                                  aGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTc                                                    >  1:3654381/1‑54 (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)
                                            gtgtAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                                       <  1:1189401/41‑1 (MQ=255)
                                            gtgtAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                                       <  1:2739327/41‑1 (MQ=255)
                                            gtgtAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCAcc                         >  1:780386/1‑71 (MQ=255)
                                               tAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                                          <  1:1286507/35‑1 (MQ=255)
                                               tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTc                                                    <  1:2778625/41‑1 (MQ=255)
                                               tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCt                      >  1:3112831/1‑71 (MQ=255)
                                                 gAAACGGGATACTGCGACAACTGTGCGAGGAtctc                                                        >  1:3384699/1‑35 (MQ=255)
                                                 gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTaa                                       >  1:3276116/1‑52 (MQ=255)
                                                 gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                                        <  1:1870742/51‑1 (MQ=255)
                                                   aaCGGGATACTGCGACAACTGTGCGAGGATCTCttt                                                     <  1:2843280/36‑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:1618018/1‑42 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                      >  1:3759899/1‑42 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                                   <  1:548376/45‑1 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                                >  1:3275765/1‑48 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                                >  1:188160/1‑48 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCg                        >  1:1738502/1‑56 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg         >  1:1534982/1‑71 (MQ=255)
                                                            cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg         >  1:1089664/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)
                                                                    |                                                                     
AGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCG  >  minE/1014881‑1015018

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