BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,186,043	A→G	G159G (GGA→GGG)	yqiJ →	DUF1449 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,186,043	0	A	G	91.9%	106.1 / 2.6	37	G159G (GGA→GGG)	yqiJ	DUF1449 family protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (19/15); total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GTGCATTATACCGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGA  >  NZ_CP009273/3185909‑3186168

gTGCATTATACCGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAAc                                                                                                                           >  1:490070/1‑139 (MQ=255)

 tGCATTATACCGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACt                                                                                                                          <  2:271560/139‑1 (MQ=255)

     ttATACCGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCAcc                                                                                                                      >  2:506764/1‑139 (MQ=255)

           cGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGAACACTCACCGATCaa                                                                                                                >  2:494937/1‑139 (MQ=255)

             gTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAAtt                                                                                                              >  2:92895/1‑139 (MQ=255)

             gTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAAtt                                                                                                              >  2:542772/1‑139 (MQ=255)

             gTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAAtt                                                                                                              >  2:416453/1‑139 (MQ=255)

                  aTAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGc                                                                                                          >  1:609888/1‑138 (MQ=255)

                         cgcCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  2:609888/139‑1 (MQ=255)

                         cgcCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  1:407140/139‑1 (MQ=255)

                         cgcCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  1:328589/139‑1 (MQ=255)

                         cgcCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGAAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  2:402098/139‑1 (MQ=255)

                                gATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACattatt                                                                                             <  1:92895/137‑1 (MQ=255)

                                gATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACattatt                                                                                             <  2:49025/137‑1 (MQ=255)

                                          gACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  2:692726/122‑1 (MQ=255)

                                          gACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  >  1:692726/1‑122 (MQ=255)

                                          gACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGc                                                                                 >  1:101029/1‑139 (MQ=255)

                                             cATAGTTCAGCTATTACAGAAGAAGAATAAATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCt                                                                               >  2:679802/1‑138 (MQ=255)

                                                      gCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGa                                                                      >  1:340735/1‑138 (MQ=255)

                                                          ttACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  >  1:103165/1‑106 (MQ=255)

                                                          ttACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACa                                                                                                  <  2:103165/106‑1 (MQ=255)

                                                                        tatatTGGCAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAg                                                      >  2:208059/1‑136 (MQ=255)

                                                                            ttGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACaa                                               <  1:506764/139‑1 (MQ=255)

                                                                             tGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACa                                                >  1:741724/1‑137 (MQ=255)

                                                                             tGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACa                                                <  2:741724/137‑1 (MQ=255)

                                                                                          ttAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTAttt                                 >  1:476656/1‑139 (MQ=255)

                                                                                                  cGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGc                               <  2:523037/133‑1 (MQ=255)

                                                                                                  cGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGc                               >  1:523037/1‑133 (MQ=255)

                                                                                                    gTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACa                                                <  1:700554/114‑1 (MQ=255)

                                                                                                    gTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACa                                                >  2:700554/1‑114 (MQ=255)

                                                                                                       atcaAGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTg                     >  2:549025/1‑138 (MQ=255)

                                                                                                          aaGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAAc                 <  2:490070/139‑1 (MQ=255)

                                                                                                           aGCCACTTCAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAAcg                <  2:325230/139‑1 (MQ=255)

                                                                                                                   cAGGTAATCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATAttt        >  2:58925/1‑139 (MQ=255)

                                                                                                                      gTAATCCCTGCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGc     >  1:401098/1‑139 (MQ=255)

                                                                                                                         aTCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCgg   <  2:340735/138‑1 (MQ=255)

                                                                                                                          tCCCTGCGAAGGGAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGa  >  1:726443/1‑138 (MQ=255)

GTGCATTATACCGGTAAAATAGTTGCGCCCTGGATTCCACGCGACCATAGTTCAGCTATTACAGAAGAAGAATATATTGGTAGCATGGCGTTAATTACCGGTCATCAAGCCACTTCAGGTAATCCCTGCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGA  >  NZ_CP009273/3185909‑3186168

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