BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,986,622	T→G	38.5%	D31A (GAT→GCT)	CJ019_03763 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,986,622	0	T	G	38.5%	9.4 / 4.0	37	D31A (GAT→GCT)	CJ019_03763	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (10/10); new base G (6/9); total (18/19)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGGTTGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  ACB122/3986491‑3986752

aGGTTGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCgg                                                                                                                                 >  3:367050/1‑136 (MQ=255)

    tGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGg‑cgc                                                                                                                             <  2:66484/135‑1 (MQ=255)

     gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGGCCGG‑ggcg                                                                                                                            >  5:256384/1‑135 (MQ=255)

          cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGCACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGG‑CGCGGGCAc                                                                                                                       >  2:319943/1‑135 (MQ=255)

           agcaccACCTGAGCCAGCAACTTAGCCGGGTACACGGCATGGGGGGCCAGCAGCTGGTCCAGTTGCTCCTGGGACAACCCCGGGGCCTTGGCTACGCGCGTGGCAGCGGACGGCGCCGGAGCCGG‑CGCGGGCACg                                                                                                                      <  6:256384/130‑1 (MQ=255)

           aTCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCGG‑CGCGGGCCCg                                                                                                                      >  1:285757/1‑135 (MQ=255)

                ccccTGCGCCAGCAACATTTCAGGGTACAAGGCAAGGGGGGCCCGCATCTGGGCCTGTGGCTCCTGGGAGAACACGGGGTCCTTGGCTAGGTGGGTGTCAGGGTAAGGCGCCGGATCCGG‑CGCGGGCACGTCTGCg                                                                                                                <  4:213417/134‑1 (MQ=255)

                cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGG‑CGCGGGCCCGGCttc                                                                                                                 >  2:139025/1‑133 (MQ=255)

                cACCTGACCCATCAACTTCTCAGGGTACACTGCAATGGGGGCCAGCATCTGCTCAAGTTGCTCCGGGGAGAACACCGGGGCCTTGGCAACTAGCTTGGAAGCGGAAGGCGACTGATCCGG‑CGCGGGCACGTCTGc                                                                                                                 <  6:16467/135‑1 (MQ=255)

                        cctgCAACTGTTCAGGGTACAAGGAAATGGGGGCCTGCATCTGGTCCATTTGCTCCTGTGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGGGCCGGATCCGG‑CGCGGGCACGTCTGCGGCCTGc                                                                                                          <  2:64032/131‑1 (MQ=255)

                        ccAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCCGCGCCGCCGG‑CGCCGGCCCGGCTGCGGGCGGCAc                                                                                                        >  3:213417/1‑136 (MQ=255)

                           gaaaCTGATCCGTTTACACGCCAATGGGGCCCAGCAGGTGGACCAGTGGCTCCGGGCAGAACCCCGGGTCCTGGGCGACGAGCGGGGCCGCGCACGGCGCCGGATCCGG‑CGCGGGCACGTCTGCGGCCTGCAccag                                                                                                     <  4:102210/134‑1 (MQ=255)

                                 gagcAGGGCACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACCCCGGGCCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGAGCCGG‑CGCGGGCACGTCTGCGGCCTGCAccaggccagg                                                                                               <  8:87423/133‑1 (MQ=255)

                                     aGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCACCCG‑CGCGGGCACGTCAGCGGCCTGCACCACGCCCGGGGc                                                                                            >  2:347643/1‑135 (MQ=255)

                                                 aaTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGG‑CGCGGGCCCGTATGCGGCCGGCCCCCGGCCAGGGGCCGCGCGCAGGc                                                                                 >  6:225966/1‑134 (MQ=255)

                                                    gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAATGCGCCGGCTCAGGACGC‑GGCACCGGAGCCGCCGGAACACGCCCAGGGGCAGCGAGCAGGCAAAgg                                                                            >  8:140162/1‑135 (MQ=255)

                                                    gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGGTCCCG‑CGCCGGCCCCGCAGCGGCCGGCCCCCCGCCACGGGCCGCGCGCACGCcaac                                                                             >  4:323449/1‑134 (MQ=255)

                                                    gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGG‑CGCGGGGGCCGCGGCGGCCGGCCCCCCGCCCCGGGCGGCGCGCCGGCCAAg                                                                             >  1:145675/1‑135 (MQ=255)

                                                      gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGG‑CGCGGGCCCGTCTGCGGCCGGCCCCAGGCCCGGGGCCGCGCGCCGGCCAAGCAg                                                                          >  8:333331/1‑136 (MQ=255)

                                                      gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGG‑CGCgg                                                                                                                           <  8:365214/87‑1 (MQ=255)

                                                      gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGG‑CGCgg                                                                                                                           >  7:365214/1‑87 (MQ=255)

                                                                    tCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGGCCGC‑GGCGGGCCGGCCGGCGGCCGGCCCCAGGCCCCGCGCCGCGCGCAGGCCCCGCAGTGCCGAcccca                                                               >  8:2715/1‑129 (MQ=255)

                                                                     ccAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGGGCCCGGGCCGG‑CGCGGGGACGGCGGCGGCCGGCCCCCCGCCCAGGGCCGCCCGCCGGCAAATCAGTTGCGCCAGGaaa                                                             >  2:175303/1‑134 (MQ=255)

                                                                                        cccggggTCCGTGGCTCCGCGGGTGGCCGCGGACGGCGCCGGAGCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCAt                                         <  1:175303/131‑1 (MQ=255)

                                                                                             gggcgttGGCTACGCGGGTGGCCGGGGATGGCGCCGGATCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                    <  1:9665/130‑1 (MQ=255)

                                                                                                 gttGGCGACGCGCGGGCAAGGGGATGGCCGCGGAGCCGG‑CGCGGGCTCGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATcc                                <  5:105382/134‑1 (MQ=255)

                                                                                                 ctgggCTACGCGGGTGGCCGGGGAAGGCGGCGTTTCCGG‑GGCGGGCACGTCTGCGGCGTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCTTCGCTTATcc                                <  2:90904/132‑1 (MQ=255)

                                                                                                           aGCGTGGCAGCGGAAGGCGCGGGCTCAGG‑GGCCGGCCCGGCTGCGGCCGGCGCCCCGCCCGGGGCCGCGCGCACGCCACGCAGGGCCCGCACGCCAGCCAACAAAGACATACGCCCCGCGTATCCCCAACAGGTCa                     >  4:190870/1‑136 (MQ=255)

                                                                                                             cGTGGCAGCGGAAGGCGCGGGAACCCG‑CGCGGGCACGGCTGCGGCCTGCCCCACGCCCGGGGCAGCGCTCAGGCAAACCAGTTCCCAAACGAAAGCAAACAAATGAAAGCGCAACGCTTATCCCCATCAGGTCAg                    >  6:166021/1‑135 (MQ=255)

                                                                                                               tGGCCGCGGACGGCGCCGGCTCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  7:385800/136‑1 (MQ=255)

                                                                                                                ggctggggACGGCGCCGGCGCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCGAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  1:139025/129‑1 (MQ=255)

                                                                                                                ggCTGCGGAAGGCGCGGGAGCGGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACcg                <  8:232720/136‑1 (MQ=255)

                                                                                                                ggCAGCGGACGGCGCCGGATCCGG‑CGCGGGCACGGCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  2:169790/135‑1 (MQ=255)

                                                                                                                 cccgCGGAGGGCCCGGGAGCGGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACcg                <  2:285757/132‑1 (MQ=255)

                                                                                                                     cggccggCGCCGGAGCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg           <  3:369664/131‑1 (MQ=255)

                                                                                                                      ggcaGGCGCCGTAGCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg           <  6:233800/132‑1 (MQ=255)

                                                                                                                        aaGGCGCCGGGTCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACCGGAACGCAAACAATGGCATGCGCGTCGCTTATCCCCATGAGGTCAGACCGCGgggggt        >  7:314331/1‑134 (MQ=255)

                                                                                                                               gggctCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGt  <  1:347643/131‑1 (MQ=255)

AGGTTGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGG‑CGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  ACB122/3986491‑3986752

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