BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,978,150	G→A	15.4%	noncoding (2651/2928 nt)	EO946_RS15430 ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP034943	2,978,150	0	G	A	15.4%	1.2 / 10.5	39	noncoding (2651/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (7/26); new base A (2/4); total (9/30)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

AGATGCTTTCAGCACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAA  >  NZ_CP034943/2978016‑2978280

aaatccTTTCAGCACAAATCCCGTCCGCACAAAGCTACACAGAGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATACAGGTCCTCACGTACTAAGTACACCTCCTCTCAAATATCCGGCGCCCGCGACg                                                                                                                               >  2:120340/6‑140 (MQ=1)

   tggtttCAGAACTTATCCCGTCCGCACATAGCAACCCAGCGTGGCCGTTGGAAGAACACGGGGTACACCCGGGGTGCTTCAATCCGGGGCCTATCTAACTACGCACAGCTCCTTAAAATTTTCCTGCGCCGGCGACgttt                                                                                                                            <  2:275598/137‑4 (MQ=0)

   tgcgtgcAGCACTTCGCCCGTCCGCACATAGCGACCCCGCGAGGCCCGGGGCAGAACAACTGGGACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGAt                                                                                                                            <  1:3843185/134‑1 (MQ=1)

             actgagccccTCCGCACCGAGCTCCCCAGCGAGGCCCCGGGCAGAACAACGCCTACACCAGCGGTGCGTCCCTCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                                  <  2:3453295/134‑1 (MQ=1)

             acgtATCCCGTCCGCACATCGCTACCCAGCGATGCCCTTCGCAGAACAACTGGTACACCAGCGGTGCGTCCAGCCCGGTCCTCTCCTACTAAGCACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                                  <  2:4477630/137‑1 (MQ=1)

             acgtATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGGACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                                  <  2:1527800/137‑1 (MQ=1)

              ctgaTCCCGGCCGCACCTAGCGACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAAc                                                                                                                 <  1:2201495/137‑1 (MQ=0)

                             cagaGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTc                                                                                                  <  1:3031604/137‑1 (MQ=1)

                                 gcgaCCCAGCGATGCCCTTGGCAGAACAACTGGTACACCCGCGCTGCGTCCATCCCGGTCCGCTCGTACTAAGGACAGCTCCTCTCAAAATTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGaa                                                                                              <  1:168528/137‑1 (MQ=1)

                                 gcgaCCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGaa                                                                                              <  1:4430876/137‑1 (MQ=1)

                                          cgctGCCCGTGGCAGAACACCTGGTACACCAGCGGTGCGTCCAGCCCGGGCGTCTCGTGCTAAGGAGAGCGGCGCGCAACTTGCCTGCGCGCGCGACGGATCGGGCCCGCACTGTCTCACGACGTGCGGACCCCAGCTcg                                                                                     >  2:4186972/4‑140 (MQ=0)

                                              gcgcTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCGCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTGTGAACCCAGCTCGCGTa                                                                                 <  2:3952262/137‑1 (MQ=0)

                                              gcccggggCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTa                                                                                 <  2:4388865/134‑1 (MQ=1)

                                               ccgttGGCAGAACAACTGGTACAGCAGCGGTGCGTCCATCCCGGTCGTCTCGTACTAAGGACCGCTCCTCGCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCAAGACGTTCGGAACCCAGCTCGCGTAc                                                                                >  1:2968287/4‑140 (MQ=1)

                                                         aaaaacTTGTACACCAGCGGTTCGTCCATAAAGGAACTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTACTGAACCCAGCTCGCGTACCGCTTTAATg                                                                      >  1:640263/4‑140 (MQ=1)

                                                               tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGCCGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  2:4661319/137‑1 (MQ=1)

                                                               tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  1:5649845/137‑1 (MQ=1)

                                                               tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  2:4640963/137‑1 (MQ=1)

                                                               tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  1:467677/137‑1 (MQ=1)

                                                               tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGGACCGCTTTAATGGGCGaa                                                                <  2:4020454/137‑1 (MQ=1)

                                                               tgcgaccccAGCGCGGCGTCCATCCCCGTCCTCACGTACCAAGGACAGCTCCTCTCAAATTTCCTGCCCCCGCGACCGATAGGGACCGAACTGTCACACGACGTTCTGAACCCAGCTCGCGAACCGCTTTAATGGGCGaa                                                                <  1:3042538/133‑1 (MQ=1)

                                                               cgggacacCCGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATGTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  2:3300131/136‑1 (MQ=0)

                                                               acggaccccAGCGGTGCGTCCATCCCGCTCCGCTCGTACTAAGCACACCTCCTCGCAAAGTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                                <  2:3358332/133‑1 (MQ=1)

                                                                      cctgCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                         <  2:3512047/137‑1 (MQ=1)

                                                                      ccggCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                         <  2:3751760/137‑1 (MQ=1)

                                                                      cccgCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                         <  2:737406/137‑1 (MQ=1)

                                                                             ggcctccctcccGCTCCTCTCGGACTCAGGACAGCTGCTCTCAAATTTCCTGAGCCCGCGACGGATAGGGACCGAACTGGCTCACGACGTTCTGACCCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTg                                                  <  2:4352478/136‑1 (MQ=1)

                                                                                 tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGaa                                                                >  2:3903664/7‑122 (MQ=0)

                                                                                 tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGaa                                                                <  1:3903664/116‑1 (MQ=0)

                                                                                 tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGAc                                              >  1:1733360/7‑140 (MQ=0)

                                                                                 gccctCCCGGTCCTCTCGTACTCAGGACAGCTCCTCTCAAATGTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGAc                                              <  1:541993/136‑1 (MQ=0)

                                                                                            ccgctcGTACTAAGCACAGCTCCTCTCAAAGTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGCGACCGACTACAGcc                                   <  2:4627007/137‑1 (MQ=0)

                                                                                                            catctcctcTCTAATTTCCTGCGCCCGCGACGGATCGGGCCCGAACTGTCTCACGCCGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACGACAGCCCCAGGATGCGATGAGc                   <  1:5029822/137‑1 (MQ=1)

                                                                                                               ctgctctCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCACAGGATGCGATGAGCCGa                >  1:639076/4‑140 (MQ=1)

                                                                                                               ctgctcgCAAGTTTCCTGCGCCCGCGACGGGTAGGGACCGAACTGTCTCACGACGGTCTGCACCCAGCTCGCGTACCGCTTTAATGGGCGACCAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGa                >  2:1447534/4‑140 (MQ=1)

                                                                                                                    ctaaaaTTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTATCACGACGTTCTGAACCCATCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGATCCGACATCg           >  2:3707978/4‑140 (MQ=0)

                                                                                                                             ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa  <  2:1796060/137‑1 (MQ=1)

                                                                                                                             ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa  <  2:1462213/137‑1 (MQ=1)

                                                                                                                             ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa  <  2:1225802/137‑1 (MQ=1)

AGATGCTTTCAGCACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAA  >  NZ_CP034943/2978016‑2978280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`