BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,729,168	A→G	10.2%	F525F (TTT→TTC)	U712_14070 ←	UvrABC system protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,729,168	0	A	G	10.2%	103.3 / 3.5	39	F525F (TTT→TTC)	U712_14070	UvrABC system protein C
Reads supporting (aligned to +/- strand): ref base A (13/22); new base G (2/2); total (15/24)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAAGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAGTAAATTT  >  CP006881/2729036‑2729279

cTCAAGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACt                                                                                                            <  7:238256/138‑1 (MQ=255)

cTCAAGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACt                                                                                                            >  3:219428/1‑138 (MQ=255)

cTCAAGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACt                                                                                                            >  5:151868/1‑138 (MQ=255)

   aaGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGa                                                                                                          <  8:277195/137‑1 (MQ=255)

    aGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTg                                                                                                       <  6:77924/139‑1 (MQ=255)

      gcttgcttCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCa                                                                                                     <  6:189516/139‑1 (MQ=255)

           cttccttcATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCa                                                                                                     <  2:220823/134‑1 (MQ=255)

           cttccttcATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCa                                                                                                     >  1:220823/1‑134 (MQ=255)

                    attttttAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCt                                                                                        <  7:110868/137‑1 (MQ=255)

                            aCGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGAAGTCGTCCAAAACGGAATGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGaa                                                                               <  1:292663/139‑1 (MQ=255)

                              ggcaCCGACATGCTTTAACAGCATTTTCTTTCTTTTCTCTCAGATACGCGGGATGGCGTCTAAAACTGATTGAAACGCGCTTTTTCCGCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATg                                                                             <  4:130818/136‑1 (MQ=255)

                              ggAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATg                                                                             <  4:192506/139‑1 (MQ=255)

                                       aTGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCagg                                                                    <  5:171339/139‑1 (MQ=255)

                                          cTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCaggagg                                                                 >  8:53797/1‑139 (MQ=255)

                                          cTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTCGCCTGTGAAAACTGa                                                                                                          <  6:132401/98‑1 (MQ=255)

                                          cTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTCGCCTGTGAAAACTGa                                                                                                          >  5:132401/1‑98 (MQ=255)

                                             tAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTa                                                               <  7:264043/138‑1 (MQ=255)

                                                            tcttttctCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGtt                                                <  4:104682/138‑1 (MQ=255)

                                                             cttttctCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTc                                              >  6:71129/1‑139 (MQ=255)

                                                             cttttctCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACATCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTc                                              >  4:165346/1‑139 (MQ=255)

                                                                        gATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATAc                                   >  4:258225/1‑139 (MQ=255)

                                                                           aCCCGGGATGTCGTCCAACACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGcc                                <  3:262154/139‑1 (MQ=255)

                                                                            cccGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGcc                                <  8:203226/138‑1 (MQ=255)

                                                                            cccGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGcc                                <  2:160699/138‑1 (MQ=255)

                                                                               gggATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCt                            >  7:49316/1‑139 (MQ=255)

                                                                                ggATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTc                           <  7:181371/139‑1 (MQ=255)

                                                                                    gtcgtcCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGc                       >  6:265960/1‑139 (MQ=255)

                                                                                       gtcCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGGAAACTGa                                                                                                          <  7:162231/53‑1 (MQ=255)

                                                                                       gtcCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGGAAACTGa                                                                                                          <  7:166659/53‑1 (MQ=255)

                                                                                       gtcCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGGAAACTGa                                                                                                          >  8:166659/1‑53 (MQ=255)

                                                                                       gtcCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGGAAACTGa                                                                                                          >  8:162231/1‑53 (MQ=255)

                                                                                             aaCGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCAcc               >  2:222424/1‑138 (MQ=255)

                                                                                               cGGATTGAATCGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTGGTCCTGAATGCGCTGCAGGAGGTTAAATTCCTGGCTGTTTCGTTCCAGATACGCCCCCTCCAGCGGATCCCCGAt            >  4:31475/1‑139 (MQ=255)

                                                                                               cGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGAt            <  2:143630/139‑1 (MQ=255)

                                                                                               cGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGAt            <  4:219428/139‑1 (MQ=255)

                                                                                                  aTTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAg         <  3:165346/139‑1 (MQ=255)

                                                                                                  aTTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAg         <  5:71129/139‑1 (MQ=255)

                                                                                                      aaaCGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAAGTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAGTaa      <  7:53797/138‑1 (MQ=255)

                                                                                                           cgcTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAATCGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAGTAAAttt  >  2:218581/1‑137 (MQ=255)

CTCAAGGCTTGCTTCCTTCATTTTTTTAACGGAACCGAAATGCTTTAACAGCATTTTCTTTCTTTTCTCTCCGATACCCGGGATGTCGTCCAAAACGGATTGAAACGCGCTTTTTCCCCGGATTTGCCTGTGAAAACTGATTGCAAAACGGTGCACCTCGTCCTGAATGCGCTGCAGGAGGTAAAATTCCTGGCTGTTTCGTTCCAGATACGCCACCTCCAGCGGATCACCGATCAGTAAATTT  >  CP006881/2729036‑2729279

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