Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,737,004 A→G 15.0% T202A (ACC→GCC)  pheA → fused chorismate mutase P and prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,737,0040AG15.0% 89.7 / 9.0 40T202A (ACC→GCC) pheAfused chorismate mutase P and prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (2/4);  total (17/23)
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.65e-01

TCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGTTATCCGC  >  W3110S.gb/2736935‑2737073
                                                                     |                                                                     
tCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATccac                                                                      >  1:3185358/1‑71 (MQ=255)
   aTTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCAt                                                                   >  1:3817285/1‑71 (MQ=255)
   aTTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCAt                                                                   >  1:3259089/1‑71 (MQ=255)
    ttgttgGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATcca                                                                       >  1:3850085/1‑66 (MQ=255)
    ttgttgGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATcca                                                                       >  1:1755724/1‑66 (MQ=255)
     tgttgGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATc                                                                  >  1:995845/1‑70 (MQ=255)
       ttgGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCaa                                                                <  1:2575714/70‑1 (MQ=255)
         ggCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATa                                                              <  1:1745685/70‑1 (MQ=255)
         ggCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATAc                                                             >  1:316142/1‑71 (MQ=255)
          gCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACg                                                            <  1:2121369/71‑1 (MQ=255)
             agaTGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACg                                                            <  1:3506410/68‑1 (MQ=255)
                tGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATAGGGTCTa                                                       <  1:496262/70‑1 (MQ=255)
                  aCGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAg                                                    >  1:554599/1‑71 (MQ=255)
                         cTATCGATCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCGCCATCAATACGGTCTACAGCCATCCg                                             <  1:1461934/71‑1 (MQ=255)
                         cTATCGATCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCGCCATCAATACGGTCTACAGCCATCCg                                             <  1:2772689/71‑1 (MQ=255)
                                    tgtttgttGGTCTCCGGCACTACTGATTTAACCACCATCAATACGGTCTACAGCCATCCGCAGCCATTc                                    >  1:1091509/1‑69 (MQ=255)
                                      tttgttGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGc                                <  1:3655242/71‑1 (MQ=255)
                                      tttgttGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGc                                <  1:237691/71‑1 (MQ=255)
                                      tttgttGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGc                                <  1:3337272/71‑1 (MQ=255)
                                        tgttgGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCaa                              >  1:1129586/1‑71 (MQ=255)
                                        tgttgGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCaa                              >  1:2489354/1‑71 (MQ=255)
                                         gttgGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCGATTc                                    <  1:669060/64‑1 (MQ=255)
                                          ttgGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGACATTcc                                   <  1:154/64‑1 (MQ=255)
                                            ggTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATg                            >  1:2239408/1‑69 (MQ=255)
                                            ggTCTCCGGCACTACTGATTTATCCACCATCAAAACGGTCTACAGCCATCCGCAGCCATTCCAGCAATg                            >  1:1984502/1‑69 (MQ=255)
                                             gTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATgca                          <  1:619514/70‑1 (MQ=255)
                                                tcCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATgcagca                       >  1:2038798/1‑70 (MQ=255)
                                                  cGGCACTACTGATTTATCCGCCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCaaa                     <  1:1214744/70‑1 (MQ=255)
                                                  cGGCACTACTGATTTATCCGCCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCaaa                     <  1:754917/70‑1 (MQ=255)
                                                  cGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATgcagca                       <  1:3556853/68‑1 (MQ=255)
                                                  cGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAAt                    >  1:3253124/1‑71 (MQ=255)
                                                      actactGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCt                <  1:311427/71‑1 (MQ=255)
                                                        tactGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGc                                         <  1:411611/44‑1 (MQ=255)
                                                         actGATTTATCCGCCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGc                                >  1:3653526/1‑52 (MQ=255)
                                                         actGATTTATCCGCCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGc                                >  1:2272482/1‑52 (MQ=255)
                                                         actGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAAt                    >  1:3614110/1‑64 (MQ=255)
                                                          ctGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTa              <  1:2017596/69‑1 (MQ=255)
                                                          ctGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTa              <  1:783263/69‑1 (MQ=255)
                                                          ctGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAAt            >  1:2837325/1‑71 (MQ=255)
                                                             aTTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGt         <  1:3824836/71‑1 (MQ=255)
                                                             aTTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGt         <  1:2099481/71‑1 (MQ=255)
                                                              tttATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGtt        <  1:2838378/71‑1 (MQ=255)
                                                                tATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGTTAt      >  1:2168320/1‑71 (MQ=255)
                                                                    caccaTCAATACGGTCTACAGCCATCCGCAGCCATTcc                                   <  1:606650/38‑1 (MQ=255)
                                                                    caccaTCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGTTATCCGc  <  1:644507/71‑1 (MQ=255)
                                                                     |                                                                     
TCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTACAGCCATCCGCAGCCATTCCAGCAATGCAGCAAATTCCTTAATCGTTATCCGC  >  W3110S.gb/2736935‑2737073

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