Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,861 C→A 14.1% R403L (CGG→CTG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,8610CA14.1% 344.3 / 28.4 120R403L (CGG→CTG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (52/51);  new base A (5/12);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 135 total aligned reads displayed.

TTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGA  >  Exported/1014792‑1014931
                                                                     |                                                                      
ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                                       >  1:3490139/1‑71 (MQ=255)
ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                                       >  1:2099960/1‑71 (MQ=255)
ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                                       >  1:1742759/1‑71 (MQ=255)
     cGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                  <  1:386899/71‑1 (MQ=255)
      ggACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                                 <  1:4819662/71‑1 (MQ=255)
         cccGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAg                                                              >  1:488525/1‑71 (MQ=255)
         cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                              >  1:4683493/1‑71 (MQ=255)
               tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                                <  1:5550914/63‑1 (MQ=255)
                 gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                                      >  1:591435/1‑71 (MQ=255)
                 gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                                   <  1:2701980/58‑1 (MQ=255)
                 gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                                   <  1:669181/58‑1 (MQ=255)
                 gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                              >  1:4561844/1‑63 (MQ=255)
                  taAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                      <  1:5614760/70‑1 (MQ=255)
                    aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                      >  1:2767455/1‑68 (MQ=255)
                    aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                   <  1:3039072/71‑1 (MQ=255)
                    aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                   >  1:3063084/1‑71 (MQ=255)
                     aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                         <  1:571955/64‑1 (MQ=255)
                     aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                         <  1:1891405/64‑1 (MQ=255)
                      tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                                  <  1:5657746/54‑1 (MQ=255)
                      tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                                         <  1:5341876/63‑1 (MQ=255)
                       cGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCaa                                                >  1:76916/1‑71 (MQ=255)
                        ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                              >  1:3956383/1‑56 (MQ=255)
                        ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                              >  1:4045622/1‑56 (MQ=255)
                        ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                                <  1:2865223/70‑1 (MQ=255)
                        ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                               <  1:5541178/71‑1 (MQ=255)
                          atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                                             <  1:401998/71‑1 (MQ=255)
                           tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                                          <  1:845235/57‑1 (MQ=255)
                           tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                                          <  1:793295/57‑1 (MQ=255)
                           tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATc                                            <  1:3127483/71‑1 (MQ=255)
                           ggatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                              >  1:5177737/2‑69 (MQ=255)
                            gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATc                                            >  1:2108383/1‑70 (MQ=255)
                             atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                                  <  1:3620500/47‑1 (MQ=255)
                             atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                                  <  1:1608258/47‑1 (MQ=255)
                             atCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAatat                                             <  1:1186710/68‑1 (MQ=255)
                              tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                           >  1:1096361/1‑53 (MQ=255)
                              tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                           >  1:1742229/1‑53 (MQ=255)
                              tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                               <  1:5089026/65‑1 (MQ=255)
                              tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                                         >  1:3013531/1‑71 (MQ=255)
                              tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                                         >  1:4376413/1‑71 (MQ=255)
                               cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                                >  1:1976161/1‑47 (MQ=255)
                               cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                                >  1:3973381/1‑47 (MQ=255)
                               cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                                >  1:1393894/1‑47 (MQ=255)
                                 ttGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgacc                                      >  1:1243314/1‑71 (MQ=255)
                                  tGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                         <  1:1014956/67‑1 (MQ=255)
                                   gATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccc                                     <  1:3442358/70‑2 (MQ=255)
                                    aTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccca                                    <  1:2628690/70‑3 (MQ=255)
                                       tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACcgacg                                  >  1:2855539/1‑69 (MQ=255)
                                       tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                                <  1:1044178/71‑1 (MQ=255)
                                       tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                <  1:4260998/71‑1 (MQ=255)
                                       tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                <  1:2416040/71‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                           <  1:2640018/42‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                           <  1:4970889/42‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                           <  1:3486280/42‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                                <  1:3857005/53‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcga                                        >  1:1522583/1‑61 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATaa                              >  1:4373141/1‑71 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATa                               <  1:2024512/70‑1 (MQ=255)
                                         cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATa                               >  1:4221716/1‑70 (MQ=255)
                                           gCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGAGTGAGCAat                                               >  1:5339659/1‑52 (MQ=255)
                                            cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCg                           >  1:4918468/1‑71 (MQ=255)
                                             ggCAGCGCTTCACCGTTATCCATCCGCTCTTTCAGt                                                             <  1:3237745/36‑1 (MQ=255)
                                             ggCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt                                                             <  1:3230830/36‑1 (MQ=255)
                                               cAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                       >  1:5531401/1‑56 (MQ=255)
                                                 gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                                <  1:2042683/61‑1 (MQ=255)
                                                 gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg                      >  1:359039/1‑71 (MQ=255)
                                                 gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg                      >  1:1479363/1‑71 (MQ=255)
                                                 gCCCTTCACCGTAATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                                <  1:5290922/61‑1 (MQ=255)
                                                  cccTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                                         <  1:1365618/35‑1 (MQ=255)
                                                    cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAGTCGTGc                        >  1:1366287/1‑66 (MQ=255)
                                                    cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc                        >  1:3192557/1‑66 (MQ=255)
                                                    cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt                    >  1:5052641/1‑70 (MQ=255)
                                                    cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt                    >  1:1772105/1‑70 (MQ=255)
                                                     ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgacc                                      <  1:16583/51‑1 (MQ=255)
                                                      tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa                   <  1:5148124/69‑1 (MQ=255)
                                                      tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg                 >  1:4027284/1‑71 (MQ=255)
                                                      tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCAACGATAATCGTGCCGTTAAg                 >  1:174298/1‑71 (MQ=255)
                                                       cACCGTTATCCATCCGCTCTTTCAGTTTGGCTTGAGCAATATCAcgaga                                      >  1:3351006/1‑47 (MQ=39)
                                                       cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                                          >  1:2372162/1‑45 (MQ=255)
                                                       cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                                          >  1:2737524/1‑45 (MQ=255)
                                                        aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                   >  1:5164418/1‑35 (MQ=255)
                                                        aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCg               <  1:2110683/71‑1 (MQ=255)
                                                         ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                               >  1:4881699/1‑38 (MQ=255)
                                                         ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                               >  1:883362/1‑38 (MQ=255)
                                                          cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                         <  1:3509224/43‑1 (MQ=255)
                                                            ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                       >  1:4098501/1‑43 (MQ=255)
                                                            ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                       >  1:3658776/1‑43 (MQ=255)
                                                            ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATaa           <  1:4427120/71‑1 (MQ=255)
                                                            ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATaa           <  1:2487399/71‑1 (MQ=255)
                                                             tATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAg          >  1:4278265/1‑71 (MQ=255)
                                                               tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg        <  1:2951692/68‑1 (MQ=255)
                                                               tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg        <  1:792600/68‑1 (MQ=255)
                                                               tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg        <  1:1676986/68‑1 (MQ=255)
                                                                 cacccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgt       <  1:5103627/67‑1 (MQ=255)
                                                                   tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                       >  1:4689288/1‑52 (MQ=255)
                                                                   tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                       >  1:1322693/1‑52 (MQ=255)
                                                                   tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg      >  1:3549813/1‑69 (MQ=255)
                                                                     cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAg   >  1:933886/1‑70 (MQ=255)
                                                                     cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa  <  1:3308024/71‑1 (MQ=255)
                                                                     cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa  <  1:2569817/71‑1 (MQ=255)
                                                                     cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa  <  1:1470266/71‑1 (MQ=255)
                                                                     |                                                                      
TTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGA  >  Exported/1014792‑1014931

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: