BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,234	G→A	G129G (GGG→GGA)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,234	0	G	A	98.8%	268.2 / NA	81	G129G (GGG→GGA)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (34/46); minor base T (1/0); total (35/46)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCC  >  NC_000913/3765184‑3765323

cGCTATTATGTCCCGCTGCGGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGAATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:1584199/93‑1 (MQ=11)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:354005/93‑1 (MQ=9)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:604879/93‑1 (MQ=9)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:4593610/93‑1 (MQ=9)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                                >  1:2715537/1‑93 (MQ=12)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                               <  1:3524977/93‑1 (MQ=255)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                              <  2:590258/93‑1 (MQ=255)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                              <  2:590257/93‑1 (MQ=255)

   tattatGACCCGCTGCTGGGGCGCTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                              <  2:3390030/93‑1 (MQ=255)

     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                                            >  1:643003/1‑93 (MQ=255)

     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                                            >  1:643004/1‑93 (MQ=255)

     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                                            >  1:724484/1‑93 (MQ=255)

      tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa                                           >  1:2587753/1‑93 (MQ=255)

       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa                                          <  2:1950190/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                         <  2:2873528/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                         <  1:1446562/93‑1 (MQ=255)

         gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                                        <  2:388373/93‑1 (MQ=255)

         gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                                        <  2:2963096/93‑1 (MQ=255)

           cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                                      >  2:2245881/1‑93 (MQ=255)

           cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                                      >  2:749588/1‑93 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                                     <  1:3303703/93‑1 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                                     >  1:3615237/1‑93 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                                     <  2:2686297/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  2:4741205/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    <  1:4508384/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    <  2:3609880/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  1:3647122/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  1:3647123/1‑93 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                                  <  2:215508/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 <  1:4200790/93‑1 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                <  2:312131/90‑1 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                >  1:4540419/4‑93 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                >  2:3362612/4‑93 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                <  2:1419237/90‑1 (MQ=255)

                  ctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCg                               >  1:4690296/3‑93 (MQ=255)

                    ggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                             >  2:3703513/1‑93 (MQ=255)

                    ggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                             <  2:1866310/93‑1 (MQ=255)

                    ggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                             >  2:2937529/1‑93 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:2599631/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:1668493/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:1489552/93‑1 (MQ=255)

                      ggaggTATATTGCTCGGGAGCCGAGTGGACTGAAGGGGGTATGGAATTTTTATCAGTATCCGTTGAATCCTGTCAGAAATGTAGATCCGCAAg                           <  1:1291273/90‑1 (MQ=17)

                        cGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGt                         <  2:276881/93‑1 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                        <  1:1807411/93‑1 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGAAAGGtt                        <  1:3158823/93‑1 (MQ=17)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       >  1:3637928/1‑93 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       >  2:4402923/1‑93 (MQ=255)

                           tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg                      <  1:3492968/93‑1 (MQ=255)

                           tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg                      >  1:1841070/1‑93 (MQ=255)

                           tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg                      >  2:456919/1‑93 (MQ=255)

                           gatatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg                      <  1:773181/92‑1 (MQ=255)

                            atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg                     >  1:737646/1‑93 (MQ=255)

                             tatCACTCAGGATCCGCTTGGACTGATGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt                    <  1:4214175/93‑1 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 <  2:4321429/80‑1 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 >  1:4321429/1‑80 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt                    <  1:2572985/93‑1 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt                    >  2:3147493/1‑93 (MQ=255)

                              atCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGtt                   <  1:2700550/93‑1 (MQ=255)

                              atCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGtt                   <  1:2156067/93‑1 (MQ=255)

                                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGat                >  1:3775734/1‑93 (MQ=255)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata               >  1:3901879/1‑93 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGCTGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAAATCCGCAAGGTTTGGTTGataa              >  2:1270615/1‑92 (MQ=21)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              <  2:2721682/93‑1 (MQ=255)

                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTATTATCATTATGCGTTGAATCCGGTCATAAATGTAGATCAGCAAGGTTTGGTTGATATaa            <  1:1151722/93‑1 (MQ=25)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa           <  1:3630348/93‑1 (MQ=255)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa           <  2:1571589/93‑1 (MQ=255)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa           >  2:4115990/1‑93 (MQ=255)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa           >  1:3102571/1‑93 (MQ=255)

                                      ggATACGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATACATGTAGATCCGCAAGGTTTGGTTGATATaaa           <  2:4432716/93‑1 (MQ=255)

                                       gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAt          <  2:2814949/93‑1 (MQ=255)

                                       gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAt          >  1:1594532/1‑93 (MQ=255)

                                       gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAt          <  2:3227052/93‑1 (MQ=255)

                                        aTCCGATTGGTCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCATGGCTTGGTTGATATAAAtt         >  1:4197238/1‑93 (MQ=255)

                                        aTCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAtt         <  1:3028981/93‑1 (MQ=255)

                                         tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCCTAAATGTAGATCCGCAAGGTTTGGTTGATATAAAttt        >  2:2519715/1‑93 (MQ=255)

                                         tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAttt        >  2:2532928/1‑93 (MQ=255)

                                         tCCGATTGGACAGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCTGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAttt        <  2:1970351/93‑1 (MQ=255)

                                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTta       <  2:1721674/93‑1 (MQ=255)

                                             aTTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAc    <  1:773138/93‑1 (MQ=255)

                                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAcc   >  2:163006/1‑93 (MQ=255)

                                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGGCATAAATGTAGAGCCGCAAGGTTTGGTTTATATAAATTTATAcc   >  2:163007/1‑93 (MQ=25)

                                               tGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATcccc  <  2:4078554/93‑1 (MQ=255)

                                               tGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAccc  <  1:1819111/93‑1 (MQ=255)

                                               tGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAccc  <  1:3593372/93‑1 (MQ=255)

                                               tGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAccc  >  1:3037785/1‑93 (MQ=255)

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCC  >  NC_000913/3765184‑3765323

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