BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,756	G→A	7.9%	V221I (GTC→ATC) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,756	0	G	A	7.9%	305.3 / 11.4	101	V221I (GTC→ATC) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (42/51); new base A (5/3); total (47/54)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 101 total aligned reads displayed.

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCC  >  minE/2570690‑2570824

gAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTTAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                  <  1:1771725/71‑1 (MQ=255)

gAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                  >  1:3201044/1‑71 (MQ=255)

gAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                  >  1:613860/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

           gATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGa                                                       <  1:2273811/71‑1 (MQ=255)

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

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

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

                     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                   <  1:2994323/49‑1 (MQ=255)

                     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                   <  1:3690524/49‑1 (MQ=255)

                     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                   <  1:3686092/49‑1 (MQ=255)

                     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                   <  1:3169225/49‑1 (MQ=255)

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

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

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

                         gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                  <  1:1871820/46‑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:1909125/1‑60 (MQ=255)

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

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

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

                             gcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTc                                     <  1:457105/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)

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

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

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

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

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

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

                                      tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                            >  1:46339/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)

                                        aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                          >  1:2765907/1‑70 (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:203839/1‑58 (MQ=255)

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

                                            aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAg                                   >  1:3827925/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)

                                                     ccTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                <  1:3372093/52‑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:1923624/71‑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)

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCC  >  minE/2570690‑2570824

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