BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	821,445	T→C	38.5%	S7G (AGC→GGC)	ldrA ←	toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	821,445	0	T	C	38.5%	30.3 / 38.7	39	S7G (AGC→GGC)	ldrA	toxic polypeptide, small
Reads supporting (aligned to +/- strand): ref base T (13/11); new base C (7/8); total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGAC  >  minE/821395‑821511

cGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGt                                                <  1:2078650/71‑1 (MQ=255)

 gCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTc                                               <  1:673433/71‑1 (MQ=34)

     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCa                                              >  1:2619757/1‑68 (MQ=35)

     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCa                                              >  1:829238/1‑68 (MQ=35)

     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:1722756/1‑70 (MQ=35)

       aGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGc                                          >  1:2178173/1‑70 (MQ=35)

        gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCa                                         >  1:3702546/1‑70 (MQ=35)

        gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCa                                         >  1:973467/1‑70 (MQ=35)

          cgACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAAc                                      <  1:2772487/70‑1 (MQ=19)

           gACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACt                                     <  1:3640022/71‑1 (MQ=21)

             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   >  1:807236/1‑71 (MQ=34)

             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTg                                    >  1:2803204/1‑70 (MQ=35)

             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   >  1:591217/1‑71 (MQ=21)

               aTCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGt                                                <  1:933498/56‑1 (MQ=33)

               aTCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                 <  1:1447074/71‑1 (MQ=21)

                  gctgctGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGa                               <  1:3317541/70‑1 (MQ=25)

                   cTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGc                             >  1:467586/1‑71 (MQ=34)

                     gCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGt                           >  1:358348/1‑71 (MQ=255)

                      cGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGt                                                >  1:3782639/1‑49 (MQ=255)

                        gTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCat                        <  1:2510512/71‑1 (MQ=34)

                         taataaTTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCata                       <  1:3015779/71‑1 (MQ=34)

                         taataaTTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCata                       <  1:1158666/71‑1 (MQ=34)

                             aatTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:3704152/1‑46 (MQ=38)

                             aatTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:1257121/1‑46 (MQ=38)

                              atTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGcgcg                                <  1:233042/57‑1 (MQ=35)

                               ttCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   <  1:1125714/53‑1 (MQ=35)

                               ttCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAgcc                  >  1:1607271/1‑70 (MQ=35)

                                   cGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGTCAAACTGCGCGAGCGTCATATAGCCGCCtt             <  1:3484011/71‑1 (MQ=21)

                                   cGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCtt             <  1:1269293/71‑1 (MQ=34)

                                        ggATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCTCGAGCGTCATATAGCCGCCTtgttgt        <  1:1617364/71‑1 (MQ=255)

                                        ggATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCttgttg         >  1:1801539/1‑70 (MQ=35)

                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTtgttgt        >  1:2965086/1‑70 (MQ=25)

                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTtgttgt        >  1:2130256/1‑70 (MQ=25)

                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTa       <  1:1069279/71‑1 (MQ=21)

                                          aTCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTaa      <  1:383146/71‑1 (MQ=21)

                                            cGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                 >  1:3932904/1‑42 (MQ=255)

                                            cGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCtt             >  1:542941/1‑62 (MQ=25)

                                             ggTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATg    <  1:135627/70‑1 (MQ=25)

                                             ggTGCCGCCAGGTCGTGCCAGAAAGACATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTCAtt    <  1:3483857/70‑2 (MQ=11)

                                               tgctgcCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGac  >  1:1841920/1‑70 (MQ=35)

CGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGAC  >  minE/821395‑821511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`