Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,013 T→C 9.3% K352K (AAA→AAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0130TC 9.3% 204.8 / 11.8 75K352K (AAA→AAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (4/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGG  >  minE/1014945‑1015082
                                                                    |                                                                     
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)
                 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)
                   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)
                                             ttcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                        >  1:1674801/3‑71 (MQ=255)
                                              gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                       <  1:2405387/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)
                                              gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                       <  1:2237799/71‑1 (MQ=255)
                                                aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTg                                        <  1:865198/52‑1 (MQ=255)
                                                aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                     >  1:99078/1‑71 (MQ=255)
                                                aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                     <  1:869609/71‑1 (MQ=255)
                                                 ccGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGata                                                  <  1:3080140/41‑1 (MQ=255)
                                                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCg                    >  1:423275/1‑70 (MQ=255)
                                                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                   >  1:2705091/1‑71 (MQ=255)
                                                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                   >  1:3842882/1‑71 (MQ=255)
                                                    cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                                                 <  1:4109325/39‑1 (MQ=255)
                                                    cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                                                 <  1:3500533/39‑1 (MQ=255)
                                                    cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                 <  1:3260101/71‑1 (MQ=255)
                                                    cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                 <  1:3447399/71‑1 (MQ=255)
                                                       cGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa               <  1:326395/70‑1 (MQ=255)
                                                        gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                   <  1:1672012/49‑1 (MQ=255)
                                                          cTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCt                                         >  1:2978340/1‑41 (MQ=255)
                                                           ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                     <  1:4140745/44‑1 (MQ=255)
                                                           ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                 <  1:3598271/64‑1 (MQ=255)
                                                           ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcc          >  1:274429/1‑71 (MQ=255)
                                                            tctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATccc         >  1:417230/1‑71 (MQ=255)
                                                             ctcCCGCTTTGCGGAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                      >  1:3477920/1‑57 (MQ=255)
                                                             ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                      >  1:3289995/1‑57 (MQ=255)
                                                             ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                      >  1:727207/1‑57 (MQ=255)
                                                             ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc        >  1:2102682/1‑71 (MQ=255)
                                                                ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                  >  1:1201972/1‑58 (MQ=255)
                                                                ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                 <  1:1035314/59‑1 (MQ=255)
                                                                ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                 <  1:1437471/59‑1 (MQ=255)
                                                                ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTg      <  1:1254966/70‑1 (MQ=255)
                                                                 cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc    >  1:1971624/1‑71 (MQ=255)
                                                                 cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc    >  1:1914875/1‑71 (MQ=255)
                                                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTTGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg  <  1:4101072/71‑1 (MQ=255)
                                                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                       <  1:3314699/50‑1 (MQ=255)
                                                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg  <  1:1286753/71‑1 (MQ=255)
                                                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg  <  1:2770773/71‑1 (MQ=255)
                                                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg  <  1:3294437/71‑1 (MQ=255)
                                                                    tttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg  <  1:2962426/70‑1 (MQ=255)
                                                                    |                                                                     
GACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGG  >  minE/1014945‑1015082

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: