BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	572,252	A→G	29.5%	A165A (GCT→GCC)	aqpZ ←	aquaporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	572,252	0	A	G	29.5%	48.9 / 22.1	54	A165A (GCT→GCC)	aqpZ	aquaporin
Reads supporting (aligned to +/- strand): ref base A (24/14); new base G (13/3); total (37/17)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAACCTGC  >  minE/572189‑572320

gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                                                                 <  1:2581928/69‑1 (MQ=255)

gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                                                                 <  1:1157276/69‑1 (MQ=255)

gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                                                                 <  1:1428513/69‑1 (MQ=255)

 ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGAt                                                                 >  1:3208379/1‑68 (MQ=255)

 ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                                                >  1:333286/1‑69 (MQ=255)

            gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCAGCGCAAAAcc                                                     >  1:1980053/1‑69 (MQ=255)

                        aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa                                         <  1:2623768/69‑1 (MQ=255)

                        aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa                                         <  1:1435480/69‑1 (MQ=255)

                        aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa                                         <  1:1948636/69‑1 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa                                         >  1:948065/1‑68 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa                                         >  1:3009856/1‑68 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa                                         >  1:2063640/1‑68 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa                                         >  1:1992461/1‑68 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:3360277/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:444566/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:1480765/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:1178516/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:857902/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:2804589/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:2518255/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:1209389/1‑69 (MQ=255)

                         aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt                                        >  1:1186120/1‑69 (MQ=255)

                            cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTg                                     >  1:2918117/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:968939/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:3070121/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:2978187/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:24404/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:2215747/1‑69 (MQ=255)

                                  aaGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGtt                               >  1:1619151/1‑69 (MQ=255)

                                           aGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTgg                      <  1:2581630/69‑1 (MQ=255)

                                           aGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTgg                      <  1:607893/69‑1 (MQ=255)

                                           aGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTgg                      <  1:2956806/69‑1 (MQ=255)

                                           aGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTgg                      <  1:3081152/69‑1 (MQ=255)

                                           aGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTgg                      <  1:3288176/69‑1 (MQ=255)

                                                tAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCa                  >  1:2528056/1‑68 (MQ=255)

                                                tAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCa                  >  1:2894779/1‑68 (MQ=255)

                                                tAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCa                  >  1:1781348/1‑68 (MQ=255)

                                                tAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGAACa                  >  1:1633817/1‑68 (MQ=255)

                                                      cAGACCTATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACa            >  1:439926/1‑68 (MQ=255)

                                                      cAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGcc                          >  1:597574/1‑54 (MQ=255)

                                                      cAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAAc             >  1:2222968/1‑67 (MQ=255)

                                                      cAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACa            >  1:1042165/1‑68 (MQ=255)

                                                      cAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACa            >  1:1685382/1‑68 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaa        <  1:2304519/69‑1 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaa        <  1:772804/69‑1 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaa        <  1:778045/69‑1 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaa        <  1:1338563/69‑1 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaa        <  1:2693000/69‑1 (MQ=255)

                                                         aCCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGGTGCGCCGTGGATCACCAACAGGaa        <  1:909248/69‑1 (MQ=255)

                                                          ccAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGTaaa       >  1:2183793/1‑69 (MQ=255)

                                                          ccAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGaaa       >  1:3258252/1‑69 (MQ=255)

                                                             aTAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAAcc     >  1:3139111/1‑68 (MQ=255)

                                                             aTAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAAcc     >  1:90364/1‑68 (MQ=255)

                                                             aTAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAAc      >  1:1071822/1‑67 (MQ=255)

                                                               aGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAACCTGc  >  1:3008202/1‑69 (MQ=255)

GTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGATCACCAACAGGAAACCTGC  >  minE/572189‑572320

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