Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 349,140 C→T 12.5% intergenic (+344/‑96) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048349,1400CT12.5% 132.9 / 9.5 48intergenic (+344/‑96)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (4/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTA  >  AP009048/349080‑349209
                                                            |                                                                     
aCGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtga                                                             <  1:3993243/71‑1 (MQ=31)
  gCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCTTACGGGGtgagc                                                           >  1:71033/1‑69 (MQ=9)
      aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactgaa                                                       >  1:2523963/1‑71 (MQ=34)
         gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGaactgaactg                                                    <  1:4599333/71‑1 (MQ=34)
            gCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAg                                                 >  1:3944498/1‑71 (MQ=34)
             cATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg                                                <  1:1352971/71‑1 (MQ=34)
                  gACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGa                                           <  1:4544203/71‑1 (MQ=34)
                    cAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGAt                                          >  1:1589325/1‑70 (MQ=255)
                    cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATa                                         <  1:76233/71‑1 (MQ=34)
                    cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATa                                         <  1:2359312/71‑1 (MQ=34)
                       ccGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGa                                      <  1:4238364/71‑1 (MQ=34)
                         gATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtga                                                             >  1:545251/1‑46 (MQ=32)
                           tGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCa                                  <  1:2215149/71‑1 (MQ=34)
                           tGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCa                                  >  1:3591795/1‑71 (MQ=34)
                            gCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaact                                                          >  1:3073493/1‑46 (MQ=38)
                             ccTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATa                                <  1:3615339/71‑1 (MQ=34)
                             ccTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATa                                <  1:1153485/71‑1 (MQ=34)
                               tGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGc                                   <  1:349169/66‑1 (MQ=35)
                               tGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATa                                <  1:1011521/69‑1 (MQ=35)
                                gATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg                             >  1:4398223/1‑71 (MQ=34)
                                 aTGCGCCGCTGACGCGACTTATCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg                             <  1:1691203/70‑1 (MQ=25)
                                 aTGCGCCGCTGACGCGACTTATCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGt                            >  1:4180688/1‑71 (MQ=21)
                                 aTGCGCCGCTGACGCGACTTATCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGt                            >  1:3549995/1‑71 (MQ=21)
                                 aTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAg                                                 >  1:3268360/1‑50 (MQ=38)
                                 aTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg                             <  1:1856628/70‑1 (MQ=35)
                                   gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt                                                   <  1:3557227/46‑1 (MQ=38)
                                   gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                          >  1:1776299/1‑71 (MQ=34)
                                     gccgcTGACGCGACTTATCAGGCCTACGGGGtgaactg                                                         <  1:532251/38‑1 (MQ=40)
                                     gccgcTGACGCGACTTATCAGGCCTACGGGGtgaactg                                                         <  1:1047033/38‑1 (MQ=40)
                                     gccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                        <  1:2023738/71‑1 (MQ=34)
                                      ccgcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGa                                           >  1:3148038/1‑51 (MQ=35)
                                        gcTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                     <  1:1845842/71‑1 (MQ=34)
                                         cTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                     >  1:2917295/1‑70 (MQ=35)
                                           gACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                    <  1:3843111/69‑1 (MQ=35)
                                           gACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                   >  1:4593603/1‑70 (MQ=35)
                                            aCGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac                  <  1:1254062/70‑1 (MQ=35)
                                              gcgACTTATCAGGCCTACGGGGTGAACTGATCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGacaa                >  1:3410956/1‑70 (MQ=25)
                                                gACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac                  <  1:2345099/66‑1 (MQ=35)
                                                gACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAcaacaa             >  1:4257139/1‑71 (MQ=38)
                                                gACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAcaacaa             >  1:2067432/1‑71 (MQ=38)
                                                  cTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAtc           >  1:4483802/1‑71 (MQ=40)
                                                   ttATCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg                             <  1:2857563/52‑1 (MQ=255)
                                                   ttATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAtct          >  1:3486763/1‑71 (MQ=255)
                                                   ttATCAGACCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAtct          >  1:2505537/1‑71 (MQ=255)
                                                     aTCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                      >  1:528617/1‑57 (MQ=255)
                                                     aTCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAtctc         >  1:392334/1‑70 (MQ=255)
                                                      tCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGa       <  1:734680/71‑1 (MQ=255)
                                                        aGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGAcc     <  1:3152602/71‑1 (MQ=255)
                                                            cTACGGGTTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTa  >  1:35022/1‑70 (MQ=255)
                                                            cTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTa  >  1:1261690/1‑70 (MQ=255)
                                                            |                                                                     
ACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTA  >  AP009048/349080‑349209

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