BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,935,427	+G	noncoding (260/1542 nt)	BW25113_RS19520 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,935,423	1	.	G	86.8%	66.1 / 12.6	38	noncoding (256/1542 nt)	BW25113_RS19520	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base . (1/4); new base G (16/17); total (17/21)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

TCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCA  >  NZ_CP009273/3935290‑3935523

tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑ggggg                                                                                                 >  1:653362/1‑138 (MQ=17)

    ggAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑gggggt                                                                                                <  2:244409/136‑3 (MQ=17)

       aaCTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCt                                                                                          <  2:789129/139‑1 (MQ=32)

         cTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCt                                                                                          >  2:531336/1‑138 (MQ=17)

                   gAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGcga                                                                               <  2:140425/139‑1 (MQ=17)

                    aGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGcgac                                                                              <  1:961186/139‑1 (MQ=17)

                    aGGGGGATAACTACTGGAAACAGTAGCTAATACCGCATATCGTCGCAAGACCAAAGAGGGGGACCTCCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGcgac                                                                              <  2:1101325/139‑1 (MQ=11)

                          aTAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATccc                                                                        >  1:417941/1‑139 (MQ=17)

                             actactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAg                                                                     <  2:1168885/139‑1 (MQ=17)

                                    gAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCt                                                              <  1:103480/139‑1 (MQ=17)

                                      aaCGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑gggggt                                                                                                <  2:303322/102‑3 (MQ=17)

                                      aaCGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑gggggt                                                                                                >  1:303322/1‑100 (MQ=17)

                                         ggTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGagag                                                         >  1:275101/1‑139 (MQ=17)

                                           tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGa                                                       <  1:1188565/139‑1 (MQ=17)

                                                aaTACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGAc                                                   >  2:736271/1‑137 (MQ=32)

                                                   aCCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               >  2:284678/1‑139 (MQ=17)

                                                    ccGCATATCGTCGCAAGACCAAAGAGGGTGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATTGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               <  2:342778/138‑1 (MQ=11)

                                                    ccGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCGTCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               <  2:205484/138‑1 (MQ=12)

                                                      gCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCca                                             >  2:228903/1‑138 (MQ=17)

                                                       cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGt                                                                                                       >  1:953492/1‑79 (MQ=34)

                                                       cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGt                                                                                                       <  2:953492/79‑1 (MQ=34)

                                                       cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCcaca                                           >  2:569179/1‑139 (MQ=17)

                                                         tAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                                         <  1:569179/139‑1 (MQ=17)

                                                            cGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTgg                                       >  1:1138821/1‑138 (MQ=17)

                                                                 cAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACt                                   >  1:938451/1‑137 (MQ=17)

                                                                 cAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACt                                   >  1:355391/1‑137 (MQ=17)

                                                                    gACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                             <  1:736271/139‑1 (MQ=32)

                                                                     aCCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                             <  1:230261/139‑1 (MQ=17)

                                                                     aCCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                             <  2:275101/139‑1 (MQ=17)

                                                                     aCCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                             <  1:506396/139‑1 (MQ=17)

                                                                         aaGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGt                         >  1:612927/1‑139 (MQ=17)

                                                                         aaGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGt                         >  1:471382/1‑139 (MQ=17)

                                                                          agagGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               <  1:661963/116‑1 (MQ=12)

                                                                          agagGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               >  2:661963/1‑116 (MQ=12)

                                                                            agGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCa                      <  1:52729/139‑1 (MQ=17)

                                                                                   ccTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTcc                >  2:526701/1‑138 (MQ=17)

                                                                                     ttCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               >  2:811523/1‑105 (MQ=12)

                                                                                     ttCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                               <  1:811523/105‑1 (MQ=12)

                                                                                      tCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg            >  2:1127025/1‑139 (MQ=17)

                                                                                      tCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg            >  2:970067/1‑139 (MQ=17)

                                                                                       cGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg            <  2:386064/137‑1 (MQ=32)

                                                                                                 gCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGGG‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcagca  <  2:613119/138‑1 (MQ=12)

                                                                                                 gCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcagca  <  2:612927/138‑1 (MQ=17)

                                                                                                 gCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGTGGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcag    <  2:355391/136‑1 (MQ=17)

TCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTA‑GGT‑GGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCA  >  NZ_CP009273/3935290‑3935523

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