Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 131,770 G→A 27.3% E21E (GAG→GAA yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE131,7700GA27.3% 142.3 / 49.2 124E21E (GAG→GAAyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base A (18/16);  total (64/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

ATGGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACA  >  minE/131708‑131834
                                                              |                                                                
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:1939107/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:1845210/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:3278095/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:1213602/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:2088492/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:1231873/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:2541823/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:528409/69‑1 (MQ=255)
atggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:2627549/69‑1 (MQ=255)
atggatgCACCGAATACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:3240228/69‑1 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:209955/68‑1 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:272638/68‑1 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGAGGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:3173365/1‑69 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:1024726/1‑69 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:1304315/1‑69 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:1297278/1‑69 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:888217/1‑69 (MQ=255)
 tggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTc                                                           >  1:815963/1‑69 (MQ=255)
          cGAGTACCACACCGCATGATGCGGTGTTTAAACAATTTTTAATGCATGCGGAAACGGCt                                                            <  1:1227878/59‑1 (MQ=255)
              tACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGga                                              <  1:2557168/69‑1 (MQ=255)
              tACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGga                                              <  1:1114132/69‑1 (MQ=255)
                ccACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGgaga                                            >  1:2531890/1‑69 (MQ=255)
                ccACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGgaga                                            >  1:2301120/1‑69 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:1260283/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:3207648/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:3249280/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:3317996/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:3201596/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:2174412/69‑1 (MQ=255)
                  acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:2504371/69‑1 (MQ=255)
                   cacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATAc                                         <  1:187021/69‑1 (MQ=255)
                              gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTc                                                  >  1:2987223/1‑49 (MQ=255)
                              gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATa                                          <  1:2072164/57‑1 (MQ=255)
                              gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAg                                >  1:2955972/1‑67 (MQ=255)
                              gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGt                               >  1:1291631/1‑68 (MQ=255)
                              gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGt                               >  1:1712509/1‑68 (MQ=255)
                              gcgGTATTAAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGt                               >  1:2598818/1‑68 (MQ=255)
                               cgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:417042/69‑1 (MQ=255)
                               cgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:487169/69‑1 (MQ=255)
                                  tATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGga                                              <  1:1774922/49‑1 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACa                                        >  1:2880286/1‑53 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACa                                        >  1:857179/1‑53 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACa                                        >  1:2681327/1‑53 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACa                                        >  1:1448731/1‑53 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACa                                        >  1:1769960/1‑53 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTa                        >  1:548356/1‑69 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTa                        >  1:640944/1‑69 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTa                        >  1:3422/1‑69 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTa                        >  1:742827/1‑69 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTa                        >  1:1108568/1‑69 (MQ=255)
                                    tttAAACAATTTTTAATGCATGCGGAGACGGCTAGCGACTTTCTGGAGATACa                                        >  1:3001233/1‑53 (MQ=255)
                                      tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCTAGTGGAATTAc                       >  1:1394498/1‑68 (MQ=255)
                                      tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTAc                       >  1:1917820/1‑68 (MQ=255)
                                      tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTAc                       >  1:703078/1‑68 (MQ=255)
                                      tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTAc                       >  1:729324/1‑68 (MQ=255)
                                      tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTAc                       >  1:1719027/1‑68 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTAcgcg                    >  1:1049712/1‑68 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:2528525/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:925206/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:1116148/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:843757/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:1389845/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:1730922/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:183090/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:1993921/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:2082546/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:3148557/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:2533567/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:2242715/1‑69 (MQ=255)
                                         aCAATTTTTAATGCATGCGGAAACGGCTCGCGACTTTCTGGATATCCATTTGCCAGCGGAGCTACGCGa                   >  1:242290/1‑69 (MQ=255)
                                          cAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGaa                  <  1:1649389/69‑1 (MQ=255)
                                          cAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGaa                  <  1:3028309/69‑1 (MQ=255)
                                          cAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGaa                  <  1:1733937/69‑1 (MQ=255)
                                          cAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGaa                  <  1:1975085/69‑1 (MQ=255)
                                          cAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGaa                  <  1:2722997/69‑1 (MQ=255)
                                                   atgcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGa          <  1:2123194/68‑1 (MQ=255)
                                                   atgcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGa          <  1:2761824/68‑1 (MQ=255)
                                                   atgcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGa          <  1:3328158/68‑1 (MQ=255)
                                                    tgcatgcTGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGa          <  1:545258/67‑1 (MQ=255)
                                                     gcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGa          <  1:2917277/66‑1 (MQ=255)
                                                      catgcGGAGACGTCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:1360801/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCTCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:1294690/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:2213473/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:3430250/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:2834071/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:918231/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:167058/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:554332/46‑1 (MQ=255)
                                                      catgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:1610151/46‑1 (MQ=255)
                                                         gcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTgg                             <  1:3010223/43‑1 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAAc                 >  1:178153/1‑54 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:680082/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:1672151/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:790539/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:1820682/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:3250290/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTtgtg           >  1:2155199/1‑60 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGCCCTCAaca  >  1:2122064/1‑69 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAaca  >  1:154800/1‑69 (MQ=255)
                                                          cGGAGACGGCTCGCGACTTTCTCGAGATACATTTGCCa                                 >  1:170302/1‑38 (MQ=255)
                                                              |                                                                
ATGGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACA  >  minE/131708‑131834

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: