BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,173,752	G→A	24.8%	V70I (GTT→ATT)	yrbL →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,173,752	0	G	A	24.8%	163.2 / 45.3	128	V70I (GTT→ATT)	yrbL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (53/43); new base A (24/8); total (77/51)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 128 total aligned reads displayed.

ACTACGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTC  >  minE/2173690‑2173816

actacGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTgaga                                                            >  1:496365/1‑69 (MQ=255)

actacGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTgaga                                                            >  1:2490651/1‑69 (MQ=255)

actacGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTgaga                                                            >  1:1166081/1‑69 (MQ=255)

actacGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTgaga                                                            >  1:2924910/1‑69 (MQ=255)

    cgcgCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGa                                                        >  1:1841141/1‑69 (MQ=255)

    cgcgCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGa                                                        >  1:824578/1‑69 (MQ=255)

    cgcgCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGa                                                        >  1:1193294/1‑69 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:1972689/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:1607020/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:3310723/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:207499/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:2873693/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:3174926/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:909131/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:610279/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:1078560/1‑68 (MQ=255)

        cATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTg                                                     >  1:2441420/1‑68 (MQ=255)

         aTCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGc                                                    >  1:2965254/1‑68 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:1318597/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:2759427/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:1062640/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:373879/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:1737978/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:613609/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:672397/68‑1 (MQ=255)

            ttGGTCGCCGGTTAAAAGACTGAAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:1103244/68‑1 (MQ=255)

             tGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:912153/67‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:698461/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:2309732/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:3299401/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:2313200/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:2686774/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:1904796/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:406223/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:1490019/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:2863078/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:830919/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:1204251/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:1184878/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:2978401/69‑1 (MQ=255)

                 cGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGt                                           <  1:1018485/69‑1 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGt                                      >  1:1167642/1‑68 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:457022/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:437425/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:3128044/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:2535137/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:2284630/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:624427/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:1908027/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:410294/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:1573859/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:343061/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:76630/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:1289079/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:3433390/1‑69 (MQ=255)

                       ttAAAAGACTGGAGTGGAATACAGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:2947403/1‑69 (MQ=255)

                        tAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTc                                     >  1:2458904/1‑68 (MQ=255)

                        tAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCt                                    >  1:3032436/1‑69 (MQ=255)

                        tAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCt                                    >  1:315657/1‑69 (MQ=255)

                        tAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCt                                    >  1:1140928/1‑69 (MQ=255)

                          aaaGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGa                                                 <  1:2020307/54‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:2151909/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:3222124/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:3187812/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:2861512/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:2563707/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:518410/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:1853967/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:1775331/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGAtgt                             >  1:910104/1‑59 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:2660375/69‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:2158632/69‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:1983698/69‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:733450/69‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:1054980/69‑1 (MQ=255)

                                         aTACCGCGCTATCACGGTACCATTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCCGAt                   <  1:380372/69‑1 (MQ=255)

                                            ccGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGAtttt                >  1:987927/1‑69 (MQ=255)

                                            ccGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGAtttt                >  1:1365063/1‑69 (MQ=255)

                                                 cTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACg            >  1:3314430/1‑68 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTTCGATGTGATTGCAGATTTTGACGGCAAg       <  1:136503/68‑1 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAg       <  1:1010586/68‑1 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAg       <  1:1938705/68‑1 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAg       <  1:3076502/68‑1 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:2377864/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:2405022/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:2592560/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:977736/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:933193/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:1479996/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:1127893/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:851694/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:1229382/1‑69 (MQ=255)

                                                      aCGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGc      >  1:1164628/1‑69 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:3403216/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:917847/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:1297478/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:2627435/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:1575893/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:1840731/69‑1 (MQ=255)

                                                          tACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTc  <  1:1982944/69‑1 (MQ=255)

ACTACGCGCATCTTGGTCGCCGGTTAAAAGACTGGAGTGGAATACCGCGCTATCACGGTACCGTTGAGACGGATTGCGGAACCGGTTATGTCTACGATGTGATTGCAGATTTTGACGGCAAGCCTTC  >  minE/2173690‑2173816

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`