Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,717,846 C→A 15.8% V249V (GTG→GTT paaG ← 2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,717,8460CA15.8% 90.6 / 10.0 38V249V (GTG→GTTpaaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base A (3/3);  total (18/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 = 4.77e-01

TGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGT  >  NC_002947/3717707‑3717991
                                                                                                                                           |                                                                                                                                                 
tGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcg                                                                                                                                          <  5:111497/149‑1 (MQ=255)
tGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTca                                                                                                                                          <  2:56494/149‑2 (MQ=255)
      gccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAg                                                                                                                                    >  2:33805/1‑149 (MQ=255)
            gCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCg                                                                                                                              <  1:33805/149‑1 (MQ=255)
              cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGAGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCt                                                                                                                            >  8:211085/1‑149 (MQ=255)
                       ggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTg                                                                                                                    <  2:338656/148‑1 (MQ=255)
                          gATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAg                                                                                                                <  6:319810/149‑1 (MQ=255)
                                     ccATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGt                                                                                                     <  2:207/149‑1 (MQ=255)
                                                           cGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgcgc                                                                               <  6:330156/149‑1 (MQ=255)
                                                           cGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTg                                                                                  <  4:55933/146‑1 (MQ=255)
                                                           cGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTg                                                                                  >  3:55933/1‑146 (MQ=255)
                                                                     acaTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGTTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCa                                                                     <  7:211085/149‑1 (MQ=255)
                                                                          gcttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGt                                                                 <  6:13240/148‑1 (MQ=255)
                                                                              gcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGACTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGaaa                                                            >  1:215479/1‑149 (MQ=255)
                                                                                     cGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagc                                                                                          >  3:292235/1‑112 (MQ=255)
                                                                                     cGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagc                                                                                          <  4:292235/112‑1 (MQ=255)
                                                                                      ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                    >  5:348018/1‑149 (MQ=255)
                                                                                      ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatg                                                                                                                                        <  1:205022/65‑4 (MQ=38)
                                                                                      ggCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatg                                                                                                                                        >  2:205022/1‑62 (MQ=38)
                                                                                        ccATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                    >  7:221247/1‑147 (MQ=255)
                                                                                          aTTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAgg                                                <  8:221247/149‑1 (MQ=255)
                                                                                            tttCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCt                                              >  4:140816/1‑149 (MQ=255)
                                                                                                 gcgcCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCt                                         >  7:296205/1‑149 (MQ=255)
                                                                                                       cTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGAt                                     >  6:37000/1‑147 (MQ=255)
                                                                                                              gCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTgg                                                        >  6:248257/1‑121 (MQ=255)
                                                                                                               cTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAgg                           <  3:140816/149‑1 (MQ=255)
                                                                                                                tGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGc                          >  5:309084/1‑149 (MQ=255)
                                                                                                                        gtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtg                  <  8:296205/149‑1 (MQ=255)
                                                                                                                            taccgcgAAGGCGTGAACGCCTTCACGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggct              >  7:111952/7‑149 (MQ=255)
                                                                                                                             ttCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctg             >  8:270969/1‑149 (MQ=255)
                                                                                                                                 cgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGcc                         >  1:20895/2‑133 (MQ=255)
                                                                                                                                 cgAAGGCGTGAACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGcc                         <  2:20895/132‑1 (MQ=255)
                                                                                                                                    aGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCt                                              >  6:161928/1‑109 (MQ=255)
                                                                                                                                    aGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCt                                              <  5:161928/109‑1 (MQ=255)
                                                                                                                                     ggCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGa     <  2:215479/149‑1 (MQ=255)
                                                                                                                                      gCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAg    <  7:270969/149‑1 (MQ=255)
                                                                                                                                      gCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAg    <  8:101485/149‑1 (MQ=255)
                                                                                                                                        gTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGt  >  4:329315/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCCTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGT  >  NC_002947/3717707‑3717991

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