BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_004302	13,763	T→C	V122A (GTT→GCT)	r1tp24 →	ORF24

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_004302	13,763	0	T	C	95.7%	53.0 / ‑2.6	23	V122A (GTT→GCT)	r1tp24	ORF24
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (19/3); total (20/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAAGATGAAGTGAATACAAATCCTAATAAGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGT‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATCCACTT  >  NC_004302/13693‑13830

tggcaTGAATTGAATACAAATCCTAATAAGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCa                                                                 >  1:11134736/5‑75 (MQ=255)

 tttttgaattgaatACAAATCCTAATAAGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCAt                                                                >  1:7494234/5‑75 (MQ=255)

               cctattgcTAAGAAGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGaacaa                                                  >  1:12564568/8‑75 (MQ=255)

                    gtgaaataatACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCt                                             >  1:13823121/5‑75 (MQ=255)

                    aagtaataaGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCt                                             >  1:10510840/4‑75 (MQ=255)

                      tctgtaaGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTAt                                           >  1:9950901/5‑75 (MQ=255)

                       taataaGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAacgacaccctcc                                          <  1:6673501/75‑7 (MQ=255)

                       aggagaGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACACCTAtt                                          >  1:2977171/6‑75 (MQ=255)

                       acaggggACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTAtt                                          >  1:7517544/7‑75 (MQ=255)

                        aataaGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACAtaacaac                                         <  1:12323933/75‑7 (MQ=255)

                             gttgagcACTATGTATCAAGAAGCTACAAGCAAAGGGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACaa                                    >  1:12345311/7‑75 (MQ=255)

                                tctcACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACatcaa                                 >  1:13814490/1‑71 (MQ=255)

                                 cttaatATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAaatat                                >  1:17520208/6‑75 (MQ=255)

                                 ctcACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGCCaaaata                                <  1:17615278/75‑4 (MQ=255)

                                  tgattgagTATCAAGAAGCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAaatata                               >  1:16211623/8‑75 (MQ=255)

                                           atctgttagCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATg                      >  1:13615865/8‑75 (MQ=255)

                                           agtaattagCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTACTAGACAAATATAATATCGATg                      >  1:14340490/8‑75 (MQ=255)

                                                 aagCTACAAGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTtata                >  1:15912480/1‑75 (MQ=255)

                                                      gataGCAAAGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACaa           >  1:9221422/4‑75 (MQ=255)

                                                      atgttttaaGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACaa           >  1:17590237/8‑75 (MQ=255)

                                                      aaatcgtaaGAGTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACaa           >  1:16989653/8‑75 (MQ=255)

                                                          ttagagagTTAGC‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATc       >  1:2586971/5‑75 (MQ=255)

                                                            tcagtattaga‑tGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATCCa     >  1:2083495/7‑75 (MQ=255)

                                                            gtatagttagc‑tGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATCCa     >  1:7182931/5‑75 (MQ=255)

                                                            gccgagTTAGTCTGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATcc      >  1:6846281/4‑75 (MQ=255)

                                                               tggtcagt‑tGCATTAAAACTTGAACAACATCTATTAGTCAAATATAATATCGATGTTTATAAACAAAATCCACtt  >  1:6119960/6‑75 (MQ=255)

                                                               tggtcagt‑tGCATTAAAACTTGAACAACATCTATTAGTCAAATATAATATCGATGTTTATAAACAAAATCCACtt  >  1:6139735/6‑75 (MQ=255)

AAAGATGAAGTGAATACAAATCCTAATAAGACTCTCACTATGTATCAAGAAGCTACAAGCAAAGAGTTAGT‑TGCATTAAAACTTGAACAACATCTATTAGACAAATATAATATCGATGTTTATAAACAAAATCCACTT  >  NC_004302/13693‑13830

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