Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 669,698 T→C 31.2% H333H (CAT→CAC pqiA → paraquat‑inducible membrane protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE669,6980TC31.2% 85.5 / 49.4 99H333H (CAT→CACpqiAparaquat‑inducible membrane protein A
Reads supporting (aligned to +/- strand):  ref base T (20/48);  new base C (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.91e-01

Only 100 of 108 total aligned reads displayed.

CGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGAT  >  minE/669649‑669766
                                                 |                                                                    
cGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGaa                                                    >  1:299452/1‑68 (MQ=255)
 gtcaaaGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:1153909/66‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGTATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:1759031/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:1685469/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:2688021/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:2648239/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:907021/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:801203/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:1020597/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:484482/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:2145731/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:729277/69‑1 (MQ=255)
 gTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAg                                                  <  1:675905/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:1820986/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:2171152/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:1044364/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:49198/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:1076731/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:1545706/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:765944/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:129507/69‑1 (MQ=255)
   tAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGaa                                                <  1:742798/69‑1 (MQ=255)
          gatCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCAttt                                         <  1:2089453/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                                   <  1:159255/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                                   <  1:2886377/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                                   <  1:2887952/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                                   <  1:2911093/69‑1 (MQ=255)
                cATTGCGTGGCTGTGTTGGGATGCCAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTa                                   <  1:1278908/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:1566057/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                               >  1:2808333/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                               >  1:456275/1‑69 (MQ=255)
                    gcgTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                               >  1:1839707/1‑69 (MQ=255)
                    gcctGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGaa                               >  1:1143362/4‑69 (MQ=255)
                         gCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAgttg                           >  1:1909977/1‑68 (MQ=255)
                         gCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAgttg                           >  1:2578675/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:2956013/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt                    >  1:3127929/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt                    >  1:2196040/1‑69 (MQ=255)
                               ttGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGtt                    >  1:2659028/1‑69 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:2550494/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:961843/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:2798498/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:3241600/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:3406148/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:2781074/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:1110045/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:1870055/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:1322682/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:1516604/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:2378695/43‑1 (MQ=255)
                                 gggATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCa                                            <  1:2369080/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:567063/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:3428992/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:526993/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:1887734/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:1271362/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:1055189/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:967392/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:3327417/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:2652719/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:2599393/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:2413448/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:2913422/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:3045213/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:216100/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:3244345/1‑69 (MQ=255)
                                       ccAAAGGACACGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGcc            >  1:1910939/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)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:2167518/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:588887/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:2909986/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:135704/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:742304/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:1990884/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCgat  <  1:2615706/69‑1 (MQ=255)
                                                 tGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTACGCCGCTGGTCgat  <  1:3056112/69‑1 (MQ=255)
                                                 |                                                                    
CGTTAAAGATGATCGCCATCGCGTGGCTGTGTTGGGATGCCAAAGGGCATGGCAAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGAT  >  minE/669649‑669766

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: