Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,953 T→G 7.8% A372A (GCA→GCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,9530TG 7.8% 158.9 / 9.7 64A372A (GCA→GCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (34/25);  new base G (3/2);  total (37/27)
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.89e-01

CAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCG  >  minE/1014883‑1015018
                                                                      |                                                                 
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:3102014/1‑71 (MQ=255)
     tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:1233175/1‑71 (MQ=255)
     tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:307299/1‑71 (MQ=255)
     tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:3304851/1‑71 (MQ=255)
     tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                              >  1:144169/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:1596212/71‑1 (MQ=255)
             gacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCtt                                                      >  1:2447373/1‑71 (MQ=255)
                gaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCa                                                   >  1:2998540/1‑71 (MQ=255)
                  tAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCg                                                                <  1:2524493/56‑1 (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)
                             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:3914917/1‑65 (MQ=255)
                                   gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                      >  1:1494838/1‑65 (MQ=255)
                                   gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                      >  1:1300614/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:2739327/41‑1 (MQ=255)
                                          gtgtAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                                       <  1:1189401/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:3759899/1‑42 (MQ=255)
                                                          cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAg                                      >  1:1618018/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)
                                                                   cTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAGGCGCACCGCTTCGCCTTCTCCCGCTTTgcg  >  1:2612827/1‑69 (MQ=255)
                                                                    tgtgCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCtt      >  1:193981/1‑64 (MQ=255)
                                                                      |                                                                 
CAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCG  >  minE/1014883‑1015018

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: