Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,527 A→G 100% G354G (GGT→GGC)  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,685,5270AG98.1% 189.5 / ‑4.3 53G354G (GGT→GGC) fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/23);  total (30/23)
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.65e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGA  >  NC_000913/1685381‑1685661
                                                                                                                                                  |                                                                                                                                      
ttGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGcc                                                                                                                                      >  2:321824/1‑149 (MQ=255)
ttGAGTAATGGATTGATTCGCTCACCGTCCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACCTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTTTGGATCACCATTGGACGGAAGACGTTCAGCTCAAAGTTGcc                                                                                                                                      >  4:667433/1‑149 (MQ=255)
   aGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGa                                                                                                                                   >  2:584250/1‑149 (MQ=255)
      aaTTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAgc                                                                                                                                <  2:460032/149‑1 (MQ=255)
      aaTTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAgc                                                                                                                                >  1:460032/1‑149 (MQ=255)
          gattgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGccc                                                                                                                            >  3:648520/1‑149 (MQ=255)
                       aCGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCg                                                                                                               >  3:135330/1‑149 (MQ=255)
                          aTTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCa                                                                                                            <  4:648520/149‑1 (MQ=255)
                              ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtc                                                                                                        <  4:125272/149‑1 (MQ=255)
                                ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcgt                                                                                                      >  3:393027/1‑149 (MQ=255)
                                    aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                  <  2:243440/149‑1 (MQ=255)
                                        ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTCGATCGCCACGTCGTTCCCCATc                                                                                              >  1:52246/1‑149 (MQ=255)
                                           cTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAcc                                                                                           <  1:400336/149‑1 (MQ=255)
                                               gcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                                                       <  1:406747/149‑1 (MQ=255)
                                                 aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                     <  3:590079/149‑1 (MQ=255)
                                                 aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                     >  2:231807/1‑149 (MQ=255)
                                                     tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                 >  2:178708/1‑149 (MQ=255)
                                                        gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                              >  4:668611/1‑149 (MQ=255)
                                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                             >  1:56036/1‑149 (MQ=255)
                                                            aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACGTGACAGCAGAGCATGGt                                                                          >  2:369441/1‑149 (MQ=255)
                                                              aCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                        >  2:339718/1‑149 (MQ=255)
                                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                       <  2:159179/149‑1 (MQ=255)
                                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                       <  4:491717/149‑1 (MQ=255)
                                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                                                       >  3:204070/1‑149 (MQ=255)
                                                                  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                    >  1:422124/1‑149 (MQ=255)
                                                                  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGc                                                                    >  3:570637/1‑149 (MQ=255)
                                                                   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                   <  1:648069/149‑1 (MQ=255)
                                                                   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                   <  1:231807/149‑1 (MQ=255)
                                                                    cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                  >  1:106610/1‑149 (MQ=255)
                                                                     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                 >  3:146402/1‑149 (MQ=255)
                                                                        ccATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                              <  1:584250/149‑1 (MQ=255)
                                                                                  gCAAGCACACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                                                    >  4:165712/1‑149 (MQ=255)
                                                                                     aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTc                                                 <  1:137351/149‑1 (MQ=255)
                                                                                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCa                                                >  4:640605/1‑149 (MQ=255)
                                                                                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                   >  1:162354/1‑149 (MQ=255)
                                                                                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCAt                                   <  3:640605/149‑1 (MQ=255)
                                                                                                         ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                             <  3:165712/149‑1 (MQ=255)
                                                                                                         ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                             >  4:283220/1‑149 (MQ=255)
                                                                                                         ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCACGGCATGATTGa                             <  4:167213/149‑1 (MQ=255)
                                                                                                         ttGTGGATCACCATTGGACAGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                             >  2:400745/1‑149 (MQ=255)
                                                                                                            tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                       <  1:356541/120‑1 (MQ=255)
                                                                                                            tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                       >  2:356541/1‑120 (MQ=255)
                                                                                                                tCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGccc                      >  4:268071/1‑149 (MQ=255)
                                                                                                                      ttGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                >  3:225989/1‑149 (MQ=255)
                                                                                                                             gAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg         <  4:393027/149‑1 (MQ=255)
                                                                                                                             gAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg         <  4:135330/149‑1 (MQ=255)
                                                                                                                             gAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg         <  4:83495/149‑1 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                     <  2:158716/98‑1 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                     >  1:158716/1‑98 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                           <  4:364330/108‑1 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                           >  3:364330/1‑108 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga  <  4:93393/149‑1 (MQ=255)
                                                                                                                                    ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGACGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga  >  2:543694/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
TTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGA  >  NC_000913/1685381‑1685661

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