BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,758	C→A	8.6%	V221V (GTC→GTA) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,758	0	C	A	8.6%	333.2 / 8.6	103	V221V (GTC→GTA) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (43/51); new base A (6/3); total (49/54)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 103 total aligned reads displayed.

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTG  >  minE/2570690‑2570827

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:2968019/1‑71 (MQ=255)

        gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                             <  1:313053/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:3169225/49‑1 (MQ=255)

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

                     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                      <  1:3686092/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:2551507/1‑58 (MQ=255)

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

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

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

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

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTG  >  minE/2570690‑2570827

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