Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,677,417 T→G 15.1% R171R (CGA→CGC yhgA ← predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,677,4170TG15.1% 166.6 / 31.6 86R171R (CGA→CGCyhgApredicted transposase
Reads supporting (aligned to +/- strand):  ref base T (41/32);  new base G (2/11);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

TGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  Exported/2677349‑2677485
                                                                    |                                                                    
tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:5485432/1‑71 (MQ=255)
tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1750307/1‑71 (MQ=255)
tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1580988/1‑71 (MQ=255)
tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:4488858/1‑71 (MQ=255)
tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGAGCCATACGTCg                                                                    >  1:4893237/1‑71 (MQ=255)
  ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:422764/70‑1 (MQ=255)
  ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:2703192/70‑1 (MQ=255)
  ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:1545122/70‑1 (MQ=255)
   tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:140882/1‑71 (MQ=255)
       ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2235050/1‑70 (MQ=255)
       ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2848944/1‑70 (MQ=255)
       ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:15778/1‑70 (MQ=255)
         aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:4615940/1‑65 (MQ=255)
         aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:2862998/1‑65 (MQ=255)
          atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                           <  1:2602926/68‑1 (MQ=255)
           tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa                                                         >  1:5350984/2‑69 (MQ=255)
             ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:1087263/71‑1 (MQ=255)
             ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:4987377/71‑1 (MQ=255)
                tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCa                                                            >  1:4551637/1‑63 (MQ=255)
                tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                    >  1:2979993/1‑71 (MQ=255)
                 cAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                    >  1:4320673/1‑70 (MQ=255)
                  aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:3072811/1‑70 (MQ=255)
                  aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:3514677/1‑70 (MQ=255)
                  aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:1993096/1‑70 (MQ=255)
                    gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                                                                      <  1:4713278/49‑1 (MQ=255)
                    gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                <  1:2660088/71‑1 (MQ=255)
                       gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:1036944/70‑2 (MQ=255)
                       gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:4776253/70‑2 (MQ=255)
                       gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:4977092/70‑2 (MQ=255)
                       gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   >  1:5460968/1‑49 (MQ=255)
                        cgctGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCgg                                            <  1:5353350/68‑1 (MQ=255)
                          ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                     <  1:886975/44‑1 (MQ=255)
                          ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                          <  1:5202700/71‑1 (MQ=255)
                           tGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         >  1:3925169/1‑71 (MQ=255)
                            gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         >  1:1599291/1‑70 (MQ=255)
                             aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         <  1:921756/69‑1 (MQ=255)
                             aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGa                                       <  1:5650889/71‑1 (MQ=255)
                               gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACg                                     >  1:2702693/1‑71 (MQ=255)
                               gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3023090/50‑2 (MQ=255)
                               gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:5350784/50‑2 (MQ=255)
                               gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3992153/50‑2 (MQ=255)
                                aaTATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 <  1:2917339/42‑1 (MQ=255)
                                 atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:5192344/1‑71 (MQ=255)
                                 atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   <  1:2240814/71‑1 (MQ=255)
                                 atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1906093/1‑71 (MQ=255)
                                 atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1891867/1‑71 (MQ=255)
                                 atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:665995/1‑71 (MQ=255)
                                  tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGa                                                         >  1:1334354/1‑48 (MQ=255)
                                  tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   <  1:889489/70‑1 (MQ=255)
                                    tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                                <  1:389629/71‑1 (MQ=255)
                                    tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                                <  1:1823691/71‑1 (MQ=255)
                                    tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGAATGATTTCATCATCAGGCATGACGGTGAc                                <  1:2451865/71‑1 (MQ=255)
                                     gTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACa                               >  1:1055009/1‑71 (MQ=255)
                                      tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                             <  1:5547609/56‑1 (MQ=255)
                                      tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                             <  1:3642943/56‑1 (MQ=255)
                                       ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                              >  1:2630057/1‑70 (MQ=255)
                                       ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                              >  1:2246251/1‑70 (MQ=255)
                                        tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatc                             <  1:3539959/70‑1 (MQ=255)
                                        tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatc                                               >  1:336148/3‑52 (MQ=255)
                                        tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCCGCACCACGGTAATatc                             <  1:5120915/68‑1 (MQ=255)
                                         ctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatcc                            <  1:1557292/69‑2 (MQ=255)
                                          tGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaat                          >  1:2548112/1‑71 (MQ=255)
                                           gAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatc                             <  1:2017244/67‑1 (MQ=255)
                                            aaTTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaatca                        <  1:5275029/71‑1 (MQ=255)
                                               tAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                >  1:3982781/1‑44 (MQ=255)
                                                aaCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCg                     >  1:4759618/1‑70 (MQ=255)
                                                 aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                   <  1:2503480/71‑1 (MQ=255)
                                                   tCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaa                  >  1:3936979/1‑70 (MQ=255)
                                                    ccAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAAt                >  1:3397653/1‑71 (MQ=255)
                                                     cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaaa                 >  1:4642543/1‑69 (MQ=255)
                                                     cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:3098297/1‑71 (MQ=255)
                                                     cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:391025/1‑71 (MQ=255)
                                                      agcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:4353244/1‑70 (MQ=255)
                                                        cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCgg            <  1:4105215/71‑1 (MQ=255)
                                                        cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg             >  1:4473905/1‑70 (MQ=255)
                                                         agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:494097/71‑1 (MQ=255)
                                                         agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:4718488/71‑1 (MQ=255)
                                                         agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:3833914/71‑1 (MQ=255)
                                                         agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:3949063/71‑1 (MQ=255)
                                                           cgcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg         <  1:4749342/71‑1 (MQ=255)
                                                            gcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGc        <  1:4308832/71‑1 (MQ=255)
                                                               aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa     <  1:5114754/71‑1 (MQ=255)
                                                               aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa     <  1:1096705/71‑1 (MQ=255)
                                                                tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa    <  1:1019130/71‑1 (MQ=255)
                                                                tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa    <  1:5646981/71‑1 (MQ=255)
                                                                  cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:1007875/1‑71 (MQ=255)
                                                                  cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:4786541/1‑71 (MQ=255)
                                                                    |                                                                    
TGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  Exported/2677349‑2677485

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: