BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,524,660	A→T	H73Q (CAT→CAA)	xapA ←	xanthosine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,524,660	0	A	T	97.3%	120.4 / ‑3.5	37	H73Q (CAT→CAA)	xapA	xanthosine phosphorylase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (23/13); total (24/13)
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.87e-01

CAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCTAAAATAAA  >  NC_000913/2524538‑2524791

cAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTc                                                                                                                    <  2:154967/140‑1 (MQ=255)

 aGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTcc                                                                                                                   <  1:427804/140‑1 (MQ=255)

     cAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCg                                                                                                               <  2:437918/140‑1 (MQ=255)

        tAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCTCCCGCAt                                                                                                            >  1:441461/1‑140 (MQ=255)

        tAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCAt                                                                                                            >  2:100849/1‑140 (MQ=255)

         aGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATg                                                                                                           <  1:416432/140‑1 (MQ=255)

            tCGCAGCCCAGCAGCATAAAGGTACTAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGAcc                                                                                                        <  1:194825/140‑1 (MQ=255)

                aGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGt                                                                                                    >  1:546318/1‑140 (MQ=255)

                         gcTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACtt                                                                                           <  2:30534/140‑1 (MQ=255)

                                 gTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATGGCCCAGCACTTACTCCCCCGCATGAACATGTAGCGTACTTACCGGGaa                                                                                   >  1:209723/1‑140 (MQ=255)

                                           cGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGc                                                                         >  2:95423/1‑140 (MQ=255)

                                                        ggTCATTCCACGACCTTCGTAGAAACGTCCGCGACCTGTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGaaa                                                            <  1:322705/140‑1 (MQ=255)

                                                           cATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTg                                                         >  2:136198/1‑140 (MQ=255)

                                                           cATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTg                                                         >  1:159886/1‑140 (MQ=255)

                                                            aTTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGc                                                        >  1:482974/1‑140 (MQ=255)

                                                            aTTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGc                                                        >  2:431211/1‑140 (MQ=255)

                                                             ttCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCg                                                       >  1:473073/1‑140 (MQ=255)

                                                                 aCGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   <  2:159886/140‑1 (MQ=255)

                                                                 aCGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   <  1:493601/140‑1 (MQ=255)

                                                                      cTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   <  1:167638/135‑1 (MQ=255)

                                                                      cTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   >  2:167638/1‑135 (MQ=255)

                                                                        tCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTtctc                                            >  2:354981/1‑140 (MQ=255)

                                                                        tCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTtctc                                            >  1:429906/1‑140 (MQ=255)

                                                                                   tCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   >  2:138706/1‑122 (MQ=255)

                                                                                   tCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAg                                                   <  1:138706/122‑1 (MQ=255)

                                                                                   tCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCgg                                 >  2:612650/1‑140 (MQ=255)

                                                                                   tCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCgg                                 >  2:649354/1‑140 (MQ=255)

                                                                                            tttCATACATACCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCCTTTCGTAGAAAATTGCCACAGCGTTATCATACTGATCGGCCAgcgcgc                        >  2:233387/1‑140 (MQ=255)

                                                                                                     tacCACCGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGAAAATTGCGACAGGGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCgg               >  2:265889/1‑140 (MQ=255)

                                                                                                      accaccGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGa              >  1:408943/1‑140 (MQ=255)

                                                                                                      accaccGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGa              <  1:612650/140‑1 (MQ=255)

                                                                                                      accaccGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACACCGTTCTCTATCTGACCGGCCAGCGCGCCGAGCCCGGa              >  1:300201/1‑140 (MQ=255)

                                                                                                        caccGGTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAAc            >  1:36264/1‑140 (MQ=255)

                                                                                                          ccGGTACGCCCTGGAGTTGGCCCAGCACTAACTCTCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCt          <  2:441461/140‑1 (MQ=255)

                                                                                                             gTACGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCTaaa       >  1:465620/1‑140 (MQ=255)

                                                                                                                cGCCCTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCTAAAATa    >  1:484332/1‑140 (MQ=255)

                                                                                                                  cccTGGAGTTGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCTAAAATaaa  >  1:584154/1‑140 (MQ=255)

CAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATGGCCCAGCACTAACTCCCCCGCATGACCATGTACCGTACTTACCGGGAAGCCCGGCAGCTTTTCGTAGGAAATTGCGACAGCGTTCTCAATCTGATCGGCCAGCGCGCCGAGCCCGGAACCTAAAATAAA  >  NC_000913/2524538‑2524791

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