BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	311,500	G→A	G176D (GGC→GAC)	yahG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	311,500	0	G	A	98.8%	1034.3 / ‑0.7	259	G176D (GGC→GAC)	yahG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (161/95); total (164/95)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 267 total aligned reads displayed.

GACTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGGCAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCGT  >  REL606/311427‑311572

gACTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACa                                                                        >  1:1803045/1‑76 (MQ=255)

gACTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACa                                                                        >  1:2375122/1‑76 (MQ=255)

  cTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAAc                                                                      >  2:5623913/1‑76 (MQ=255)

  cTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAAc                                                                      >  1:2904412/1‑76 (MQ=255)

    gCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGacaacac                                                                    >  2:6006968/1‑75 (MQ=255)

      gTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCg                                                                  <  2:9199606/76‑1 (MQ=255)

      gTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCg                                                                  <  1:5197752/76‑1 (MQ=255)

      gTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCg                                                                  <  1:3170589/76‑1 (MQ=255)

        gggATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCt                                                                <  1:5625718/76‑1 (MQ=255)

           aTCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTtat                                                             <  1:7968775/76‑1 (MQ=255)

           aTCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTtat                                                             <  1:5087886/76‑1 (MQ=255)

              gATGGCGGGCGGTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAAAAAAACCTACGACCACCTCGCTtataac                                                          >  2:7793408/1‑74 (MQ=255)

                tGGCGGGCGTTAACTCGGCCTAGATGTTTATGCACATCGTGAAAAACAAAACCTACGATAAAACCGCTTTTACCaa                                                        >  2:9363515/1‑76 (MQ=255)

                    gggcgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATg                                                    >  1:7438970/1‑76 (MQ=255)

                    gggcgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATg                                                    >  1:4146653/1‑76 (MQ=255)

                    gggcgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATg                                                    >  1:859446/1‑76 (MQ=255)

                    gggcgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATg                                                    >  1:9392264/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagt                                                 >  1:7232684/1‑75 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:8185176/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:8421689/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:6700951/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:4388634/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:8497826/1‑76 (MQ=255)

                       cgTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagc                                                 >  1:619298/1‑76 (MQ=255)

                        gTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGGCAACATCGCTTATACCAATATgagcg                                                >  2:7995090/1‑76 (MQ=255)

                        gTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagcg                                                >  2:1143311/1‑76 (MQ=255)

                          tACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagcgag                                              <  2:6452003/76‑1 (MQ=255)

                          tACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagcgag                                              >  1:4947902/1‑76 (MQ=255)

                          tACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATgagcgag                                              >  2:3281578/1‑76 (MQ=255)

                             cTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAg                                           <  2:4037453/76‑1 (MQ=255)

                             cTCGGCCGCGATGTTTATGCACATCTTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAg                                           <  2:8667196/76‑1 (MQ=255)

                              tCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGa                                          <  1:3744985/76‑1 (MQ=255)

                              tCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGa                                          <  1:4569437/76‑1 (MQ=255)

                              tCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGa                                          <  1:7689399/76‑1 (MQ=255)

                               cGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGAt                                         <  2:1244567/76‑1 (MQ=255)

                               cGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGAt                                         <  2:5081800/76‑1 (MQ=255)

                               cGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGAt                                         <  2:5348995/76‑1 (MQ=255)

                                ggCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATg                                        <  2:6286191/76‑1 (MQ=255)

                                ggCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATg                                        >  2:4207803/1‑76 (MQ=255)

                                    tCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGa                                    <  1:8066983/76‑1 (MQ=255)

                                    tCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGa                                    >  2:1462533/1‑76 (MQ=255)

                                    tCGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGa                                    >  2:3467938/1‑76 (MQ=255)

                                    tCGATGTTTATGCAAAGCGTGAAAAACAAACCCTACGACAACATCGCTTATCCCAATATGAGCGAGCAGATGGCGa                                    <  1:6151473/76‑1 (MQ=255)

                                     cGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGaa                                   >  1:9467418/1‑76 (MQ=255)

                                     cGATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGaa                                   >  1:2190757/1‑76 (MQ=255)

                                      gATGTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAg                                  >  1:9326056/1‑76 (MQ=255)

                                         gTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAtt                               >  2:8446469/1‑76 (MQ=255)

                                         gTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAtt                               <  1:8832573/76‑1 (MQ=255)

                                         gTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAtt                               >  1:4586796/1‑76 (MQ=255)

                                         gTTTATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAtt                               <  1:9534427/76‑1 (MQ=255)

                                           ttATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAtttt                             <  1:7556862/76‑1 (MQ=255)

                                            tATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAGGATTTTg                            >  1:7814356/1‑76 (MQ=255)

                                            tATGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTg                            <  2:7684930/76‑1 (MQ=255)

                                             aTGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGc                           >  1:4104538/1‑76 (MQ=255)

                                              tGCACATCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCg                          >  1:4845465/1‑76 (MQ=255)

                                                  caTCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATg                      >  1:6827815/1‑76 (MQ=255)

                                                  caTCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATg                      >  1:5451270/1‑76 (MQ=255)

                                                    tCGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATggg                    >  1:5836359/1‑76 (MQ=255)

                                                     cGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGt                   <  1:5397052/76‑1 (MQ=255)

                                                     cGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGt                   <  2:8257492/76‑1 (MQ=255)

                                                     cGTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGt                   <  1:3281578/76‑1 (MQ=255)

                                                      gTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAttt                              >  2:4577213/1‑64 (MQ=255)

                                                      gTGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGAttt                              >  2:5746144/1‑64 (MQ=255)

                                                       tGAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGc                 <  1:4441527/76‑1 (MQ=255)

                                                        gAAAAACAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCt                >  1:1021885/1‑76 (MQ=255)

                                                         aaaaaCAAAACCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTa               >  2:3152467/1‑76 (MQ=255)

                                                          aaaacaaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTaa              <  1:3781193/76‑1 (MQ=255)

                                                          aaaacaaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTaa              >  2:8965056/1‑76 (MQ=255)

                                                             acaaaccCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGa           <  1:4006215/76‑1 (MQ=255)

                                                              caaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGAc          <  1:358780/76‑1 (MQ=255)

                                                              caaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGAc          >  1:1888727/1‑76 (MQ=255)

                                                               aaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGAcc         <  1:3315288/76‑1 (MQ=255)

                                                               aaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGAcc         >  2:7076265/1‑76 (MQ=255)

                                                               aaaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGAcc         <  2:4084733/76‑1 (MQ=255)

                                                                aaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCa        >  2:8760645/1‑76 (MQ=255)

                                                                aaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCa        >  2:1684912/1‑76 (MQ=255)

                                                                aaacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCa        >  2:2593552/1‑76 (MQ=255)

                                                                 aacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCag       <  1:7051275/76‑1 (MQ=255)

                                                                 aacCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCag       >  2:1609502/1‑76 (MQ=255)

                                                                  acCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCaga      <  1:1206094/76‑1 (MQ=255)

                                                                  acCTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCaga      >  1:6261344/1‑76 (MQ=255)

                                                                   ccTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCagag     <  1:6586741/76‑1 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:9177121/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:8299255/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:6690519/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:2266953/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:6547940/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:260203/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:2872932/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:9123053/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:4501568/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:602929/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:5758710/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:5292954/1‑76 (MQ=255)

                                                                    cTACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGc    >  1:4295584/1‑76 (MQ=255)

                                                                     tACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCg   >  2:2474247/1‑76 (MQ=255)

                                                                     tACGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCg   <  1:1640488/76‑1 (MQ=255)

                                                                      aCGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTTCTAACGACCCAAGCgg  >  2:2924173/1‑75 (MQ=255)

                                                                      aCGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCGt  >  2:1193254/1‑76 (MQ=255)

                                                                      aCGACAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCGt  >  2:163812/1‑76 (MQ=255)

GACTGCGTGGGATCGATGGCGGGCGTTACCTCGGCCTCGATGTTTATGCACATCGTGAAAAACAAAACCTACGGCAACATCGCTTATACCAATATGAGCGAGCAGATGGCGAAGATTTTGCGTATGGGTGCTAACGACCAGAGCGT  >  REL606/311427‑311572

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