Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 760,485 G→A 100% G594S (GGC→AGC)  sucA → subunit of E1(0) component of 2‑oxoglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913760,4850GA99.2% 39.5 / ‑7.1 18G594S (GGC→AGC) sucAsubunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/9);  minor base T (0/1);  total (8/10)
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.93e-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.

GTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTA  >  NC_000913/760350‑760621
                                                                                                                                       |                                                                                                                                        
agaatgaTGTAGCGTGTGTAGGAGCAGGCGAAACGCAGCAGCACGGTGCCGGAAGCAGTTGAAATGCAGGCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGc                                                                                                                                   <  1:269682/138‑1 (MQ=255)
 ttGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGCAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGACCCTGTTCGACTGGAGCGGTGCg                                                                                                                                  >  1:35609/1‑143 (MQ=255)
          aGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAGATGCAGTCTCGCGTTGCCACGCTTTATGGCGACCGCCAGGCGACGGCTGCCGGCGTGAAACTGTTCGACAGGAGCGGTGCGGGACAccac                                                                                                                         >  2:178789/1‑141 (MQ=255)
           gcgcCTGCGGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCg                                                                                                                        <  2:348977/143‑1 (MQ=255)
                tGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTAc                                                                                                                   <  2:151562/143‑1 (MQ=255)
                                      cAGCACGGTGGCGGAAGCAGTTGAAATGCAGGCTCGCGGTGCCAAGATTTATGGCGAGCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCa                                                                                             <  2:194084/143‑1 (MQ=255)
                                         cACGGTGCCGGAAGCAGCTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGCAACTGTTCGACTGGAGCGGTGCGGAAGACCTCGCTTACGCCACGCTGGCTGGTGAGGGCAGTc                                                                                          >  2:225390/1‑143 (MQ=255)
                                                 cGGAAGCAGGTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGc                                                                                  <  2:345175/143‑1 (MQ=255)
                                                  ggAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGAGCGGAAAACCTCGCTTACGCCACGCTGGTTGAGGAAGGCATTCCGGtgaccc                                                                                 >  2:39973/1‑138 (MQ=255)
                                                  ggAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCAACAACAGAGAAACTGCTCCACTGGAGCGGAGCGGAAAACCTCGCTTACGCCCCGCTGGCAGATGAAGGCACACCGGtcggca                                                                                 >  2:238220/1‑138 (MQ=255)
                                                             aaaTGCGGTCTCGCGTTGCCAAGATTTATGGCTATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGaa                                                                      <  2:163752/143‑1 (MQ=255)
                                                                    gTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACGCCg                                                               >  2:351949/1‑143 (MQ=255)
                                                                          cgTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTcgcg                                                         >  1:171885/1‑143 (MQ=255)
                                                                                cAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGATCGGTGCGGAAAACCTCACTTACGCCACGCTGATTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACct                                                   >  1:220404/1‑143 (MQ=255)
                                                                                                 gCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCGCGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCAcgcg                                  <  2:289870/143‑1 (MQ=255)
                                                                                                    gggCGGTGGCTGCCGGTGGGACGCGTTTCTGCGGGTGCTGTGCCGAAATCCGGGCTTGCGCCACGCTGGGTGATGGAGGCATTCCGGTTCGCCTGTCGGGTGAGGCCTCCAGGCGCGGTACCTTCTTCCACCGCCActttttc                               <  2:298386/142‑7 (MQ=255)
                                                                                                          tGGCTCCCGGGCAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCATGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCcac                         <  1:46918/143‑1 (MQ=255)
                                                                                                                                 gACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTa  <  2:171885/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTA  >  NC_000913/760350‑760621

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