BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,759	A→G	8.1%	I222V (ATT→GTT) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,759	0	A	G	8.1%	280.8 / 12.0	99	I222V (ATT→GTT) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (44/47); new base G (5/3); total (49/50)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 105 total aligned reads displayed.

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCC  >  minE/2570690‑2570829

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

    aaTCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                                   <  1:1478451/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)

           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)

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

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

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

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

                                                              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)

GAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCC  >  minE/2570690‑2570829

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