BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	139,696	T→A	19.0%	V157V (GTA→GTT) ‡	htrE ←	predicted outer membrane usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	139,696	0	T	A	19.0%	237.6 / 44.0	121	V157V (GTA→GTT) ‡	htrE	predicted outer membrane usher protein
Reads supporting (aligned to +/- strand): ref base T (62/36); new base A (10/13); total (72/49)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAA  >  minE/139637‑139756

gCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGc                                                      >  1:2432812/1‑68 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1281586/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1624846/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:2823959/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:2907534/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:2946808/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:2232853/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1677882/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               <  1:2070033/69‑1 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:36701/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1213071/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1180834/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:1905314/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               >  1:2347335/1‑69 (MQ=255)

      ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                               <  1:100206/69‑1 (MQ=255)

        aaTTCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:2277998/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAg                                                        >  1:1267584/1‑58 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTg                                                   >  1:2939106/1‑63 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                              >  1:1844352/1‑68 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                              >  1:3014777/1‑68 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:2285136/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:2611467/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:2443457/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:1627833/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:2376289/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:3224452/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:1364024/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:62521/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:654654/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:90514/1‑69 (MQ=255)

        aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                             >  1:1027762/1‑69 (MQ=255)

           gCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTCGAGGAACGTc                                          <  1:1053444/69‑1 (MQ=255)

            ccGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTc                                          >  1:933159/1‑68 (MQ=255)

                       tAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                              <  1:1644478/69‑1 (MQ=255)

                       tAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                              <  1:1648354/69‑1 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCa                                  >  1:2778558/1‑63 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            <  1:2229606/69‑1 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:2285845/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:2221525/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:2900362/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:3056430/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:3239482/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:3416418/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:871359/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:1346478/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:1340056/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:1603957/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:1947245/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                            >  1:1816690/1‑69 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                             >  1:2103966/1‑68 (MQ=255)

                         aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                             >  1:1158495/1‑68 (MQ=255)

                         aCGATGGATAAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                             >  1:2943867/1‑68 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:2735954/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:2649183/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1818288/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1512616/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1930711/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1932299/68‑1 (MQ=255)

                           gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:3380809/68‑1 (MQ=255)

                            aTGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1819848/67‑1 (MQ=255)

                            aTGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1424009/67‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:2216036/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:3035118/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:60607/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:3301159/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:3286438/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:2654917/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:78405/65‑1 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATc                        >  1:1038818/1‑68 (MQ=255)

                              ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATc                        >  1:138425/1‑68 (MQ=255)

                                 tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                              >  1:1737039/1‑59 (MQ=255)

                                 tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                              >  1:778831/1‑59 (MQ=255)

                                 tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGtt                     <  1:2715569/68‑1 (MQ=255)

                                 tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGtt                     <  1:78475/68‑1 (MQ=255)

                                  cAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGtt                     <  1:1273878/67‑1 (MQ=255)

                                    aCATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                           <  1:1129147/59‑1 (MQ=255)

                                      aTAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACg                 >  1:756191/1‑67 (MQ=255)

                                        aGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCAt             <  1:467084/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGTGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:749256/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:1157484/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:1444888/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:1022105/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:3377166/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:205781/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:1542301/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:1607980/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:2716367/69‑1 (MQ=255)

                                         gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:2508012/69‑1 (MQ=255)

                                         gTTTTGATAATTCATCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:2322446/69‑1 (MQ=255)

                                          ttttGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa            <  1:2030272/68‑1 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:2141195/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:2063249/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:866857/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:2365083/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:558154/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:3069268/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:3082853/2‑68 (MQ=255)

                                               atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa       >  1:3187997/2‑68 (MQ=255)

                                                    ttttCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGaa  <  1:2434839/68‑1 (MQ=255)

                                                    ttttCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGaa  <  1:1530033/68‑1 (MQ=255)

GCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAA  >  minE/139637‑139756

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