BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,614,805	G→A	16.9%	E392E (GAG→GAA)	PP_4085 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,614,805	0	G	A	16.9%	80.8 / 7.5	37	E392E (GAG→GAA)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand): ref base G (18/12); new base A (3/3); total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAG  >  NC_002947/4614668‑4614947

gCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCGTCCTCCAGCTGGAAGAAGGCCaag                                                                                                                                      <  1:19294/148‑3 (MQ=255)

  ccGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCg                                                                                                                                   >  7:37443/1‑149 (MQ=255)

           acaTCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCg                                                                                                                             >  2:47706/1‑146 (MQ=255)

                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                                 <  5:22944/119‑1 (MQ=255)

                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                                 >  6:22944/1‑119 (MQ=255)

                                        aCATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCTTGGACGACTTGGCCTTATTCCAGCTGGAGGACGGCCAGACCGTTTTAGAAg                                                                                                                          >  3:26735/1‑120 (MQ=255)

                                        aCATCGAAAACGACGGTAACGGACTGGAACTGGGCCGCCTAAGTGCCGGTGATGCCTTCATCAGTGCAGTGGACTGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAg                                                                                                                          <  4:26735/120‑1 (MQ=255)

                                         cATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCttt                                                                                                        <  1:132433/137‑1 (MQ=255)

                                         cATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCttt                                                                                                        >  2:132433/1‑137 (MQ=255)

                                               cAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGat                                                                                                                        <  7:78519/115‑1 (MQ=255)

                                               cAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGat                                                                                                                        >  8:78519/1‑115 (MQ=255)

                                                                 ggAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAAc                                                                    >  2:210607/1‑149 (MQ=255)

                                                                  gAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGAGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTATTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAAACCGAAACg                                                                   >  4:238370/1‑149 (MQ=255)

                                                                      tGGGCCGCCTCAGTGCCGGTGATGCCTTCATGAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCa                                                               >  1:91743/1‑149 (MQ=255)

                                                                            gccTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCt                                                          >  5:293363/1‑148 (MQ=255)

                                                                                     ccGGTGATGCCTTCATCAGTGCCGTGGACGGCTGGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTCTCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGcgc                                                <  6:159849/149‑3 (MQ=255)

                                                                                      cGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCg                                               >  2:191653/1‑149 (MQ=255)

                                                                                         tGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGc                                             >  2:270046/1‑148 (MQ=255)

                                                                                           aTGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGAc                                          <  1:191653/149‑1 (MQ=255)

                                                                                               cTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGc                                       >  4:273775/1‑148 (MQ=255)

                                                                                                  catcaGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAAt                                    <  1:177293/148‑1 (MQ=255)

                                                                                                                aCGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGAc                     <  3:273775/149‑1 (MQ=255)

                                                                                                                  ggCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACct                   >  8:92121/1‑149 (MQ=255)

                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                        <  5:24794/140‑1 (MQ=255)

                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                        >  6:24794/3‑142 (MQ=255)

                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc                 >  1:224613/3‑149 (MQ=255)

                                                                                                                    cTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc                 >  6:2071/1‑149 (MQ=255)

                                                                                                                     ttGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgct                <  1:47706/149‑1 (MQ=255)

                                                                                                                        gCCTTCATCCAGCTGGACGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGt              <  4:189957/148‑1 (MQ=255)

                                                                                                                        gCCGTCCTCGAGCTGTAAGAAGGCCAGACCGTTGTCGTAGATATCTATGTAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAATCGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGcgt                         <  4:46862/137‑3 (MQ=255)

                                                                                                                        gCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTAGGCGACCGCAATCTAAATTGCCt                         >  3:46862/1‑137 (MQ=255)

                                                                                                                          cttcttCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacg           >  3:249166/1‑149 (MQ=255)

                                                                                                                              ttcCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacgacg        >  3:182364/1‑148 (MQ=255)

                                                                                                                                  aGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACaa   <  7:92121/149‑1 (MQ=255)

                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCg                                        >  5:45471/1‑111 (MQ=255)

                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCg                                        <  6:45471/111‑1 (MQ=255)

                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAg  <  8:37443/149‑1 (MQ=255)

                                                                                                                                       gaggaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAAt                              <  3:81426/117‑1 (MQ=255)

                                                                                                                                       gaggaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGAAACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAAt                              >  4:81426/1‑117 (MQ=255)

GCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAG  >  NC_002947/4614668‑4614947

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