Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 669,695 G→A 31.2% G332G (GGG→GGA pqiA → paraquat‑inducible membrane protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE669,6950GA31.2% 92.3 / 47.8 99G332G (GGG→GGApqiAparaquat‑inducible membrane protein A
Reads supporting (aligned to +/- strand):  ref base G (20/48);  new base A (24/7);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 9.97e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCG  >  minE/669649‑669757
                                              |                                                              
cGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGaa                                           >  1:299452/1‑68 (MQ=255)
 gtcaaaGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:1153909/66‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGTATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:1759031/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:1020597/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:1685469/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:2145731/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:2648239/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:2688021/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:484482/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:675905/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:729277/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:801203/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:907021/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCGGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:2783259/69‑1 (MQ=255)
  ttAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:49389/68‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGTCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGat                                       <  1:2231598/69‑2 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:49198/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:742798/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:765944/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:2171152/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:1820986/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:969862/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:1044364/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:1545706/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:129507/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                       <  1:1076731/69‑1 (MQ=255)
          gatCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCAttt                                <  1:2089453/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:1278908/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:159255/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:2911093/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:2887952/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:2886377/69‑1 (MQ=255)
                 attgcgTGGCTGTGTTTGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:57782/65‑1 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGa                              >  1:401248/1‑61 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:1914710/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:1839707/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:2808333/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:3224450/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:2129525/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:1395622/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:456275/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:1566057/1‑69 (MQ=255)
                    gcctGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                      >  1:1143362/4‑69 (MQ=255)
                         gCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAgttg                  >  1:2578675/1‑68 (MQ=255)
                         gCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAgttg                  >  1:1909977/1‑68 (MQ=255)
                           tgtgTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:1406689/69‑1 (MQ=255)
                           tgtgTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:2871893/69‑1 (MQ=255)
                           tgtgTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:2121707/69‑1 (MQ=255)
                            gtgtTGGGATGCCAAAGGGCATGGCAAGCGCTACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:1536223/68‑1 (MQ=255)
                            gtgtTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:2240594/68‑1 (MQ=255)
                            gtgtTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:2415553/68‑1 (MQ=255)
                            gtgtTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:2741853/68‑1 (MQ=255)
                            gtgtTGGGATGCCAAAGGACACGGGAAGCGCGACAGTGAAAGAATGCATTTGATTTa                          <  1:3327731/57‑1 (MQ=255)
                            gtgtCGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGttgttg               <  1:3254826/68‑1 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                         <  1:1685791/39‑1 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt           >  1:2659028/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt           >  1:1749283/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt           >  1:3127929/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt           >  1:2196040/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt           >  1:2956013/1‑69 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2798498/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2369080/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:3406148/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2378695/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:1322682/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:3241600/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:961843/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2781074/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:1870055/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2485556/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:1110045/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:2550494/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                   <  1:1516604/43‑1 (MQ=255)
                                       ccAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:158068/1‑69 (MQ=255)
                                       ccAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1845772/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGa                       >  1:2094960/1‑49 (MQ=39)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:561327/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1055189/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:526993/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1271362/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:967392/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:567063/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:3327417/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:169483/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1663111/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1640263/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1887734/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:1910939/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:216100/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:2413448/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:2599393/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:2652719/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:2913422/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:3045213/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:3244345/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc   >  1:3428992/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGc    >  1:3173473/1‑68 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTAATGAAGTTGTTGAGTTTGTAGGcc   >  1:2483967/1‑69 (MQ=255)
                                       ccAAAGGACACGGAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCg  >  1:1078164/1‑69 (MQ=255)
                                              |                                                              
CGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCG  >  minE/669649‑669757

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: