BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,968,928	C→T	16.2%	Q338* (CAG→TAG)	PP_3498 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,968,928	0	C	T	16.2%	72.1 / 9.0	37	Q338* (CAG→TAG)	PP_3498	transposase
Reads supporting (aligned to +/- strand): ref base C (19/12); new base T (4/2); total (23/14)
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.46e-01

CAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCG  >  NC_002947/3968810‑3969074

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                                                                                              <  1:896540‑M1/65‑1 (MQ=255)

                          gaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                              >  2:896540‑M1/84‑148 (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                                            <  2:954613/136‑1 (MQ=32)

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

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

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

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

                                                                                                 cAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGc                    >  2:1624439/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                                                                <  2:719891‑M1/94‑1 (MQ=255)

                                                                                                          gcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGc                                                                >  1:719891‑M1/4‑97 (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)

CAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCG  >  NC_002947/3968810‑3969074

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: `‑`