BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,124	C→T	7.8%	I343I (ATC→ATT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,124	0	C	T	7.8%	350.8 / 13.3	115	I343I (ATC→ATT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (62/44); new base T (6/3); total (68/47)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 119 total aligned reads displayed.

GCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAA  >  minE/2571054‑2571191

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:778516/1‑71 (MQ=255)

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:1084205/1‑71 (MQ=255)

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     <  1:484090/71‑1 (MQ=255)

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:920613/1‑71 (MQ=255)

 cACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    <  1:130901/71‑1 (MQ=255)

     gacggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGc                                                                >  1:1474017/1‑71 (MQ=255)

      acggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGcc                                                               >  1:2809320/1‑71 (MQ=255)

       cggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGccc                                                              <  1:2890056/71‑1 (MQ=255)

            tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTg                                                         >  1:3348802/1‑71 (MQ=255)

            tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTg                                                         >  1:1394229/1‑71 (MQ=255)

            ggCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTg                                                         <  1:971848/70‑1 (MQ=255)

               tGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCg                                                      >  1:1924455/1‑71 (MQ=255)

                gggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     >  1:1207196/1‑71 (MQ=255)

                 ggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:3121236/1‑57 (MQ=255)

                 ggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     <  1:1129813/70‑1 (MQ=255)

                 ggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     <  1:3638598/70‑1 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   <  1:3886982/71‑1 (MQ=255)

                   aaCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:1040565/1‑70 (MQ=255)

                     ccATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    <  1:215910/51‑1 (MQ=255)

                      cATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                               <  1:1072581/71‑1 (MQ=255)

                       aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:831212/1‑51 (MQ=255)

                       aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              <  1:2012768/71‑1 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  <  1:2137484/50‑1 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:2924826/1‑65 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                             >  1:2782638/1‑71 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              <  1:2857808/70‑1 (MQ=255)

                        ttGAACTGGACCATAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:2257528/1‑65 (MQ=255)

                         tGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGc                                                       >  1:657469/1‑60 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:2853522/1‑48 (MQ=255)

                           aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     >  1:3641952/1‑60 (MQ=255)

                           aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCAAGGTGGCgg                                                     >  1:2366098/1‑60 (MQ=255)

                            aCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAcc                                         <  1:1789316/71‑1 (MQ=255)

                            aCTGGACACTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                          >  1:801584/1‑70 (MQ=255)

                             cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCa                                        >  1:2231903/1‑71 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:1893495/1‑71 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                      >  1:2102634/1‑70 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAt                                       <  1:181142/69‑1 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTAACCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGaa                                          >  1:1334145/1‑65 (MQ=255)

                                  cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:2810103/1‑69 (MQ=255)

                                  cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:2571257/1‑69 (MQ=255)

                                   ccTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    >  1:1573982/1‑37 (MQ=255)

                                   ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGttt                                                  <  1:2134404/55‑1 (MQ=255)

                                   ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                  >  1:1721624/1‑71 (MQ=255)

                                   ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                  >  1:3908733/1‑71 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    >  1:2301272/1‑35 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCa                                                             <  1:149089/42‑1 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:1274267/1‑57 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCg                                            <  1:593424/59‑1 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                 >  1:2467528/1‑70 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                 >  1:3405967/1‑70 (MQ=255)

                                      aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               >  1:3696805/1‑71 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                           >  1:3162418/1‑42 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                           >  1:1069144/1‑42 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                              >  1:2296564/1‑71 (MQ=255)

                                        cGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCa                                   >  1:828070/1‑65 (MQ=255)

                                         gatgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGt                            >  1:1919951/1‑71 (MQ=255)

                                           tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                          <  1:1883885/71‑1 (MQ=255)

                                               atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:3848879/62‑1 (MQ=255)

                                               atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:2926085/62‑1 (MQ=255)

                                               atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                        <  1:1803250/69‑1 (MQ=255)

                                                tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:504025/61‑1 (MQ=255)

                                                tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      <  1:1753831/70‑1 (MQ=255)

                                                 cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:3213189/1‑48 (MQ=255)

                                                 cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATg                             >  1:3206653/1‑62 (MQ=255)

                                                 cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    >  1:3917161/1‑71 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    <  1:445920/70‑1 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    <  1:2318931/70‑1 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    <  1:1608881/70‑1 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                   <  1:1329041/71‑1 (MQ=255)

                                                   aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:3485451/58‑1 (MQ=255)

                                                    tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:1718118/1‑42 (MQ=255)

                                                    tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:194779/1‑45 (MQ=255)

                                                    tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:353665/1‑45 (MQ=255)

                                                       tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:1066394/2‑36 (MQ=25)

                                                       tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:122854/2‑36 (MQ=25)

                                                       tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:1826468/2‑36 (MQ=25)

                                                       tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgc               >  1:2862143/2‑69 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:3900741/1‑39 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      <  1:1740248/63‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      <  1:2898723/63‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa              <  1:707775/71‑1 (MQ=255)

                                                         tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa              >  1:3348979/1‑69 (MQ=255)

                                                         tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca              <  1:3022781/69‑3 (MQ=255)

                                                          aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa             <  1:2904217/69‑1 (MQ=255)

                                                          aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg           >  1:21810/1‑71 (MQ=255)

                                                           cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGTTGAAAg           <  1:508836/70‑1 (MQ=255)

                                                           cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg           <  1:3333110/70‑1 (MQ=255)

                                                             cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc        >  1:2823678/1‑71 (MQ=255)

                                                             agACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc        >  1:3886495/2‑71 (MQ=255)

                                                              gACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc        >  1:1792418/1‑70 (MQ=255)

                                                              gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct       >  1:1197432/1‑71 (MQ=255)

                                                              gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct       >  1:3001520/1‑71 (MQ=255)

                                                               aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                              >  1:2463630/1‑47 (MQ=255)

                                                                cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat     <  1:632030/71‑3 (MQ=255)

                                                                cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat     <  1:2438973/71‑3 (MQ=255)

                                                                cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctgc     >  1:652245/1‑71 (MQ=255)

                                                                  gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct    <  1:203597/70‑1 (MQ=255)

                                                                  gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGCtctgct    <  1:30981/70‑1 (MQ=255)

                                                                   cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct    >  1:2786236/1‑69 (MQ=255)

                                                                    aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa              <  1:3723872/58‑1 (MQ=255)

GCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAA  >  minE/2571054‑2571191

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