BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	785,969	G→A	33.2%	G321G (GGC→GGT)	nhaB ←	sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	785,969	0	G	A	33.2%	91.7 / 56.6	108	G321G (GGC→GGT)	nhaB	sodium:proton antiporter
Reads supporting (aligned to +/- strand): ref base G (44/28); new base A (22/14); total (66/42)
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 108 total aligned reads displayed.

GATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACA  >  minE/785903‑786032

gATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               <  1:2492116/69‑1 (MQ=255)

gATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               <  1:3072962/69‑1 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2810495/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2850374/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2954908/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2218113/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2999000/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               <  1:543865/60‑1 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2795945/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2191516/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:665642/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:665933/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:888233/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:1911961/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:956264/1‑60 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTcc                                                        >  1:3189043/1‑67 (MQ=255)

         cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTcc                                                        >  1:3037841/1‑67 (MQ=255)

               ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                 >  1:253054/1‑68 (MQ=255)

               ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                 >  1:2031303/1‑68 (MQ=255)

               ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                 >  1:1598776/1‑68 (MQ=255)

               ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                 >  1:2193250/1‑68 (MQ=255)

               ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                 >  1:920084/1‑68 (MQ=255)

                 tGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:2492870/68‑1 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2260600/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2145271/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2151684/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:2184387/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:227017/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:1848316/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:664365/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:876952/1‑42 (MQ=255)

                           cAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                               >  1:1035523/1‑42 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCTACTTCCGCCAGATGc                                               <  1:1759469/54‑1 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:2238912/54‑1 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:549189/54‑1 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:956952/54‑1 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:1236422/54‑1 (MQ=255)

                               gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                               <  1:1807485/54‑1 (MQ=255)

                                  gTTGCCAGAATAATGACTGACAAACCAATCAAGGCAACTTCCGCCAGATGc                                               <  1:2793286/51‑1 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:3084568/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:1019104/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:1167457/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:878528/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:726204/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:714785/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:703998/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:1387481/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:1430455/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:518958/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:377684/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:3409151/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:1553080/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:2757496/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                        >  1:2325073/1‑69 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                          >  1:2070431/1‑67 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                          >  1:929546/1‑67 (MQ=255)

                                       cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGACAGCGCAGTCAccagcca                        >  1:1340766/1‑69 (MQ=255)

                                        aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                        <  1:2974062/68‑1 (MQ=255)

                                        aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                        <  1:2766647/68‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:2022944/69‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:1386504/69‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:1866549/69‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:2034382/69‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:657516/69‑1 (MQ=255)

                                             aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:409368/69‑1 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:793690/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:463654/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  <  1:1390341/68‑1 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:1087916/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:3189941/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:1600856/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:2885790/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:1634500/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGAAACGCCAGCGCAGTCACCAGCCAGACGcc                  >  1:1447334/1‑68 (MQ=255)

                                              atGACTGACAAACCAATCAAGCAAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCAAGACGcc                  >  1:457354/1‑68 (MQ=255)

                                                 actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata              <  1:2809468/69‑1 (MQ=255)

                                                 actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata              <  1:1460705/69‑1 (MQ=255)

                                                 actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata              <  1:2940257/69‑1 (MQ=255)

                                                 actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata              <  1:891224/69‑1 (MQ=255)

                                                 actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata              <  1:565500/69‑1 (MQ=255)

                                                      cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc         >  1:2407936/1‑69 (MQ=255)

                                                      cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc         >  1:1351573/1‑69 (MQ=255)

                                                      cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc         >  1:2007468/1‑69 (MQ=255)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:2039667/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:3024253/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:2061797/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:1755950/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:2794961/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:2989402/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:566585/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:556023/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:1541287/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:539750/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:21044/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:1929087/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:251461/54‑1 (MQ=25)

                                                         aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                     <  1:3133974/54‑1 (MQ=25)

                                                                aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACa  >  1:2281766/1‑66 (MQ=255)

                                                                aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACa  >  1:3028986/1‑66 (MQ=255)

                                                                aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGAACa  >  1:2284462/1‑66 (MQ=255)

GATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACA  >  minE/785903‑786032

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