BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,761	T→G	8.4%	I222M (ATT→ATG) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,761	0	T	G	8.4%	236.2 / 14.9	95	I222M (ATT→ATG) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (42/45); new base G (5/3); total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 107 total aligned reads displayed.

AGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCG  >  minE/2570691‑2570830

aGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAtt                                                                       >  1:1030077/1‑71 (MQ=255)

   aaTCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                                    <  1:3247868/71‑1 (MQ=255)

   aaTCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                                    <  1:1478451/71‑1 (MQ=255)

    aTCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                      <  1:870369/68‑1 (MQ=255)

       gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                                >  1:2968019/1‑71 (MQ=255)

       gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                                <  1:313053/71‑1 (MQ=255)

        cAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                               <  1:2106370/71‑1 (MQ=255)

         aGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                               >  1:2891909/1‑70 (MQ=255)

               tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                        <  1:3222161/71‑1 (MQ=255)

                 aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                                      >  1:1037541/1‑71 (MQ=255)

                  gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                         >  1:2221381/1‑67 (MQ=255)

                     tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTAc                                                  >  1:2667099/1‑71 (MQ=255)

                      aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                 >  1:367369/1‑71 (MQ=255)

                       agagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAtt                                                                       <  1:3528114/48‑1 (MQ=255)

                        gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                  >  1:1455185/1‑52 (MQ=255)

                        gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCg                                                              <  1:3760123/56‑1 (MQ=255)

                        gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                                >  1:2672745/1‑70 (MQ=255)

                         agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                              >  1:95976/1‑71 (MQ=255)

                          gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                      >  1:955657/1‑46 (MQ=255)

                          gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                        >  1:1913732/1‑60 (MQ=255)

                          gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                        >  1:1909125/1‑60 (MQ=255)

                           cgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                            <  1:2765956/71‑1 (MQ=255)

                           cgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                            <  1:643878/71‑1 (MQ=255)

                             cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTc                                           <  1:1830841/70‑1 (MQ=255)

                             cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCa                                          <  1:2194854/71‑1 (MQ=255)

                              aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                                <  1:783101/64‑1 (MQ=255)

                              aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAg                                         <  1:50951/71‑1 (MQ=255)

                               aCTTACTGAAGCTAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                        >  1:3601015/1‑71 (MQ=255)

                               aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                      >  1:3701991/1‑41 (MQ=255)

                               aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                              >  1:896908/1‑65 (MQ=255)

                               aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCa                                          <  1:3559689/69‑1 (MQ=255)

                                 ttACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                      >  1:1399077/1‑71 (MQ=255)

                                  tACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTAc                                     >  1:1029350/1‑71 (MQ=255)

                                  tACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTAc                                     >  1:1561352/1‑71 (MQ=255)

                                   gcggAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACt                                    <  1:2379157/68‑1 (MQ=255)

                                    cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                    <  1:1591148/70‑1 (MQ=255)

                                    cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                    <  1:2908671/70‑1 (MQ=255)

                                    cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                    >  1:236619/1‑70 (MQ=255)

                                     tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:46339/1‑71 (MQ=255)

                                     tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:2540261/1‑71 (MQ=255)

                                     ggAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGa                                  >  1:2949883/2‑71 (MQ=255)

                                      gAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                    <  1:1056164/68‑1 (MQ=255)

                                      gAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCCAGCAATGGCt                                                    <  1:2213941/52‑1 (MQ=255)

                                       aaGCAAAATATAACCTTTAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                                      >  1:3524585/1‑49 (MQ=255)

                                       aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGAt                                 >  1:2247349/1‑70 (MQ=255)

                                       aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                                >  1:1762874/1‑71 (MQ=255)

                                        aGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATcc                               >  1:1456328/1‑71 (MQ=255)

                                        aGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTat                                                  >  1:1921474/1‑51 (MQ=255)

                                         gaaaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTc                                           <  1:1554036/56‑1 (MQ=255)

                                         gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                         >  1:1233979/1‑44 (MQ=255)

                                         gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCg                              >  1:1420436/1‑71 (MQ=255)

                                         gCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGc                                                           <  1:1056959/42‑1 (MQ=16)

                                          cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                             >  1:1267847/1‑71 (MQ=255)

                                           aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                    <  1:76420/47‑1 (MQ=255)

                                           aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAg                                         >  1:2551507/1‑58 (MQ=255)

                                           aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAg                                         >  1:3827925/1‑58 (MQ=255)

                                           aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAg                                         >  1:203839/1‑58 (MQ=255)

                                           aaaaTATAACCTTCAGCAGACTGTCATTCGCGCGCCGAGCa                                                          <  1:2016311/41‑1 (MQ=255)

                                            aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                           <  1:1325570/39‑1 (MQ=255)

                                            aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCa                                          >  1:193785/1‑56 (MQ=255)

                                            aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                           >  1:2760089/1‑71 (MQ=255)

                                            aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                           >  1:687261/1‑71 (MQ=255)

                                             aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCg                                                              <  1:2758292/35‑1 (MQ=255)

                                             aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                           <  1:3053609/38‑1 (MQ=255)

                                                 taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCa                                          <  1:2179532/51‑1 (MQ=255)

                                                 taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGt                       <  1:228450/70‑1 (MQ=255)

                                                 taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                      <  1:2530590/71‑1 (MQ=255)

                                                  atccTTGAGGAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTc                                           <  1:26174/47‑1 (MQ=255)

                                                  aaCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                    <  1:1226409/40‑1 (MQ=255)

                                                   aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGAt                                 >  1:1236603/1‑58 (MQ=255)

                                                   aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                    <  1:172741/71‑1 (MQ=255)

                                                   aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                    <  1:597022/71‑1 (MQ=255)

                                                     cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                    <  1:1013344/37‑1 (MQ=255)

                                                     cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACAt                   >  1:3429211/1‑70 (MQ=255)

                                                     cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa                  <  1:605910/71‑1 (MQ=255)

                                                        gAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg                <  1:219833/70‑1 (MQ=255)

                                                         aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGTTACATACGCa              <  1:124614/71‑1 (MQ=255)

                                                         aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATAc                 <  1:2264375/68‑1 (MQ=255)

                                                          gCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAg             >  1:3662023/1‑71 (MQ=255)

                                                           cAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAg             <  1:2766676/70‑1 (MQ=255)

                                                             gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTAc                                     <  1:2252279/44‑1 (MQ=255)

                                                             gACTGTAATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                      >  1:1486909/1‑59 (MQ=255)

                                                              aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                      <  1:2367045/42‑1 (MQ=255)

                                                              aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                         >  1:2567420/1‑55 (MQ=255)

                                                               cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGAt                                 >  1:3114791/1‑46 (MQ=255)

                                                               cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGcc        <  1:2710821/71‑1 (MQ=255)

                                                                  tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGt                       >  1:3534992/1‑53 (MQ=255)

                                                                  tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTg     >  1:2686476/1‑71 (MQ=255)

                                                                    aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc   >  1:3900848/1‑71 (MQ=255)

                                                                    aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGc    <  1:292297/70‑1 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  >  1:3757486/1‑71 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  <  1:1976035/71‑1 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  <  1:2242318/71‑1 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  >  1:1142127/1‑71 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  >  1:735532/1‑71 (MQ=255)

                                                                     ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg  <  1:3869112/71‑1 (MQ=255)

                                                                      tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:523345/1‑38 (MQ=255)

                                                                      tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:3128848/1‑38 (MQ=255)

                                                                      tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:124238/1‑38 (MQ=255)

                                                                      tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                  >  1:3298004/1‑38 (MQ=255)

AGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCG  >  minE/2570691‑2570830

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