Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,659 C→A 12.9% V470V (GTG→GTT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6590CA12.9% 169.4 / 7.7 54V470V (GTG→GTTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base A (4/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACA  >  minE/1014590‑1014728
                                                                     |                                                                     
cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                                       >  1:266624/1‑70 (MQ=255)
cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                                                      >  1:3186191/1‑71 (MQ=255)
      tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                                   >  1:3255236/1‑68 (MQ=255)
        cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                              <  1:398530/71‑1 (MQ=255)
        cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                              <  1:495580/71‑1 (MQ=255)
             gCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggc                                                             <  1:1049621/67‑1 (MQ=255)
             gCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                         <  1:3287757/71‑1 (MQ=255)
               cGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                         <  1:2535967/69‑1 (MQ=255)
                ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                         <  1:2467512/68‑1 (MQ=255)
                ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                      <  1:3375748/71‑1 (MQ=255)
                ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                      <  1:3619571/71‑1 (MQ=255)
                     gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                                                                 >  1:3591259/2‑55 (MQ=255)
                     gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCt                                                  >  1:320690/2‑70 (MQ=255)
                     accGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCtt                                                 <  1:3607258/71‑1 (MQ=255)
                      ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                                >  1:3542282/1‑71 (MQ=255)
                      ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                                >  1:1037621/1‑71 (MQ=255)
                      ccGATACTGCCAATGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                                   >  1:705339/1‑52 (MQ=255)
                      ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                                   >  1:1518011/1‑52 (MQ=255)
                       cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                               <  1:2295753/71‑1 (MQ=255)
                         aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCt                                                  <  1:1736244/66‑1 (MQ=255)
                         aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                             <  1:1788725/71‑1 (MQ=255)
                          tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                            <  1:1873513/71‑1 (MQ=255)
                            cTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAg                                          <  1:448507/71‑1 (MQ=255)
                                cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                            <  1:1702727/65‑1 (MQ=255)
                                 aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTg                                                     <  1:1357149/55‑1 (MQ=255)
                                  aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGAt                                    >  1:1908248/1‑71 (MQ=255)
                                   ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATc                                   <  1:1931665/71‑1 (MQ=255)
                                     tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctg                                                                  <  1:2254137/38‑1 (MQ=255)
                                     gaGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCgg                                                        <  1:1323498/47‑1 (MQ=255)
                                       gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                                <  1:2703189/70‑2 (MQ=255)
                                               cgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                       <  1:3049470/71‑2 (MQ=255)
                                                gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                       <  1:1691351/70‑2 (MQ=255)
                                                   ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg                   >  1:3358151/1‑71 (MQ=255)
                                                    tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                      >  1:2549758/1‑67 (MQ=255)
                                                    tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                      >  1:1178129/1‑67 (MQ=255)
                                                     ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                 >  1:2443948/1‑71 (MQ=255)
                                                     ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                 >  1:117427/1‑71 (MQ=255)
                                                     ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                 >  1:1725149/1‑71 (MQ=255)
                                                      tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc                    >  1:66361/1‑67 (MQ=255)
                                                      tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                >  1:3504993/1‑71 (MQ=255)
                                                       ttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                <  1:63543/70‑1 (MQ=255)
                                                         aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg               >  1:3463334/1‑69 (MQ=255)
                                                         aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg               >  1:1779442/1‑69 (MQ=255)
                                                             gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTg         <  1:1115746/71‑1 (MQ=255)
                                                               aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt       >  1:256349/1‑71 (MQ=255)
                                                               aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt       >  1:688963/1‑71 (MQ=255)
                                                               aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCACCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa        >  1:1528744/1‑70 (MQ=255)
                                                               aTCCGTCACCTGCTGGCTGCGGTTGCCTTTAGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt       >  1:1630231/1‑71 (MQ=255)
                                                                tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACtt                             <  1:842147/48‑1 (MQ=255)
                                                                tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc      >  1:3298645/1‑71 (MQ=255)
                                                                tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc      >  1:2316786/1‑71 (MQ=255)
                                                                  cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCg     <  1:3010248/70‑1 (MQ=255)
                                                                  cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa    <  1:96148/71‑1 (MQ=255)
                                                                   gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTc                            >  1:3073524/1‑46 (MQ=255)
                                                                    tCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt          <  1:1729560/63‑1 (MQ=255)
                                                                     cACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa  >  1:1543197/1‑70 (MQ=255)
                                                                     |                                                                     
CGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACA  >  minE/1014590‑1014728

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: