Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 349,262 C→T A194V (GCC→GTC)  prpB → 2‑methylisocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,2620CT97.8% 149.1 / ‑4.1 46A194V (GCC→GTC) prpB2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTG  >  NC_000913/349184‑349338
                                                                              |                                                                            
aTGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCa                                                                            >  2:67655/1‑81 (MQ=255)
aTGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCa                                                                            >  2:67632/1‑81 (MQ=255)
aTGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCa                                                                            >  2:67662/1‑81 (MQ=255)
   cggcgATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATg                                                                          >  2:67631/1‑80 (MQ=255)
   cggcgATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGt                                                                         <  2:67575/81‑1 (MQ=255)
    ggcgATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTa                                                                        >  2:67652/1‑81 (MQ=255)
      cgATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTAcc                                                                      <  1:67576/81‑1 (MQ=255)
          cGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGc                                                                    >  2:67648/1‑79 (MQ=255)
           gAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCa                                                                  <  2:67577/80‑1 (MQ=255)
            aGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGt                                                                >  2:67630/1‑81 (MQ=255)
                  cgcAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGcc                                                           >  2:67645/1‑80 (MQ=255)
                    cAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGAt                                                        <  1:67578/81‑1 (MQ=255)
                      ggCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGc                                                      <  2:67579/81‑1 (MQ=255)
                       gCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCg                                                     <  2:67580/81‑1 (MQ=255)
                        ccTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCgg                                                    >  1:67636/1‑81 (MQ=255)
                         cTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGt                                                   <  1:67581/81‑1 (MQ=255)
                         cTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGt                                                   >  1:67651/1‑81 (MQ=255)
                            tGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTg                                                  >  2:67650/1‑79 (MQ=255)
                            tGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCa                                                <  2:67582/81‑1 (MQ=255)
                                gAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGt                                             <  1:67583/80‑1 (MQ=255)
                                  aGCGGGTGCCGGGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGcc                                          <  2:67584/81‑1 (MQ=255)
                                     gggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGAt                                       <  2:67585/81‑1 (MQ=255)
                                      ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATc                                      >  2:67654/1‑81 (MQ=255)
                                      ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATc                                      <  1:67586/81‑1 (MQ=255)
                                      ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATc                                      >  2:67646/1‑81 (MQ=255)
                                          ccGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTc                                   <  2:67589/80‑1 (MQ=255)
                                          ccGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTc                                   <  1:67588/80‑1 (MQ=255)
                                          ccGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTc                                   <  1:67587/80‑1 (MQ=255)
                                               aTGTTGTTCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCaa                              >  1:67642/1‑80 (MQ=255)
                                                    gttCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCAc                        >  2:67665/1‑81 (MQ=255)
                                                     ttCCCGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCAc                        <  1:67590/80‑1 (MQ=255)
                                                       cccGGAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCg                      >  2:67644/1‑80 (MQ=255)
                                                          ggAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAAtt                  >  2:67657/1‑81 (MQ=255)
                                                           gAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAAttt                 >  2:67647/1‑81 (MQ=255)
                                                           gAGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAAttt                 >  1:67721/1‑81 (MQ=255)
                                                            aGGCGATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTg                >  2:67667/1‑81 (MQ=255)
                                                                gATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGc            >  2:67656/1‑81 (MQ=255)
                                                                gATTACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGc            >  2:67649/1‑81 (MQ=255)
                                                                   tACCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAc         >  1:67676/1‑81 (MQ=255)
                                                                    aCCGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAc         <  2:67591/80‑1 (MQ=255)
                                                                     ccGAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcg        <  2:67592/80‑1 (MQ=255)
                                                                       gAACTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccg     >  2:67675/1‑81 (MQ=255)
                                                                        aaCTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccgc    >  1:67643/1‑81 (MQ=255)
                                                                         aCTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCt   <  2:67594/81‑1 (MQ=255)
                                                                         aCTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCt   <  2:67593/81‑1 (MQ=255)
                                                                          cTCGTCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTg  >  1:67653/1‑81 (MQ=255)
                                                                              |                                                                            
ATGCGGCGATCGAGCGTGCGCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTG  >  NC_000913/349184‑349338

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