BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,728,110	T→C	D421G (GAT→GGT)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,110	0	T	C	80.7%	84.5 / 15.7	31	coding (1262/1455 nt)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (12/13); minor base . (3/3); total (15/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 = 9.40e-01

CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATG  >  NC_000913/3727974‑3728244

cTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAcccc                                                                                                                                      <  1:422846/140‑1 (MQ=255)

   aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACcccccc                                                                                                                                   <  1:191165/140‑1 (MQ=255)

        cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTa                                                                                                                              <  2:404541/140‑1 (MQ=255)

        cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTCGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTa                                                                                                                              <  1:69792/140‑1 (MQ=255)

          gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACg                                                                                                                            >  1:133525/1‑140 (MQ=255)

               aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAAt                                                                                                                       >  1:526855/1‑140 (MQ=255)

                                    cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACCGTAATCGAGCT‑GCTGACCGCTGatat                                                                                                 >  1:649805/1‑140 (MQ=255)

                                    cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATc                                                                                                >  1:828456/1‑140 (MQ=255)

                                       cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCCTACGGTAATCGAGCT‑GCTGACCGCTGATATCCg                                                                                              >  1:818368/1‑140 (MQ=255)

                                          ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCTTACGGAAATCGACCT‑CCTGACCCCTGATATCCGCCa                                                                                           >  1:376900/1‑140 (MQ=255)

                                          ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAg                                                                                          >  1:266738/1‑140 (MQ=255)

                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACg                                                                                 <  1:784768/140‑1 (MQ=255)

                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCCAGCT‑GCTGACCGCCGATCTCTGCCCGCATCTGCCg                                                                                 >  2:658373/1‑140 (MQ=255)

                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGGACGGTAATCGAACT‑GCTGACCGCCGCTATCCGCCAGCATCTTAcc                                                                                 >  2:306591/1‑139 (MQ=255)

                                                         cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                           >  2:144570/1‑140 (MQ=255)

                                                          ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                           >  1:771468/1‑140 (MQ=255)

                                                            caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                                         <  1:512150/140‑1 (MQ=255)

                                                              acaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                       <  1:230125/140‑1 (MQ=255)

                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccc                                                           >  1:612859/1‑139 (MQ=255)

                                                                            gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                         >  1:88734/1‑140 (MQ=255)

                                                                             aTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                        <  2:82050/140‑1 (MQ=255)

                                                                             aTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccccccc                                                        >  2:499113/1‑140 (MQ=255)

                                                                                 ctctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                   <  2:828456/140‑1 (MQ=255)

                                                                                    tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCCCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCa                                                 <  2:133525/140‑1 (MQ=255)

                                                                                                gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                    <  1:144570/140‑1 (MQ=255)

                                                                                                 aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                    <  1:658373/140‑1 (MQ=255)

                                                                                                    tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                 <  2:215781/140‑1 (MQ=255)

                                                                                                        ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑CCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                            <  2:266738/140‑1 (MQ=255)

                                                                                                          aGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                           >  2:810152/1‑140 (MQ=255)

                                                                                                                gctgcGCCCAGTGCTGGCCCCACACCCCCCCCCAAAAAGAAGCCCAGCTCTCTTA‑TACTCATATGCCCCTGCCTCCGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCgaa                     >  2:561144/1‑138 (MQ=255)

                                                                                                                    cgcCCAGTGCTGGCCCCACCCCCCCCCCCATACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                 <  1:76199/140‑1 (MQ=255)

                                                                                                                     gcCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                <  1:527084/140‑1 (MQ=255)

                                                                                                                         aGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc            >  2:512085/1‑140 (MQ=255)

                                                                                                                          gTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg           <  2:612859/140‑1 (MQ=255)

                                                                                                                                   ccACACCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg  <  2:88734/140‑1 (MQ=255)

CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCT‑GCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATG  >  NC_000913/3727974‑3728244

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