BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,174,064	G→C	H288D (CAT→GAT)	U712_10875 ←	CCA‑adding enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,174,064	0	G	C	95.5%	65.7 / NA	22	H288D (CAT→GAT)	U712_10875	CCA‑adding enzyme
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (13/8); minor base A (0/1); total (13/9)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTTGAATGTCCTTTAGTTTTTCATCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAAAGAAGTGATTTTTCCCCGGCTTTGTACATTGACATAGTGTCAAGGCTCTGGCCAAACGTATGGGAAATGTGAATCGCTTCCTTCATCACTTTTCCCGGAAGCTTCCATGCTTTCAGGAAAAGCGGCGCATCTTTCAAATCTATTTCAAGATCAATTAAAAGCGCAGCCCAAAGCTCTTCA  >  CP006881/2173929‑2174182

cgtactGAATGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATcgg                                                                                                                      >  3:94325/6‑138 (MQ=255)

 gtactGAATGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATcggca                                                                                                                    <  8:308515/135‑3 (MQ=255)

 gtactGAATGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATcggca                                                                                                                    <  7:293885/135‑3 (MQ=255)

 gtactGAATGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATcggca                                                                                                                    <  5:226415/135‑3 (MQ=255)

 gtactGAATGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATcggca                                                                                                                    <  3:266262/135‑3 (MQ=255)

         tGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAAt                                                                                                            >  6:290236/1‑139 (MQ=255)

         tGTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAAt                                                                                                            >  5:350233/1‑139 (MQ=255)

          gTCCTTCAGTTTCTTTTCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGAGGTTTTTTCCCAGCGCTGTACAGTGTCATTGCGTCACCGCTCTGGCCAACCGTATCGGCAAGGTGAAtg                                                                                                           >  1:215743/1‑138 (MQ=255)

                      cttttcatcCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTtgat                                                                                                <  7:2406/134‑1 (MQ=255)

                       ttttcatcCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTtgatca                                                                                              <  6:259403/136‑1 (MQ=255)

                                                 aaGCTGCGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTc                                                                    <  2:331649/139‑1 (MQ=255)

                                                       cGATATTTTGGCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGaa                                                               >  1:359579/1‑138 (MQ=255)

                                                                 gCCGCTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGgcgcgt                                                    >  1:42345/1‑137 (MQ=255)

                                                                     cTGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCttt                                                <  2:42345/139‑1 (MQ=255)

                                                                      tGAGAGCAGTGCTTTTTTCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTc                                               >  3:440333/1‑139 (MQ=255)

                                                                               tgcttttttCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTAtg                                      >  4:305824/4‑138 (MQ=255)

                                                                                 cttttttCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTAtgcc                                    <  3:305824/138‑4 (MQ=255)

                                                                                 cgattttCCCAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTAtgcc                                    <  8:197962/136‑4 (MQ=255)

                                                                                      ttcccAGCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATAGTCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAgg                                <  4:179756/136‑2 (MQ=255)

                                                                                           agCTCTGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGTTAAt                           >  3:386818/2‑138 (MQ=255)

                                                                                                tGTACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGttaattaa                        >  4:439254/1‑136 (MQ=255)

                                                                                                  tACATTGTCATTGCGTCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGTTAATTAAAAgag                   <  6:350233/139‑3 (MQ=255)

                                                                                                            ttgcgtCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGTTAATTAAAAGAGCAGCCCAAAg         >  7:166910/5‑139 (MQ=255)

                                                                                                            ttgcgtCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTAAAAAATATTACAAGGTTAATTAAAAGAGCAGCCCaaac         >  2:359990/5‑138 (MQ=255)

                                                                                                              gcgtCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGttaattaa                        <  2:104026/120‑1 (MQ=255)

                                                                                                              gcgtCAAGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTCTTTCAAAACTATGCCAAGGttaattaa                        >  1:104026/3‑122 (MQ=255)

                                                                                                                   aaGGCTCTGGCCAAACGTATCGGCAATGTGAATGGCTTCCTTGATCACTTTTCCCGGGAGCTTCCAGGCTTTCAGGAAAAGGGGCGCGTGTTTCAAAACTATGCCAAGGTTAATTAAAAGAGCAGCCCAAAGCTCCTCa  >  5:100129/1‑139 (MQ=255)

TATTTTGAATGTCCTTTAGTTTTTCATCATCCAGTTTTTCGTTTTGCCGAAGCTGCGATATTTTGGCCGCTGAAAGAAGTGATTTTTCCCCGGCTTTGTACATTGACATAGTGTCAAGGCTCTGGCCAAACGTATGGGAAATGTGAATCGCTTCCTTCATCACTTTTCCCGGAAGCTTCCATGCTTTCAGGAAAAGCGGCGCATCTTTCAAATCTATTTCAAGATCAATTAAAAGCGCAGCCCAAAGCTCTTCA  >  CP006881/2173929‑2174182

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