Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,991 G→T 12.0% R360S (CGC→AGC) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,9910GT12.0% 156.5 / 6.4 50R360S (CGC→AGC) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (27/17);  new base T (4/2);  total (31/19)
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.78e-01

AAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTT  >  minE/1014921‑1015061
                                                                      |                                                                      
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)
             cGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCt                                                           <  1:2180474/71‑1 (MQ=255)
                gATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTtct                                                        <  1:2101188/71‑1 (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)
                      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)
                                                      cGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTcc                                   >  1:2893468/1‑54 (MQ=255)
                                                        gTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                                             <  1:658069/42‑1 (MQ=255)
                                                        gTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                <  1:1378592/71‑1 (MQ=255)
                                                         ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAtt               <  1:1909983/71‑1 (MQ=255)
                                                          tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                            >  1:3177009/1‑57 (MQ=255)
                                                          tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa              >  1:899932/1‑71 (MQ=255)
                                                          tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa              >  1:2330732/1‑71 (MQ=255)
                                                           aaGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa              <  1:141893/70‑1 (MQ=255)
                                                            aGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTg                                               <  1:2378082/36‑1 (MQ=255)
                                                            aGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCa                                           >  1:237372/1‑40 (MQ=255)
                                                              gTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata                           >  1:1888592/1‑54 (MQ=255)
                                                                cAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTc         >  1:2154091/1‑70 (MQ=255)
                                                                 aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGcc                     <  1:2509229/57‑1 (MQ=255)
                                                                 aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                <  1:1811078/62‑1 (MQ=255)
                                                                   cgcgCACCGCTTCGCCTTCTCCAGCTTTTAGCCGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGt     >  1:68083/1‑71 (MQ=255)
                                                                    gcgcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata                           <  1:1960902/48‑1 (MQ=255)
                                                                      gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt  <  1:2405387/71‑1 (MQ=255)
                                                                      gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt  <  1:2237799/71‑1 (MQ=255)
                                                                      gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt  >  1:165342/1‑71 (MQ=255)
                                                                      gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt  <  1:1504462/71‑1 (MQ=255)
                                                                      gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt  <  1:1440185/71‑1 (MQ=255)
                                                                      |                                                                      
AAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTT  >  minE/1014921‑1015061

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: