BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,632,437	(T)_7→6	intergenic (‑71/+8)	BW25113_RS08135 ← / ← BW25113_RS25620	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,632,431	0	T	.	90.0%	60.2 / 3.9	20	intergenic (‑65/+14)	BW25113_RS08135/BW25113_RS25620	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (7/11); total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTC  >  NZ_CP009273/1632344‑1632494

tCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttt                                                               <  2:829883/90‑1 (MQ=255)

tCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttt                                                               <  1:579467/90‑1 (MQ=255)

    cATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttttc                                                           >  2:815112/1‑89 (MQ=255)

     aTTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttttca                                                          <  2:342590/90‑3 (MQ=255)

         ggATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCatgat                                                     <  1:159823/90‑1 (MQ=255)

         ggATATTCACGCATGATGGTGTCAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCatgat                                                     <  2:155228/90‑1 (MQ=255)

           atatTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGt                                                   <  2:669225/90‑1 (MQ=255)

            tatTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGtt                                                  <  1:760859/90‑1 (MQ=255)

              ttCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTa                                                 <  1:832884/90‑1 (MQ=255)

               tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTAcc                                               <  1:128366/90‑1 (MQ=255)

               tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTAcc                                               <  1:830037/90‑1 (MQ=255)

                 aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCgt                                             <  2:198586/90‑1 (MQ=255)

                 aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCgt                                             <  2:58065/90‑1 (MQ=255)

                   gCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCgtag                                           <  1:711893/90‑1 (MQ=255)

                    cATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCgtagt                                          >  1:427504/1‑90 (MQ=255)

                             gtTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTa                                 >  1:517201/1‑90 (MQ=255)

                                  cTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCAt                            >  2:332512/1‑90 (MQ=255)

                                   tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                           <  1:538123/90‑1 (MQ=255)

                                   tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                           <  1:252709/90‑1 (MQ=255)

                                     cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttt                                                               <  1:590196/53‑1 (MQ=255)

                                     cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttt                                                               >  2:590196/1‑53 (MQ=255)

                                     cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGa                         >  2:69637/1‑90 (MQ=255)

                                           tCGCTGGCACTACCACTGCTGGTTCAGAGTTGAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGac                   >  1:552954/1‑90 (MQ=255)

                                                   aCTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGc            >  1:297611/1‑90 (MQ=255)

                                                          ctgctgGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGc    >  1:775566/1‑90 (MQ=255)

                                                            gctgGTTCAGAGTTAAAAAAACTATGA‑TTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTc  >  2:367762/1‑90 (MQ=255)

TCGCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTC  >  NZ_CP009273/1632344‑1632494

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