BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	572,219	A→G	13.4%	I176I (ATT→ATC)	aqpZ ←	aquaporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	572,219	0	A	G	13.4%	266.9 / 29.1	104	I176I (ATT→ATC)	aqpZ	aquaporin
Reads supporting (aligned to +/- strand): ref base A (39/51); new base G (2/12); total (41/63)
Fisher's exact test for biased strand distribution p-value = 4.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 104 total aligned reads displayed.

CGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCG  >  minE/572152‑572288

cgccCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAAtt                                                                    <  1:555417/71‑1 (MQ=255)

       gAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGaa                                                             >  1:2978702/1‑71 (MQ=255)

       gAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGaa                                                             >  1:1462932/1‑71 (MQ=255)

       gAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGaa                                                             >  1:4474356/1‑71 (MQ=255)

       gAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGa                                                              <  1:2866329/70‑1 (MQ=255)

           aTAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAg                                                         <  1:1142657/71‑1 (MQ=255)

            tAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAg                                                         <  1:1038247/70‑1 (MQ=255)

             aGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGt                                                                >  1:5013871/1‑62 (MQ=255)

             aGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGt                                                                >  1:2569956/1‑62 (MQ=255)

             aGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGt                                                                >  1:2576055/1‑62 (MQ=255)

             aGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAggg                                                       >  1:5042732/1‑71 (MQ=255)

                    gcgGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAgg                                                 >  1:3193557/1‑70 (MQ=255)

                    gcgGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                >  1:3015408/1‑71 (MQ=255)

                      ggTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCa                                              <  1:206406/71‑1 (MQ=255)

                      ggTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCa                                              <  1:2610645/71‑1 (MQ=255)

                       gTGCTGCGCGCCGGATTGACCGAGGTATTCGTCACCGGGATGCTGATCAGGTGCACCAGCGCCAGGGCCAg                                             >  1:844260/1‑71 (MQ=255)

                       gTGCTGCGCGCCGGATTGACCGAGGTATTCGTCACCGGGATGCTGATCAGGTGCACCAGCGCCAGGGCCAg                                             >  1:4660769/1‑71 (MQ=255)

                        tactgcGCTCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:2173008/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:2965005/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:292503/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:892411/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:470284/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:3287690/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:3770345/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:2003067/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:1619369/60‑1 (MQ=255)

                        tactgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:1913942/60‑1 (MQ=255)

                         actgcGCGCCGGGTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGtt                                                     <  1:19050/60‑1 (MQ=255)

                          ctgcGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGAcc                                          <  1:4189544/71‑1 (MQ=255)

                          ctgcGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGAcc                                          <  1:772363/71‑1 (MQ=255)

                                 ccGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:4939677/71‑1 (MQ=255)

                                 ccGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:4460754/71‑1 (MQ=255)

                                 ccGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:2272864/71‑1 (MQ=255)

                                 ccGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:4708304/71‑1 (MQ=255)

                                 ccGGGTTAACAGAAGTGTTAGTAACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:1982989/71‑1 (MQ=255)

                                  caggTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:3185087/68‑1 (MQ=255)

                                      ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAg                                                                 >  1:4600797/1‑36 (MQ=255)

                                      ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   >  1:4783111/1‑66 (MQ=255)

                                      ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   >  1:2518264/1‑66 (MQ=255)

                                       tAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        >  1:4682697/1‑60 (MQ=255)

                                       tAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        >  1:3844650/1‑60 (MQ=255)

                                       tAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        >  1:103858/1‑60 (MQ=255)

                                         aCAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGAcc                                          <  1:1963094/56‑1 (MQ=255)

                                         aCAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGAcc                                          <  1:2533295/56‑1 (MQ=255)

                                         aCAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGAcc                                          <  1:2252588/56‑1 (MQ=255)

                                           aGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaa                         <  1:548209/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATTAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:351908/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:1624997/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:4971354/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:5006161/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:2279086/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:1112363/71‑1 (MQ=255)

                                            gAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:2791795/71‑1 (MQ=255)

                                               gtgtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCaaa                        <  1:1951717/68‑1 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTa                                                    >  1:4817929/1‑38 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAgg                                                 >  1:852255/1‑41 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                >  1:1092873/1‑42 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                >  1:1193132/1‑42 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                >  1:1986513/1‑42 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:2677110/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:3276812/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:3697942/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:1250000/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:3773099/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                      >  1:2490054/1‑68 (MQ=255)

                                                 gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGAGATCGGCGCAAAAc                      >  1:101099/1‑68 (MQ=255)

                                                    aGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCg                >  1:1174504/1‑71 (MQ=255)

                                                    aGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCg                >  1:256318/1‑71 (MQ=255)

                                                    aGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCg                >  1:3475265/1‑71 (MQ=255)

                                                    aGTCACCGGAATACTAATTAAGTGAATAAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCg                >  1:4054169/1‑71 (MQ=255)

                                                      tCACCGGAATACTAATTTAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:236711/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                <  1:3447665/37‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:2980214/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:4906828/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:3283910/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:4722915/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:1862375/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:4414726/45‑1 (MQ=255)

                                                      tCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCaa                                        <  1:2089923/45‑1 (MQ=255)

                                                            gAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:4121935/44‑1 (MQ=255)

                                                            gAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa         >  1:3951533/1‑70 (MQ=255)

                                                             aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCg                                   <  1:4134249/43‑1 (MQ=255)

                                                             aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAAt        >  1:3095504/1‑70 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGTCCAGACCAATAGCGATc                                <  1:3905372/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:2727638/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:592444/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:5011769/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:4057109/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:4858903/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:1327428/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:1226682/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:1256512/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                <  1:42970/42‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc   <  1:2995326/71‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc   <  1:5142811/71‑1 (MQ=255)

                                                                 cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc   <  1:2322866/71‑1 (MQ=255)

                                                                  tAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCg  >  1:1878883/1‑71 (MQ=255)

                                                                  tAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCg  >  1:4629350/1‑71 (MQ=255)

                                                                  tAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCg  >  1:3502913/1‑71 (MQ=255)

                                                                  tAATTAAGTGAATAAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCg  >  1:3176984/1‑71 (MQ=255)

CGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCG  >  minE/572152‑572288

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