Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 760,485 G→A 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,4850GA94.1% 77.5 / NA 35G594S (GGC→AGC) sucAsubunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (19/13);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACAC  >  NC_000913/760345‑760624
                                                                                                                                            |                                                                                                                                           
aCAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGCTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGCGAGAAAATGTTCGACTGgcgc                                                                                                                                           >  2:280772/1‑143 (MQ=255)
     gTTGAGAGGGGGCGCCGGCAGGGGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGc                                                                                                                                      <  2:801591/143‑1 (MQ=255)
      ttGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGCGAAACTGTTCGACTGGAGCGGGGCg                                                                                                                                     >  1:850158/1‑143 (MQ=255)
         agaTGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGCTCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGaa                                                                                                                                  >  1:544313/1‑143 (MQ=255)
         agaTGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGaa                                                                                                                                  >  2:955658/1‑143 (MQ=255)
                gcgcCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAGGGAGTTGAAATGCAGTGTCGCGTTGCCAAGATTTATGGCGATGGCCAGGCGATGGCTGCCGGTGAGACACTGTTCGACTGGGGCGGTGCGGAGAGCCTCg                                                                                                                           >  2:213894/1‑143 (MQ=255)
                 cgcCTGCATGAGCTGGCGAAACGCACCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGACACTGTTCGACTGGAGCGGTGCGGAAAACCTCGc                                                                                                                          <  2:216405/143‑1 (MQ=255)
                   cccgCAGGACCTGTCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCtt                                                                                                                        <  1:290653/140‑1 (MQ=255)
                     tGCAGGCGCTGGCGAGACGCATCAGGACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTAc                                                                                                                      <  1:552594/143‑1 (MQ=255)
                      gCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATGTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGACACTGTGCGGCAGCCCTCGCTTagg                                                                                                                     >  2:744323/1‑141 (MQ=255)
                            gCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGc                                                                                                               >  1:292043/1‑143 (MQ=255)
                                gCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGCTGCGGACAACCTCGCTTACGTCAGGCTGGt                                                                                                           >  2:908743/1‑143 (MQ=255)
                                              cACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGCAAACCTCGCTTGCGCCACGCTGGTTGGTGAAGGCCTTc                                                                                             >  1:24279/1‑143 (MQ=255)
                                              cACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGACAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGCAGGCATTc                                                                                             >  1:165982/1‑143 (MQ=255)
                                                   tGCCGGAAGCAGTGGAGATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGtt                                                                                        <  2:275250/143‑1 (MQ=255)
                                                              gTTGAAATGCAGTCTCGCGTTGCCAAGAGTTATGGCGGTCGCCAGGCGATGGCTGCCGGTGAGAACCTGTTCGACTGGAGCGGTGCGGAGAACCTCGCTTACGCCACGCGGGTTGATGAAGGCATTCCGGTTCGCGCGTCgag                                                                             >  2:587535/1‑141 (MQ=255)
                                                                tggaaTGCAGTCTCGCGTTGCCAAGATTGATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTg                                                                           <  2:629117/140‑1 (MQ=255)
                                                                 ggaaTGCAGCCGCGCGTTGCCAAGATGTATGGCGATCGCCGGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGa                                                                          <  2:689495/141‑1 (MQ=255)
                                                                 gAAATGCAGTCTCGCGTTGCCAAGCTTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGCGCGGAAACCCTCAGGTACGCCACGAGGGTTGGTGAAGGCATTAAGGTTCGCCTGTCGGGGGa                                                                          >  2:785718/1‑143 (MQ=255)
                                                                     tGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATTGCTGCCGGGGAGAAACTGTTCGACTCGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCATGTCGTGAGAAGAc                                                                      >  2:71562/1‑143 (MQ=255)
                                                                             cgcgTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTCCCGGTGAGGCACTGTTCGACTGCAGCGGTGCGGAACACCTCGCTTACGCCACGCTGGTTGATGGAGGCAGTCCGGCTCCGCTGTCGCGTGCACCCTCCGGTcg                                                              >  2:814109/1‑143 (MQ=255)
                                                                                  tGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTa                                                         <  1:666346/143‑1 (MQ=255)
                                                                                     gagggtttAGGGCGCTCGCCAGGCGATGGCGGGCGCGGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCGACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACct                                                      <  2:379357/138‑1 (MQ=255)
                                                                                             aTGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAc                                              <  1:930619/143‑1 (MQ=255)
                                                                                             aTGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAc                                              <  2:415488/143‑1 (MQ=255)
                                                                                                        cAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTGACGCCACGCTGGCTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGt                                   >  1:573264/1‑143 (MQ=255)
                                                                                                         aGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTg                                  >  1:110498/1‑143 (MQ=255)
                                                                                                             gATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATcc                              >  1:825635/1‑143 (MQ=255)
                                                                                                               tGGCTGCCGGTGAGAAGCTGTTCGACTGGAGCGGTGCGGAAATCCTCGCTGACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCcac                            <  2:633365/143‑1 (MQ=255)
                                                                                                                   tGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTACGCCTGACGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCCCGCGGTGATCCACAAcc                        >  2:373370/1‑143 (MQ=255)
                                                                                                                         tGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTGCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTa                  >  1:345546/1‑143 (MQ=255)
                                                                                                                          gagaAACTGGTCGACTGGAGCGGTGCGGAAAGCCTCGCTTACGCCACGCTGGTTGATGAGGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTaa                 <  2:178696/143‑1 (MQ=255)
                                                                                                                               aCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCCGCCGGGTGAAGACTCCGGTCGCGGTACCTTCGTCCACCGGCACGCGCGGAGCCACAACCAGGCGAACgcct            >  2:24444/1‑140 (MQ=255)
                                                                                                                                      gACTGGAGCGTTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGGAGGCGTTCCGGTTCGCCTGTCGGCTGAAGACTCCGGACGCGGTCGCGTCTTCCACCGCCACGCGGCGATCCCCACCCAGCCGAACCGTTCCAcgcc     >  2:494212/1‑140 (MQ=255)
                                                                                                                                      gACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGCAGACTCCGGTCGCGGTGCCTTCTCCCACCGCCTAGCGGTGTTCCACGACCAGTCTAGCGGTTGCACTTa     >  1:197792/1‑143 (MQ=255)
                                                                                                                                         tGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTacac  <  2:528568/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
ACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACAC  >  NC_000913/760345‑760624

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: