BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,353,513	G→A	13.3%	G227G (GGC→GGT)	CJ019_02227 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,353,513	0	G	A	13.3%	97.7 / 5.1	39	G227G (GGC→GGT)	CJ019_02227	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (19/14); new base A (3/2); total (23/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATGGCTGTGCTTCC  >  ACB143/2353396‑2353641

gcTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGcc                                                                                                                 >  3:227985/1‑135 (MQ=40)

          tGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTTCACAGAAGGGACa                                                                                                      >  8:262682/1‑136 (MQ=255)

                 gTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAt                                                                                               >  4:216378/1‑136 (MQ=40)

                 gTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCAAGCAGGGCCATTTCCAt                                                                                               >  7:102331/1‑136 (MQ=37)

                         aGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGagg                                                                                         >  7:87713/1‑134 (MQ=255)

                              ggTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCg                                                                                  >  5:111322/1‑136 (MQ=255)

                              ggTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCAATGTGAGGCAGGCCg                                                                                  >  5:1512/1‑136 (MQ=255)

                                   tCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCg                                                                             >  6:231244/1‑136 (MQ=255)

                                   tCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCg                                                                             >  2:82669/1‑136 (MQ=255)

                                   tCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCg                                                                             >  4:236400/1‑136 (MQ=255)

                                    cAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCgg                                                                            >  8:269686/1‑136 (MQ=255)

                                        cgggggcgggggACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTgg                                                                        >  4:163959/1‑136 (MQ=255)

                                         gggggcggcggACGGCGTGTGGGTCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCGGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTggg                                                                       <  5:248334/136‑1 (MQ=255)

                                             gcgggggACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGc                                                                   >  6:113351/1‑136 (MQ=255)

                                              ggggggACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGcc                                                                  <  4:254946/135‑1 (MQ=255)

                                                    aCGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGAc                                                              >  5:155718/1‑134 (MQ=255)

                                                       gCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGc                                                         <  4:246349/136‑1 (MQ=255)

                                                            tgcgCCTTCTTCCATATATCGACGGATGACTACCCACTCATCACGCAGTTCGAAGGCGCCGTGGCCTTGTTTGCCCAGCAGGGCCATTTCAATGCGCGGCAGGCCGGTGCGGGGGGGGTGCCAGACCAGGCCGAg                                                     <  8:77941/132‑1 (MQ=14)

                                                            tgGGCCTTCTTCCATATAACCACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGTCCATTTCCTAGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAg                                                     <  4:285296/135‑1 (MQ=255)

                                                               gCCTTCATCCATATTACGACGGATGAATTCCCACTGATCACTCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGGAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCcgc                                                 <  6:1512/136‑1 (MQ=17)

                                                                 cttcttCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGt                                               <  1:264813/136‑1 (MQ=255)

                                                                 cttcttCCATATAACGACGGATGAAATCCCACTCCTCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGTGGTGCCAGACCAGGCCGAGCCGCGt                                               <  1:245985/136‑1 (MQ=255)

                                                                 cttcatCCATAAAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGGCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGt                                               <  7:233734/136‑1 (MQ=255)

                                                                     ttcCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGccc                                           <  3:84786/136‑1 (MQ=255)

                                                                       ccATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGccc                                           <  8:234065/134‑1 (MQ=255)

                                                                       ccATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGATAGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCTGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGccc                                           <  1:211490/134‑1 (MQ=255)

                                                                             aaCGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCtgag                                   >  1:199557/1‑136 (MQ=255)

                                                                                  cgGATGAATTCCCCCTCATCACGCAGTTCGCTGCCACCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAAGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGcc                                                   <  5:29064/115‑1 (MQ=255)

                                                                                  cgGATGAATTCCCACTCATCACGCAGTTCGCTGCCCCCTAGGCCTTGTTTGCACAGCAGGGCCATTACCATGTGAGGCACGCCGCAGCCGCTCGGGTGCAAGACCAGGCCGACCCGCGTGCCCACGCTgagtgagt                              >  4:81916/1‑136 (MQ=18)

                                                                                  cgGATGAATTCCCACTCATCACGCAGTTCGCTGCCACCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAAGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGcc                                                   <  5:25508/115‑1 (MQ=255)

                                                                                  cgGATGAATTCCCACTCATCACGCAGTTCGCTGCCACCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAAGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGcc                                                   >  6:25508/1‑115 (MQ=255)

                                                                                  cgGATGAATTCCCACTCATCACGCAGTTCGCTGCCACCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAAGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGcc                                                   >  6:29064/1‑115 (MQ=255)

                                                                                    gATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCa                            >  2:277384/1‑136 (MQ=255)

                                                                                     aTGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCAt                           >  1:188035/1‑136 (MQ=255)

                                                                                       gAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCAtct                         >  7:230976/1‑136 (MQ=255)

                                                                                                   atcaCGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATg             <  2:296532/136‑1 (MQ=255)

                                                                                                      aCGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATGGcg          >  3:174361/1‑135 (MQ=255)

                                                                                                        gCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATGGCtgt        <  2:1635/136‑1 (MQ=255)

                                                                                                              cGCTGCCACCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATGGCTGTGCTTcc  >  1:35858/1‑136 (MQ=255)

GCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGTGCGGCTGGGGTGCCAGACCAGGCCGAGCCGCGTGCCCACGCTGAGTGAGTCATCTTCTGGGTCGATGGCTGTGCTTCC  >  ACB143/2353396‑2353641

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