Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,089,913 C→T 7.7% Q257* (CAG→TAG)  fumC‑I → class 2 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,089,9130CT 7.7% 148.6 / 3.1 52Q257* (CAG→TAG) fumC‑Iclass 2 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGAAATCGAGCTG  >  NC_002947/1089769‑1090062
                                                                                                                                                |                                                                                                                                                     
cacaTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCAGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAAc                                                                                                                                                  <  1:1628052/150‑1 (MQ=255)
   aTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTg                                                                                                                                               >  1:568317/1‑150 (MQ=255)
      gAAGGTGTGATGCAGCAAATGCAGACGCTGGCGCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTAGCTCCCCAGCTCAGTGGCCTGACCAAGCTCGAGTTCACGCCCGGCCAGAACCTGTTc                                                                                                                                            >  2:918890/1‑150 (MQ=255)
       aaGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGACGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCg                                                                                                                                           >  1:306785/1‑150 (MQ=255)
          gtgtGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCAc                                                                                                                                        <  1:1578521/150‑1 (MQ=255)
              gctgcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGAt                                                                                                                                    <  1:124394/150‑1 (MQ=255)
                  ccGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGc                                                                                                                                >  2:19898/1‑150 (MQ=255)
                        ttGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCGCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAg                                                                                                                          <  2:306785/150‑1 (MQ=255)
                                gcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGc                                                                                                                  >  2:530644/1‑150 (MQ=255)
                                   ggCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGt                                                                                                               <  2:943310/150‑1 (MQ=255)
                                      ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCACGCGGAGTTCGCTGGCGGTTTCTCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGCTCACGCCCGGCCCGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGc                                                                                                            <  2:371553/150‑1 (MQ=255)
                                        aGGGCGGAACCGCCGTGGGTACCTGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCtt                                                                                                          <  2:627421/150‑1 (MQ=255)
                                                 ccgccgTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGcc                                                                                                 >  2:1191369/1‑150 (MQ=255)
                                                           tACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGaa                                                                                                                                                   >  1:1567477/1‑90 (MQ=255)
                                                           tACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGaa                                                                                                                                                   <  2:1567477/90‑1 (MQ=255)
                                                           tACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGac                                                                                       <  1:19898/150‑1 (MQ=255)
                                                           tACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGac                                                                                       >  1:491484/1‑150 (MQ=255)
                                                              cGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCTAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGaccac                                                                                    >  2:1392971/1‑150 (MQ=255)
                                                              cGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCTAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGaccac                                                                                    >  2:1392879/1‑150 (MQ=255)
                                                              cGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGaccac                                                                                    >  2:921830/1‑150 (MQ=255)
                                                              cGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGa                                                                                        >  1:1077890/1‑146 (MQ=255)
                                                              cGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGa                                                                                        <  2:1077890/146‑1 (MQ=255)
                                                                gCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCg                                                                                  >  2:374766/1‑150 (MQ=255)
                                                                gCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCg                                                                                  >  2:422640/1‑150 (MQ=255)
                                                                          ccATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGtt                                                                        >  2:691922/1‑150 (MQ=255)
                                                                            aTCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCGCCGCTGTGGCGTtga                                                                      <  2:568317/150‑1 (MQ=255)
                                                                                      cGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCTAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGc                                                            <  1:1392971/150‑1 (MQ=255)
                                                                                      cGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCTAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGc                                                            <  1:1392879/150‑1 (MQ=255)
                                                                                      cGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGc                                                            <  2:491484/150‑1 (MQ=255)
                                                                                       gCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGcc                                                           >  2:1591192/1‑150 (MQ=255)
                                                                                          gTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGt                                                                                                                                              <  1:1461941/64‑1 (MQ=255)
                                                                                          gTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGt                                                                                                                                              >  2:1461941/1‑64 (MQ=255)
                                                                                              gTTTCGCGCGCCAGCTCAGTGGCCTGACGAAGGTCGAGTTCACGCCCGGCCAGTACCTGTTCGCACTGATCGGCTGGCAGGACACTGCCGTGGCTTTGTCTGGCCAGGTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATc                                                    <  2:749009/150‑1 (MQ=255)
                                                                                               tttCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCt                                                   >  2:903628/1‑150 (MQ=255)
                                                                                                     tcgcCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGCTCCATGTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTAGCCAGCGCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTg                                             >  1:1118275/1‑150 (MQ=255)
                                                                                                     tcgcCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTg                                             >  1:257905/1‑150 (MQ=255)
                                                                                                     tcgcCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTg                                             >  2:1218669/1‑150 (MQ=255)
                                                                                                       gcCAGCTCAGTGGCCTTACCAAGGTCGAGTTCACGCCCGGCCAGACCATGATGTCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGACCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGa                                           <  1:918890/150‑1 (MQ=255)
                                                                                                        ccAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGAt                                          <  1:422640/150‑1 (MQ=255)
                                                                                                                   gCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCt                                              >  1:1637371/1‑135 (MQ=255)
                                                                                                                   gCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCt                                              <  2:1637371/135‑1 (MQ=255)
                                                                                                                      tGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGc                            >  1:184844/1‑150 (MQ=255)
                                                                                                                          cAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGc                        >  2:294644/1‑150 (MQ=255)
                                                                                                                          cAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGc                        >  2:260604/1‑150 (MQ=255)
                                                                                                                                  aGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTgg                <  2:184844/150‑1 (MQ=255)
                                                                                                                                    ttCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTgg                >  1:1068220/1‑148 (MQ=255)
                                                                                                                                    ttCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTgg                <  2:1068220/148‑1 (MQ=255)
                                                                                                                                     tCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGc                                <  2:327312/131‑1 (MQ=255)
                                                                                                                                     tCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGc                                >  1:327312/1‑131 (MQ=255)
                                                                                                                                       aCGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGaa           >  2:729192/1‑150 (MQ=255)
                                                                                                                                        cGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGaaa          <  2:823852/150‑1 (MQ=255)
                                                                                                                                              gCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCACCGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGAAATCGAGc    <  1:530644/150‑1 (MQ=255)
                                                                                                                                                cAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGAAATCGAGCTg  >  1:1496244/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGGCCAGCTCAAGACCACCGCTGTGGCGTTGATGAAAATCGCCAACGATCTGCGCTGGATGAACTCCGGCCCGCTGGCCGGTCTGGGCGAAATCGAGCTG  >  NC_002947/1089769‑1090062

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