Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,229 C→G 10.9% A280A (GCG→GCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2290CG10.9% 149.8 / 8.2 55A280A (GCG→GCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base G (3/3);  total (28/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.96e-01

AATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTT  >  minE/1015166‑1015298
                                                               |                                                                     
aaTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGa                                                                <  1:2636938/71‑1 (MQ=255)
aaTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCGTCGa                                                                <  1:1510792/71‑1 (MQ=255)
 aTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                                                               >  1:2149742/1‑71 (MQ=255)
  tCACGCGGATGTCGTGAGCGAAGTATTTACCACCTAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                                                               >  1:1224255/1‑70 (MQ=255)
       cgGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc                                                                  >  1:255338/1‑62 (MQ=255)
           tGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCt                                                     <  1:2850582/71‑1 (MQ=255)
             tCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttt                                                   <  1:123056/71‑1 (MQ=255)
                tGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCg                                                                 <  1:2444911/54‑1 (MQ=255)
                 gAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc                                                >  1:1434402/1‑70 (MQ=255)
                   gCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAg                                             <  1:761127/71‑1 (MQ=255)
                      aaGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTgggc                                                                     >  1:1193010/1‑44 (MQ=12)
                      aaGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                          >  1:1230023/1‑71 (MQ=255)
                       aGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTcc                                         <  1:1825202/71‑1 (MQ=255)
                         tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCTCTTCGATCAGCAATTCTTTTTcc                                               <  1:3576004/63‑1 (MQ=255)
                         tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac                                       <  1:3951433/71‑1 (MQ=255)
                          aTTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcaca                                      >  1:1578601/1‑71 (MQ=255)
                          aTTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcaca                                      >  1:2323110/1‑71 (MQ=255)
                           tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtt                                                    <  1:3916735/56‑1 (MQ=255)
                             tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCa                                              >  1:3248751/1‑60 (MQ=255)
                             tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                   >  1:1425660/1‑71 (MQ=255)
                                 gccaAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAgc                                            <  1:38180/57‑2 (MQ=255)
                                  ccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttttt                                                 >  1:3846029/1‑52 (MQ=255)
                                  ccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttttt                                                 >  1:933603/1‑52 (MQ=255)
                                    aaaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACg                            <  1:1363336/71‑1 (MQ=255)
                                    aaaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACg                            >  1:3231890/1‑71 (MQ=255)
                                    aaaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCt                                                     >  1:924361/1‑46 (MQ=255)
                                     aaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGc                           >  1:2134061/1‑71 (MQ=255)
                                      aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGcc                          >  1:772568/1‑71 (MQ=255)
                                       cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                          >  1:1734605/1‑54 (MQ=255)
                                       cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                                     <  1:572389/59‑1 (MQ=255)
                                       cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                                     <  1:3433633/59‑1 (MQ=255)
                                       cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCt                         >  1:1216288/1‑71 (MQ=255)
                                        tGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCt                         >  1:31996/1‑70 (MQ=255)
                                         gcgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttt                                                   >  1:2006461/1‑43 (MQ=255)
                                           gcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAcc                     <  1:151129/71‑1 (MQ=255)
                                           gcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAcc                     <  1:164936/71‑1 (MQ=255)
                                            cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                    <  1:2768556/71‑1 (MQ=255)
                                            cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                    <  1:3950693/71‑1 (MQ=255)
                                             aCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                                     <  1:3234393/53‑1 (MQ=255)
                                               ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                 >  1:1908612/1‑71 (MQ=255)
                                                cAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                    <  1:2190500/67‑1 (MQ=255)
                                                 aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTg                        <  1:1935804/62‑1 (MQ=255)
                                                  gACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                        >  1:2516226/1‑45 (MQ=255)
                                                  gACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                        >  1:2463939/1‑45 (MQ=255)
                                                  gACAAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                        >  1:1728732/1‑45 (MQ=255)
                                                   aCCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                       >  1:1696140/1‑61 (MQ=255)
                                                       aGATTTTGCGCTTCGATCAGCAATTCTTTTTACAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTccc         >  1:3387117/1‑71 (MQ=255)
                                                        gATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                       >  1:1219349/1‑56 (MQ=255)
                                                          ttttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGaa                              >  1:2285161/1‑47 (MQ=255)
                                                          ttttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTc      >  1:1230565/1‑71 (MQ=255)
                                                           tttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTg                        <  1:3274110/52‑1 (MQ=255)
                                                           tttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTg                        <  1:920588/52‑1 (MQ=255)
                                                            ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc     <  1:691656/70‑1 (MQ=255)
                                                             tGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGa                               <  1:2296884/43‑1 (MQ=255)
                                                             tGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg    <  1:3796645/70‑1 (MQ=255)
                                                             tGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg    <  1:1863658/70‑1 (MQ=255)
                                                              gcgcTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                   <  1:617746/38‑1 (MQ=255)
                                                               cgcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGtt  <  1:3467956/70‑1 (MQ=255)
                                                               |                                                                     
AATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTT  >  minE/1015166‑1015298

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: