Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,739,325 C→T 8.1% V195I (GTT→ATT)  yieH ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,739,3250CT 8.1% 272.6 / 10.5 86V195I (GTT→ATT) yieHpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (41/38);  new base T (2/5);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

CCAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATTGT  >  AP009048/3739257‑3739392
                                                                    |                                                                   
ccAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACg                                                                    <  1:2324548/70‑1 (MQ=255)
ccAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACg                                                                    <  1:4066550/70‑1 (MQ=255)
ccAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGa                                                                   >  1:2299128/1‑71 (MQ=255)
   aCCACGCGCTTTCCACAGTTCATGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCg                                                                <  1:1614848/71‑1 (MQ=255)
   aCCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCg                                                                <  1:2673052/71‑1 (MQ=255)
      aCGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCt                                                             >  1:1909796/1‑71 (MQ=255)
        gcgcTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTa                                                           <  1:2808983/71‑1 (MQ=255)
         cgcTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTAt                                                          >  1:2449091/1‑71 (MQ=255)
           cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGc                                                              <  1:12400/65‑1 (MQ=255)
           cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATtg                                                        <  1:4116304/71‑1 (MQ=255)
           cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATtg                                                        >  1:2241519/1‑71 (MQ=255)
             ttCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATtgt                                                       >  1:523770/1‑70 (MQ=255)
             ttCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGATTATtgt                                                       >  1:2977812/1‑70 (MQ=255)
              tCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAGCGAt                                                                  >  1:1049049/1‑58 (MQ=255)
              tCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGAt                                                                  >  1:1121660/1‑58 (MQ=255)
               ccACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGc                                                     <  1:3275537/70‑1 (MQ=255)
                  caGTTGAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                 <  1:4334363/71‑1 (MQ=255)
                  caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGAt                                                                  >  1:2931608/1‑54 (MQ=255)
                  caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATc                                                                 >  1:3445892/1‑55 (MQ=255)
                  caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                 >  1:3360851/1‑71 (MQ=255)
                  caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                 <  1:1553929/71‑1 (MQ=255)
                  caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                 <  1:1032745/71‑1 (MQ=255)
                   aGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTc                                                >  1:403609/1‑71 (MQ=255)
                    gTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCg                                               <  1:4172517/71‑1 (MQ=255)
                    gTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCg                                               <  1:352902/71‑1 (MQ=255)
                     ttCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCgg                                              >  1:3581553/1‑71 (MQ=255)
                       cAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGgcg                                            >  1:161379/1‑71 (MQ=255)
                       cAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGgcg                                            >  1:3225064/1‑71 (MQ=255)
                        aGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGa                                                                   <  1:1708390/47‑1 (MQ=255)
                          gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGgcg                                            >  1:3054541/1‑68 (MQ=255)
                             aCTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAg                                         >  1:2320950/1‑68 (MQ=255)
                             aCTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGaa                                       <  1:3449386/70‑1 (MQ=255)
                               tGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCtt                                                            <  1:341405/47‑1 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAg                                      >  1:2777241/1‑67 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGc                                                              >  1:1710413/1‑43 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGc                                                              >  1:1085144/1‑43 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                     >  1:182736/1‑68 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTa                                    <  1:993431/69‑1 (MQ=255)
                                 cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                   >  1:806259/1‑70 (MQ=255)
                                  gAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGaa                                 <  1:2490934/71‑1 (MQ=255)
                                      gATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACt                             <  1:1636012/71‑1 (MQ=255)
                                       aTGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACtt                            >  1:970375/1‑71 (MQ=255)
                                       aTGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACtt                            <  1:3934194/71‑1 (MQ=255)
                                        tGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                           >  1:3175794/1‑71 (MQ=255)
                                        tGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                           >  1:3181313/1‑71 (MQ=255)
                                          ggTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACtt                            <  1:2458669/68‑1 (MQ=255)
                                           gTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCa                         <  1:972416/70‑1 (MQ=255)
                                           gTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                        <  1:3330383/71‑1 (MQ=255)
                                            tAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                        <  1:3449755/70‑1 (MQ=255)
                                                ggtggtGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGGCGGCGCAGAAGTAGAACACTTCCATACCTg                   <  1:1201073/71‑1 (MQ=255)
                                                 gtggtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGc                  <  1:2787785/71‑1 (MQ=255)
                                                   ggtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                        <  1:2108758/63‑1 (MQ=255)
                                                   ggtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                        >  1:762872/1‑63 (MQ=255)
                                                   ggtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGt                >  1:3605533/1‑71 (MQ=255)
                                                    gtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                   <  1:2890720/51‑1 (MQ=255)
                                                    gtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTg                   <  1:3064265/67‑1 (MQ=255)
                                                     tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                   <  1:3966954/50‑1 (MQ=255)
                                                     tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc               >  1:670186/1‑70 (MQ=255)
                                                     tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc               >  1:2975815/1‑70 (MQ=255)
                                                     tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCg              <  1:902477/71‑1 (MQ=255)
                                                      gACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGaa                                 >  1:3528840/1‑51 (MQ=255)
                                                      gACTTTCGGGTGAACGATCGGCTAATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGa             >  1:3030935/1‑71 (MQ=255)
                                                       aCTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGaa                                 >  1:4197121/1‑50 (MQ=255)
                                                        cTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACtt                            <  1:488866/54‑1 (MQ=255)
                                                        cTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACtt                            <  1:1608992/54‑1 (MQ=255)
                                                        cTTTCGGGTGAACGATCTGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATa           >  1:1967572/1‑71 (MQ=255)
                                                        cTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGAt            >  1:3419682/1‑70 (MQ=255)
                                                         tttCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                           <  1:3036458/54‑1 (MQ=255)
                                                          ttCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                   <  1:3397608/45‑1 (MQ=255)
                                                           tCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc          >  1:4299961/1‑69 (MQ=255)
                                                           tCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc          >  1:3421164/1‑69 (MQ=255)
                                                            cGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                     >  1:309931/1‑41 (MQ=255)
                                                            cGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                     >  1:575626/1‑41 (MQ=255)
                                                            cGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc          >  1:4108503/1‑68 (MQ=255)
                                                             gggTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGa      <  1:2711334/71‑1 (MQ=255)
                                                              ggTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGa      <  1:786392/70‑1 (MQ=255)
                                                              ggTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGAt     <  1:3755442/71‑1 (MQ=255)
                                                              ggTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGGCGATACCAGa      <  1:2268402/70‑1 (MQ=255)
                                                                tGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGAt     <  1:3112036/69‑1 (MQ=255)
                                                                tGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtg   >  1:4000101/1‑71 (MQ=255)
                                                                 gAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                     >  1:3793357/1‑36 (MQ=255)
                                                                 gAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTcc                          >  1:305184/1‑47 (MQ=255)
                                                                  aaCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                   <  1:1020322/37‑1 (MQ=255)
                                                                  aaCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtgt  <  1:4191809/70‑1 (MQ=255)
                                                                  aaCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtg   >  1:3980825/1‑69 (MQ=255)
                                                                   aCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAacac                              <  1:4158870/41‑1 (MQ=255)
                                                                    |                                                                   
CCAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATTGT  >  AP009048/3739257‑3739392

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: ‑