BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,968,935:1	+C	14.3%	coding (1019/1179 nt)	PP_3498 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,968,935	1	.	C	14.3%	99.7 / 20.1	42	coding (1019/1179 nt)	PP_3498	transposase
Reads supporting (aligned to +/- strand): ref base . (21/15); new base C (4/2); total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

CAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGG  >  NC_002947/3968810‑3969081

ccccggtggcattacccggttggtcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCt                                                                                                                            <  2:1149098‑M1/41‑1 (MQ=255)

                      gtcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAagga                                                                                                       >  1:636740‑M1/88‑150 (MQ=255)

                         agaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggc                                                                                                    >  1:1833140‑M1/85‑150 (MQ=255)

                          gaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                     >  2:896540‑M1/84‑148 (MQ=255)

                          gaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                     <  1:896540‑M1/65‑1 (MQ=255)

                             gctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGgcggcg                                                                                                <  2:133651‑M1/70‑1 (MQ=255)

                             gctcatagccgcctcccgtttttcgttcggtgtgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGa                                                                                               >  1:220649‑M1/81‑150 (MQ=255)

                                 atagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCa                                                                                            >  1:133651‑M1/77‑150 (MQ=255)

                                     ccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTTCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAggctg                                                                                       <  2:1520462‑M1/78‑1 (MQ=255)

                                            cggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTTCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATc                                                                                >  1:1520462‑M1/66‑150 (MQ=255)

                                                              gCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTatat                                                              <  2:1930977/150‑1 (MQ=31)

                                                                  gatgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGc                                                          >  1:268254/1‑150 (MQ=32)

                                                                    tgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGccc                                                        <  1:61693/150‑1 (MQ=32)

                                                                       tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGt                                                     >  1:958120/1‑150 (MQ=32)

                                                                         cctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                                   >  1:813074‑M1/37‑150 (MQ=255)

                                                                         cAGGGCGCGCCTTCCGTATTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGACGCCTCTTGCGCTGGCTATCCCGCCTTGGGCCCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                                   <  1:580247/150‑1 (MQ=255)

                                                                             gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                               <  2:173061/150‑1 (MQ=32)

                                                                             atcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                               >  1:1149098‑M1/33‑150 (MQ=255)

                                                                                    ttCCGTCTTGAGCAAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCa                                        >  2:1443203/1‑150 (MQ=17)

                                                                                       cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                                   >  1:954613/1‑136 (MQ=32)

                                                                                       cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                                   <  2:954613/136‑1 (MQ=32)

                                                                                       cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGt                                     >  2:650440/1‑150 (MQ=32)

                                                                                         tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGGTGGCTAGCCCACCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGCAACATATAAGCCCTGTTACAGAGCGCACAGGttt                                   >  1:832122/1‑150 (MQ=11)

                                                                                                 cAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGc                           >  2:1624439/1‑150 (MQ=32)

                                                                                                 cAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGc                           >  2:1624372/1‑150 (MQ=32)

                                                                                                  aGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGCAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCa                          >  2:361124/1‑150 (MQ=17)

                                                                                                   gCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCaa                         >  2:1492087/1‑150 (MQ=32)

                                                                                                     tAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTg                                                      <  1:1535347/119‑1 (MQ=33)

                                                                                                     tAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTg                                                      >  2:1535347/1‑119 (MQ=33)

                                                                                                      aaGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAcc                      <  2:1482480/150‑1 (MQ=32)

                                                                                                         gCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                   >  2:973100/1‑150 (MQ=32)

                                                                                                          gcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGc                                                                       >  1:719891‑M1/4‑97 (MQ=255)

                                                                                                          gcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGc                                                                       <  2:719891‑M1/94‑1 (MQ=255)

                                                                                                           cTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacg                 >  1:190052/1‑150 (MQ=32)

                                                                                                            tGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacga                >  2:1548588/1‑150 (MQ=32)

                                                                                                              ATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGaaa              <  1:1021932‑M1/150‑1 (MQ=255)

                                                                                                                   gATCAGCCAGA‑TGATCGCCGCCTCTTGCACGGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg         <  1:973100/150‑1 (MQ=17)

                                                                                                                        gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGTCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggc    <  2:1021932/150‑1 (MQ=17)

                                                                                                                        gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATgg                                                            >  2:1071400/1‑94 (MQ=34)

                                                                                                                        gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATgg                                                            <  1:1071400/94‑1 (MQ=33)

                                                                                                                        gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggc    >  2:822547/1‑150 (MQ=32)

                                                                                                                          cAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCTGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg  <  2:1751339/150‑1 (MQ=17)

CAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGG  >  NC_002947/3968810‑3969081

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