Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,686,586 C→A 8.0% M1M (ATG→ATT) † fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,686,5860CA 8.0% 151.3 / 3.8 50M1M (ATG→ATT) †fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (28/18);  new base A (2/2);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTAT  >  NC_000913/1686443‑1686734
                                                                                                                                               |                                                                                                                                                    
gcCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGAc                                                                                                                                                 <  4:160694/149‑1 (MQ=255)
 ccAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGAc                                                                                                                                                 <  1:145054/148‑1 (MQ=255)
 ccAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGAc                                                                                                                                                 >  2:145054/1‑148 (MQ=255)
  cAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCt                                                                                                                                               >  7:413855/1‑149 (MQ=255)
    gcgcATGACTCAGTGAGGTGCGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGc                                                                                                                                             <  4:328524/149‑1 (MQ=255)
                 tGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGAt                                                                                                                                 >  1:54898/1‑148 (MQ=255)
                      tGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAAttt                                                                                                                           >  8:262010/1‑149 (MQ=255)
                       gggCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAAtttt                                                                                                                          >  6:337416/1‑149 (MQ=255)
                         gCATTTGCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAAtttttt                                                                                                                        <  4:371378/149‑1 (MQ=255)
                           aTTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTAATGAc                                                                                                                                                 >  6:262288/1‑122 (MQ=255)
                           aTTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTAATGAc                                                                                                                                                 >  6:257656/1‑122 (MQ=255)
                           aTTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTAATGAc                                                                                                                                                 <  5:262288/122‑1 (MQ=255)
                           aTTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTAATGAc                                                                                                                                                 <  5:257656/122‑1 (MQ=255)
                                ctcCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTttctttct                                                                                                                  <  8:413855/148‑1 (MQ=255)
                                   cGTCGAAATGCGGAACTGCTCCAGCGAGCGTGGAGCGTGGGCGCCCCACAGCTGAGCTGCCGGGACATCAATCGCCCCCCTCGACTCTTTTTCGCTGCGTACGGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGt                                                                                                                <  3:117193/147‑1 (MQ=255)
                                   cGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGttt                                                                                                              <  6:37704/149‑1 (MQ=255)
                                              ggAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGtt                                                                                                    >  2:35639/1‑148 (MQ=255)
                                              ggAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGtt                                                                                                    >  8:286866/1‑148 (MQ=255)
                                                 aaTGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGc                                                                                                <  5:337416/149‑1 (MQ=255)
                                                     cTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACt                                                                                            >  3:208515/1‑149 (MQ=255)
                                                        cAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTACGCAACtt                                                                                           <  2:318632/147‑1 (MQ=255)
                                                        cAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACtttt                                                                                         <  5:407017/149‑1 (MQ=255)
                                                          gcgagcgTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAg                                                                                       >  7:221714/1‑149 (MQ=255)
                                                                    aGTTGGTGCGCCCCACAGGTTAGCTCCCGGCAAATCACTCCCCCCCTTCCAATCTTTTGCGCTGCGTACTTTTTTCATGCCCTGCCCCGCACCTTATTAATTTTTTCTTTCTGTTTGGGTTTCGTAACCCACCTTTTCCCTCACtaat                                                                              <  2:54898/148‑1 (MQ=255)
                                                                             gcCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCGCCCCTGATTAAGTTTTTCTTTGTGTTTTGCTTTCGTTAAGCACGTTTTAGCTCAGTTATTATTTACCAt                                                                    >  8:139477/1‑149 (MQ=255)
                                                                              ccccACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTGACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCAtt                                                                   <  5:288542/149‑1 (MQ=255)
                                                                              ccccACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCAtt                                                                   <  8:324273/149‑1 (MQ=255)
                                                                               cccACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCAttt                                                                  >  5:399585/1‑149 (MQ=255)
                                                                                caacaGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCAttt                                                                  >  2:176889/3‑148 (MQ=255)
                                                                                 cacaGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCAt                                                                    >  7:151461/1‑145 (MQ=255)
                                                                                  acaGCTTATCTGCCGGGCCATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGa                                                                >  7:125132/1‑148 (MQ=255)
                                                                                      cTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAAc                                                            >  1:101374/1‑148 (MQ=255)
                                                                                      cTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCATCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAAATTTTAGCTCACTTATTATTTACCATTTGATAAc                                                            >  2:225915/1‑148 (MQ=255)
                                                                                       ttATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACa                                                           >  7:219010/1‑148 (MQ=255)
                                                                                         aTCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAAt                                                        >  8:14069/1‑149 (MQ=255)
                                                                                              ccGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGttt                                                    <  7:262010/148‑1 (MQ=255)
                                                                                              ccGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTg                                                   <  1:35639/149‑1 (MQ=255)
                                                                                                gggACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTgtt                                                 >  6:126617/1‑147 (MQ=255)
                                                                                                   aCATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCtt                                              >  8:177143/1‑149 (MQ=255)
                                                                                                     aTCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTc                                            >  2:52098/1‑149 (MQ=255)
                                                                                                                      aaTCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGACATGTAAAAAAAc                           >  5:80562/1‑149 (MQ=255)
                                                                                                                          attttCGCTGCGAACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGcc                       >  4:323585/2‑149 (MQ=255)
                                                                                                                             ttCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGcccc                     >  3:243181/1‑148 (MQ=255)
                                                                                                                                gCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATTTAAAAAAACCGCCCCGaa                  <  8:125132/148‑1 (MQ=255)
                                                                                                                                      tACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGggcggc           >  6:51129/1‑149 (MQ=255)
                                                                                                                                           tATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGctct        <  8:221714/147‑1 (MQ=255)
                                                                                                                                           tATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGt      >  7:265317/1‑149 (MQ=255)
                                                                                                                                            aTTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGtt     >  3:155096/1‑149 (MQ=255)
                                                                                                                                              tCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTc                                            <  6:190167/108‑1 (MQ=255)
                                                                                                                                              tCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTc                                            >  5:190167/1‑108 (MQ=255)
                                                                                                                                               cATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTAt  >  1:276649/1‑149 (MQ=255)
                                                                                                                                               cATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTAt  >  8:415511/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTATCTGCCGGGACATCAATCGCCCCCATCGAATCTTTTTCGCTGCGTACTGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTAT  >  NC_000913/1686443‑1686734

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