BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,735,139	Δ1 bp	11.9%	coding (36/3093 nt)	PP_3302 →	RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,735,139	0	A	.	11.9%	120.4 / 11.9	42	coding (36/3093 nt)	PP_3302	RND family transporter
Reads supporting (aligned to +/- strand): ref base A (16/21); new base . (3/2); total (19/23)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTGACCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCAC  >  NC_002947/3734995‑3735282

gTGACCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCATCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAc                                                                                                                                              <  2:307234/149‑1 (MQ=255)

    ccTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcga                                                                                                                                         <  1:205097/149‑1 (MQ=255)

       gCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgagc                                                                                                                                       >  8:236458/1‑148 (MQ=255)

            cTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcg                                                                                                                                    <  4:120912/146‑3 (MQ=255)

            cTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcg                                                                                                                                    >  3:120912/1‑144 (MQ=255)

               gCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGatcaatca                                                                                                                              >  8:30840/1‑149 (MQ=255)

                      gatcgTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                     >  2:72243/1‑134 (MQ=255)

                      gatcgTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                     <  1:72243/135‑2 (MQ=255)

                           tGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGa                                                                                                                  <  1:295375/149‑1 (MQ=255)

                             gCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGAtt                                                                                                                <  7:236458/149‑1 (MQ=255)

                                   ggTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGc                                                                                                            <  2:81293/147‑1 (MQ=255)

                                   ggTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCt                                                                                                           >  5:256151/1‑148 (MQ=255)

                                    gTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGa                                                                                                         <  6:140713/149‑1 (MQ=255)

                                     tGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGAt                                                                                                        <  7:30840/149‑1 (MQ=255)

                                     tGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGAt                                                                                                        <  5:147214/149‑1 (MQ=255)

                                        gcATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGTAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGc                                                                                                     >  8:230186/1‑149 (MQ=255)

                                                tgcGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCg                                                                                             >  2:332356/1‑149 (MQ=255)

                                                tgcGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCg                                                                                             >  8:213306/1‑149 (MQ=255)

                                                         gCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTg                                                                                     <  8:207203/148‑1 (MQ=255)

                                                              gAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCttt                                                                               <  6:256151/149‑1 (MQ=255)

                                                                  tCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTc                                                                            >  5:113316/1‑148 (MQ=255)

                                                                    aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGc                                                                                       <  1:316762/135‑1 (MQ=255)

                                                                    aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGc                                                                                       >  2:316762/1‑135 (MQ=255)

                                                                    aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCttt                                                                               <  5:137922/143‑1 (MQ=255)

                                                                    aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCttt                                                                               >  6:137922/1‑143 (MQ=255)

                                                                         aCAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTggg                                                                    <  6:113316/149‑1 (MQ=255)

                                                                            ggCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCg                                                                 <  2:60053/149‑1 (MQ=255)

                                                                             gCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCCCGCTGGCCTTTTTCAAGCTGGGCCGt                                                                >  8:234612/1‑149 (MQ=255)

                                                                              cTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCg                                                                                             >  1:168977/1‑119 (MQ=255)

                                                                              cTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCg                                                                                             <  2:168977/119‑1 (MQ=255)

                                                                                cAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTg                                                               >  2:261695/1‑147 (MQ=255)

                                                                                     aGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCagaaga                                                        <  1:332356/149‑1 (MQ=255)

                                                                                          aGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTgg                                                                     <  4:342686/131‑1 (MQ=255)

                                                                                          aGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTgg                                                                     >  3:342686/1‑131 (MQ=255)

                                                                                                    gCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGTAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGt                                         <  7:230186/149‑1 (MQ=255)

                                                                                                    gCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGt                                         >  2:162960/1‑149 (MQ=255)

                                                                                                     cACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTc                                        <  1:162960/149‑1 (MQ=255)

                                                                                                            caTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGa                                 <  7:234612/149‑1 (MQ=255)

                                                                                                                                  ttATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATcc                                                     >  4:37328/1‑107 (MQ=255)

                                                                                                                                  ttATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATcc                                                     <  3:37328/107‑1 (MQ=255)

                                                                                                                                      gatcgcTCGC‑GTACTGCGAGCGATCAATCAACCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGgc       >  7:265634/4‑149 (MQ=255)

                                                                                                                                           cTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTCATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGGAGATGACGATTATCACTGCCTGGCCCGGCGCCAc  >  8:290823/1‑149 (MQ=255)

GTGACCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCAC  >  NC_002947/3734995‑3735282

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