BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	454,656	G→A	V119I (GTT→ATT)	set23 →	staphylococcal exotoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	454,656	0	G	A	96.8%	95.1 / ‑3.6	31	V119I (GTT→ATT)	set23	staphylococcal exotoxin
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (13/17); total (13/18)
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.71e-01

CTTTTTAACCAACAAAATCAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTGTTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGAAGAAGTGTC  >  CP000730/454521‑454792

cTTTTTAACCAACAAAATCAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATttagtattagt                                                                                                                                     <  1:84453/141‑1 (MQ=255)

      aaCCAACAAAATCAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTa                                                                                                                               <  1:367519/141‑1 (MQ=255)

                aTCAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGaaaa                                                                                                                     <  1:71041/141‑1 (MQ=255)

                  cAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAAc                                                                                                                   >  2:242655/1‑141 (MQ=255)

                    aaaGTTCCAAGAGTTCGTATTGGGAAAATATGAAAATAAATACAAAGAAAAAACACATGGATTATATGTCTTTGCGGTACCAGAATTAGTAGATGTAGAAGGAAGAATACTTAGTATTATTGGTGTAACACAGAAAAACGt                                                                                                                 >  2:71112/1‑141 (MQ=255)

                          cgaaGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATc                                                                                                           <  2:301304/139‑1 (MQ=255)

                          ccAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATc                                                                                                           >  1:317431/1‑141 (MQ=255)

                                gtTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATAtt                                                                                                     <  1:116406/141‑1 (MQ=255)

                                        tGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGtctc                                                                                             <  1:355576/141‑1 (MQ=255)

                                           gAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTGGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTaa                                                                                          <  1:71112/141‑1 (MQ=255)

                                            aaaaGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAg                                                                                         <  1:28201/141‑1 (MQ=255)

                                                          aataCAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTAc                                                                           >  2:207533/1‑141 (MQ=255)

                                                             aCAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAg                                                                        <  2:146799/141‑1 (MQ=255)

                                                                 aGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTGTTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTaa                                                                    <  1:358995/141‑1 (MQ=255)

                                                                      aaaCACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAGCTTACTAGTTaaaaaaa                                                               >  1:273418/1‑141 (MQ=255)

                                                                          acaTGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGa                                                           >  1:81314/1‑141 (MQ=255)

                                                                                   aGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGa                                                  >  1:368646/1‑141 (MQ=255)

                                                                                   aGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGa                                                  <  1:279764/141‑1 (MQ=255)

                                                                                      tGTCTTTCCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACgg                                               <  1:334410/141‑1 (MQ=255)

                                                                                       gtatatgtggTACCAGAATTAGTAGATTTTGATGGAACAATATATAGTATTAGTGGTGTAACAAAGAAAAACGTAAAAACAATATTTGAGTCTCTAAGAACGCCGAACATAAAAGATAAAAAAATAGACGATAAAGACGGt                                              <  2:251044/133‑1 (MQ=255)

                                                                                            tGCGGTAACAGAAAAAGTAGTCTTAGATAGCAGCATATTTGGTAGTGGTGGTGTAACAAAGAAAAACGTAAAATCAGTAATTGAGAATCTAAGAACGTCGAAGTTACTAGTTAAAAAAATAGACGAGAAAGACGGtttttt                                         >  2:868/1‑140 (MQ=255)

                                                                                                gTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTAtt                                     <  2:317431/141‑1 (MQ=255)

                                                                                                    cAGAATTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATtgatg                                 <  2:289904/141‑1 (MQ=255)

                                                                                                        aTTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAAtt                             >  2:50667/1‑141 (MQ=255)

                                                                                                        aTTAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAAtt                             >  2:238202/1‑141 (MQ=255)

                                                                                                         ttAGTAGATTTAGATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAAttt                            <  1:158067/141‑1 (MQ=255)

                                                                                                                  ttAGATGGAAGAATATTTAGGATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACAAAATAAATAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTAtt                   <  2:129389/141‑1 (MQ=255)

                                                                                                                     gATGGAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCaa                >  1:62422/1‑141 (MQ=255)

                                                                                                                       tGGAAGAATATTTAGTAGTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCaaaa              <  2:125906/141‑1 (MQ=255)

                                                                                                                        ggAAGAATATTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAg             >  1:205453/1‑141 (MQ=255)

                                                                                                                               tatTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGaagaag      <  2:10509/141‑1 (MQ=255)

                                                                                                                                atTTAGTATTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGAAGAAgt     >  1:175785/1‑141 (MQ=255)

                                                                                                                                  ttagtattagtgGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGAAGAAgtgt   >  2:112950/1‑141 (MQ=255)

                                                                                                                                   tagtattagtgGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGAAGAAGTGTc  >  2:321677/1‑141 (MQ=255)

CTTTTTAACCAACAAAATCAAAAGTTCCAAGTGTTTTTATTGGGAAAAGATGAAAATAAATACAAAGAAAAAACACATGGTTTAGATGTCTTTGCGGTACCAGAATTAGTAGATTTAGATGGAAGAATATTTAGTGTTAGTGGTGTAACAAAGAAAAACGTAAAATCAATATTTGAGTCTCTAAGAACGCCGAACTTACTAGTTAAAAAAATAGACGATAAAGACGGTTTTTCTATTGATGAATTTTTCTTTATTCAAAAGGAAGAAGTGTC  >  CP000730/454521‑454792

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