BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,978,137	A→T	23.3%	noncoding (2664/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,137	0	A	T	23.3%	‑1.7 / 5.3	25	noncoding (2664/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (3/16); new base T (3/3); total (6/19)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGC  >  NZ_CP034943/2978029‑2978263

actgagcccGTCCGCACATAGCGACCCAGCGATGCCCTTGGCAGAACCACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCTCGACGGATAGGGACCGaa                                                                                                 <  1:2215988/134‑1 (MQ=1)

         gtacGCACATAGCTACCCAGCGATCCCCTTGGCAGAACAACTGGAACACCAGCGGTGCTTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCATCCGCCCGCGACGCATAGGGACCGACCTCTCTCac                                                                                        >  1:1426265/4‑140 (MQ=1)

            cgcatttAGCTACCGGGCAGTGCCATTGGCATGACAACCCGAACACCAGTGATGCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCacgac                                                                                     <  1:2885979/134‑1 (MQ=0)

            cgaacaTAGATACCCAGCGATGCCCTTGGCAGAACAACTGGTTCACCAGCGGTGCGTCCATACCGGTCCTATAGTACTAAGGACAGGTCCTCTCAATTTTCCTGCGCGCGCGAAGGATAGGTTCCGTACTGTCTCaagtc                                                                                     <  2:1279781/137‑5 (MQ=1)

                    gcgaCCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGaa                                                                             <  2:2567700/137‑1 (MQ=1)

                           agagATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCGCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAGCTGTCGCACGACGGTCGGACCCCCGCt                                                                      >  2:4638408/4‑140 (MQ=1)

                             cgttGCCCTTGGCAGAACAACTGGTAACCCAGCGGTGTGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTGGCAATTTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGACCCCCGCTcg                                                                    >  1:3848139/4‑140 (MQ=1)

                                 gcgcTTGGCAGAACACCTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGAAAGCTCCTCTGAAATTTCCTGCGCCCGCGACGGATAGGGACCGAAATGTCTCACTACGTTCTTTACGCAGCTCGCGTa                                                                <  2:4251235/137‑1 (MQ=1)

                                    ctgggCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTATGAACCCAGCTCGCGTACCg                                                             <  2:4011821/137‑1 (MQ=1)

                                                acgggTACACCAGCGGTGCGGACATCCCGCGCCTCTCGTACTAAGGACAGCTCCTCTCAAATGTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCg                                                 <  1:2861763/137‑1 (MQ=1)

                                                 ctcgTACACCAGAGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTATAAAATTTCCTGAGCCCGCGACGGATAGGTACCGAACTGTCTCACGACGTTCAGAACCCAGCTCGCGTACCGCTTTAATGGGCGa                                                <  2:2989813/137‑1 (MQ=1)

                                                  tgttACACCAGGGGTGGTTCATTCCCGGTGCTTTGGTACTAAGGACAGCTCGTCTCATTTTTCCTGCGCCGGCGACTGTTTGGGACCGAACTGTCTCCCGACTTTCTGAACCCAGCTGGCGTACCGCTTTAATGGGCGac                                               <  2:396643/137‑2 (MQ=0)

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

                                                  tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTAATAAGGACAGCTCCTATCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCACGTACCGCTTTAATGGGCGaa                                               <  1:2323008/137‑1 (MQ=0)

                                                  tgctACACCAGCGGTGCGTCCATCCCGCTCCTCTCGTACTAAGGACAGCTCCTCACAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                               <  1:2705954/137‑1 (MQ=0)

                                                         cccgCGGTGCGTCCATCCCGGTCCTCTCGGACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                        <  1:756484/137‑1 (MQ=1)

                                                               gttcGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGTACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCtt                                  >  2:1259980/4‑140 (MQ=1)

                                                                tgggTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCGTGCGCCCGCGACGGATAGGGACCGAACTGTCTCTCGACGTTCTGAACCCAGCGCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCtgg                                 <  2:4866648/137‑3 (MQ=1)

                                                                   gtgcATCCCGGTCCTCTCGTACTCAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGATCGCTTAGCTGTTTAATGTTCGAACAGCCCAACCCTTggtt                              <  2:2063325/137‑3 (MQ=1)

                                                                   gtacATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGa                              <  2:3265890/137‑1 (MQ=1)

                                                                    tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGAc                             >  2:2640480/7‑140 (MQ=0)

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

                                                                                 tcgcGTACTAAGGACAGCTCCTCTCAAAGTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGcccc                <  2:1303182/137‑1 (MQ=1)

                                                                                            ccaccgCTCCTCACATAGTCCCTGCGCCCGCGACGGATAGGGACCGAACTCACTCCCCCCGATCTGAACCCAGCTCGCGAACCGCTTGAATGGGCGAACAGCCCAACCCTGCGGACCGACTACCGCCCCAGGAGGCGATg     <  2:3007722/135‑1 (MQ=1)

                                                                                               catctcctcTCAAATTTCCTGCGCCAGCGACGGATAGGGACCGAACTGTCTCACTACGTTCTGAACCCAGCTCGCGGACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGc  <  2:2968639/137‑1 (MQ=0)

ACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGC  >  NZ_CP034943/2978029‑2978263

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