BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,968,634	G→A	G305G (GGC→GGT)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,968,634	0	G	A	98.2%	193.2 / NA	56	G305G (GGC→GGT)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (25/29); minor base C (0/1); total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGA  >  REL606/2968565‑2968705

aaGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCttt                                                                   >  2:1233849/1‑76 (MQ=255)

 aGGCGTGCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCtttt                                                                  <  2:1644117/76‑1 (MQ=255)

  ggCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCttttt                                                                 <  2:1010868/76‑1 (MQ=255)

     gTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGTGCATCGCATCTGGTGCTTACTCACCCTTTTTcgc                                                              >  2:572326/1‑76 (MQ=255)

       tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTcgccg                                                            >  1:546175/1‑76 (MQ=255)

        cACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTcgccgc                                                           >  2:1702854/1‑76 (MQ=255)

         aCGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCgccgcc                                                          <  2:822713/76‑1 (MQ=255)

          cgccgcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCt                                                         <  1:1171858/76‑1 (MQ=255)

          cgccgcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCt                                                         <  2:1640122/76‑1 (MQ=255)

           gccgcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTg                                                        <  2:408589/76‑1 (MQ=255)

                aTCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGTTTGAatat                                                   >  2:157225/1‑76 (MQ=255)

                aTCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGTCTGAatat                                                   <  2:1485908/76‑1 (MQ=255)

                aTCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAatat                                                   <  2:41356/76‑1 (MQ=255)

                aTCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAatat                                                   >  1:1756162/1‑76 (MQ=255)

                  ccGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATAc                                                 <  1:1303535/76‑1 (MQ=255)

                  ccGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATAc                                                 >  1:496894/1‑76 (MQ=255)

                       aTTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCAt                                            <  1:1742292/76‑1 (MQ=255)

                       aTTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCAt                                            >  2:129011/1‑76 (MQ=255)

                        ttAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCAtt                                           <  1:1747140/76‑1 (MQ=255)

                         tAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCAttt                                          <  2:1562220/76‑1 (MQ=255)

                          aaTGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTc                                         <  1:1392483/76‑1 (MQ=255)

                          aaTGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTc                                         >  2:573458/1‑76 (MQ=255)

                           aTGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTcc                                        <  2:530083/76‑1 (MQ=255)

                           aTGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTcc                                        >  1:488186/1‑76 (MQ=255)

                            tGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTcca                                       <  2:1603461/76‑1 (MQ=255)

                            tGGCACATTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTcca                                       <  1:1341271/76‑1 (MQ=255)

                             ggCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccag                                      <  2:1099846/76‑1 (MQ=255)

                                 cGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgcc                                  <  1:67747/76‑1 (MQ=255)

                                   ttttACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgccag                                >  2:726917/1‑76 (MQ=255)

                                    tttACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgccagc                               <  2:238770/76‑1 (MQ=255)

                                     ttACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgccagcg                              >  1:1363405/1‑76 (MQ=255)

                                     ttACCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgccagcg                              >  1:23858/1‑76 (MQ=255)

                                       aCCCGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGtt                            >  1:1400440/1‑76 (MQ=255)

                                          cGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGcc                         <  1:759983/76‑1 (MQ=255)

                                          cGTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGcc                         <  1:1478223/76‑1 (MQ=255)

                                           gTGCGCATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCg                        >  1:1470165/1‑76 (MQ=255)

                                             gcgcATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCt                      <  2:1514513/76‑1 (MQ=255)

                                             gcgcATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCt                      <  1:1514289/76‑1 (MQ=255)

                                             gcgcATCGCATCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCt                      <  2:1164919/76‑1 (MQ=255)

                                              cgcatcgcatCTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTg                     >  2:1559393/1‑76 (MQ=255)

                                                catcgcatcTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGcc                   >  2:70380/1‑76 (MQ=255)

                                                  tcgcatcTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCag                 <  2:1576137/76‑1 (MQ=255)

                                                   cgcatcTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCaga                <  1:510073/76‑1 (MQ=255)

                                                    gcgtcTGGTGCTTACTCCCCCTGTGTCGCCGCCTGAATCTACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCagag               <  2:145901/73‑1 (MQ=255)

                                                     catcTGGTGCTTACTCTCCGTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGc              <  2:82250/76‑1 (MQ=255)

                                                      atcTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCa             >  2:614983/1‑76 (MQ=255)

                                                      atcTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCa             >  1:1066815/1‑76 (MQ=255)

                                                        cTGGTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCCGt           >  2:359068/1‑76 (MQ=255)

                                                          ggTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGa         <  1:1265994/76‑1 (MQ=255)

                                                           gTGCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGAt        >  2:1667180/1‑76 (MQ=255)

                                                             gCTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGAttt      <  1:1086427/76‑1 (MQ=255)

                                                              cTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc     <  1:1163450/76‑1 (MQ=255)

                                                              cTTACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc     >  1:1527600/1‑76 (MQ=255)

                                                                tACTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCg   >  2:10169/1‑76 (MQ=255)

                                                                 ncTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGa  >  1:558539/2‑76 (MQ=255)

                                                                 aCTCACCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGa  >  1:1187996/1‑76 (MQ=255)

AAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGA  >  REL606/2968565‑2968705

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`