BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,764	C→T	8.2%	R223R (CGC→CGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,764	0	C	T	8.2%	289.6 / 11.4	97	R223R (CGC→CGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (45/44); new base T (5/3); total (50/47)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 106 total aligned reads displayed.

AATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCT  >  minE/2570694‑2570832

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

aaTCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                                      <  1:3247868/71‑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)

       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)

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

                   aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                   >  1:367369/1‑71 (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)

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

                       gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                          >  1:1909125/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)

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

                                 cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                      <  1:1591148/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)

                                    aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                                  >  1:2765907/1‑70 (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)

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

                                                aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                      <  1:172741/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)

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

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

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

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

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

                                                                    cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  >  1:3179459/1‑71 (MQ=255)

                                                                    cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  >  1:3438704/1‑71 (MQ=255)

                                                                     gcgcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  <  1:664843/70‑1 (MQ=255)

                                                                     gcgcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  <  1:481804/70‑1 (MQ=255)

AATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCT  >  minE/2570694‑2570832

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