BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,634	A→C	T108P (ACC→CCC)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,634	0	A	C	97.5%	106.4 / NA	40	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (21/18); total (22/18)
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.99e-01

CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAT  >  NC_000913/3041494‑3041769

cccATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGCGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTGTTCTCTAAAGTGccc                                                                                                                                       >  1:80997/1‑143 (MQ=255)

     tGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGCCGGAGAATAACTTGTTTGGTTTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATcgc                                                                                                                                  <  1:580992/143‑1 (MQ=255)

     ggCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGAGAAAACTATTTGGTGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATcgc                                                                                                                                  <  1:581855/142‑1 (MQ=255)

       cGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCgcgc                                                                                                                                <  2:407971/143‑1 (MQ=255)

           gCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAGGTGCCCATCGCGCCAga                                                                                                                            >  1:103130/1‑143 (MQ=255)

           gCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATACGGAATTCTCTGAAGTGCCCCTCGCGACAga                                                                                                                            >  2:445160/1‑143 (MQ=255)

              aaaGGTGAAGTGTGGGCCAATTTACGTCGGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATGCTGTAAAGGCCCCATCGCGCCAgacga                                                                                                                         <  2:699948/143‑1 (MQ=255)

                       gtgtgGAGCAATTTACGTCTGTTCCGCGCCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCgtgt                                                                                                                <  2:221914/142‑1 (MQ=255)

                             aGCAATTTACGGCTGTTCCGCGACGGCGATGGCTTTGCATGGATGGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                          <  2:689748/143‑1 (MQ=255)

                               cAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                        >  1:433177/1‑143 (MQ=255)

                                        tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                               >  1:433959/1‑143 (MQ=255)

                                        gcTGTTCCGCGCCGGCGATGGCTGTGCATGGATTGGACGGCGGAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                               <  2:712574/142‑1 (MQ=255)

                                          tGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGACCTGAAACAATATGCGGTATTCTCTACAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGt                                                                                             >  1:309394/1‑143 (MQ=255)

                                            ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGGGCTGCTTGGTGTTGCCGGttt                                                                                           >  2:224486/1‑143 (MQ=255)

                                              ccGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                         >  1:179810/1‑143 (MQ=255)

                                              ccGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                         >  1:738233/1‑143 (MQ=255)

                                                                 gCATGGCTTGACCGGCTCACCGTGCGTGAACCGCAGCTGACTGGTCTGGAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGc                                                                      <  2:558166/143‑1 (MQ=255)

                                                                 gCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCCTCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGGTTTCAGGCGCGCGCCGCGCTGGc                                                                      >  2:220966/1‑143 (MQ=255)

                                                                   aTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCGTGCGCGCGCCGAGGTGGCac                                                                    >  1:244649/1‑142 (MQ=255)

                                                                    tGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaaa                                                                   <  1:469443/143‑1 (MQ=255)

                                                                     ggATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGCGCGTGTGCTGCTTGGTGTTGCCGGCTTTCAGGCGCGCGCCGCGCTGGCaaaa                                                                  >  2:250130/1‑142 (MQ=255)

                                                                        ttGAACGGCGCAGCGTGCGTGAACCGCAGCTGAATGAACTGAAAAACTATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCATGGGGATGACAGTTGTCAAGCGCGCGCCGCGATGGCAAAtctc                                                               >  2:411260/1‑143 (MQ=255)

                                                                                cgcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTGTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa                                                       >  1:342825/1‑143 (MQ=255)

                                                                                cgcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGACAAATCTCGTTAGCgg                                                       >  2:36861/1‑142 (MQ=255)

                                                                                             aaCCGCAGCTGACTGAACTGCAAAAAGATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                          <  1:42919/143‑1 (MQ=255)

                                                                                              aCCCCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaag                                         <  1:268188/143‑1 (MQ=255)

                                                                                                   aGCTGACGGAACGGAAAAAATATGCGGTCTTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCTGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGaaaaa                                    <  1:99668/143‑1 (MQ=255)

                                                                                                      tcactgaactgcAACAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAg                                 <  1:195977/141‑1 (MQ=255)

                                                                                                       gactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGCCGATCGCGCCAGACGACGGGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAGACCTCTTTAGCGAACTGCCTTCGAAAGAAACACAgg                                >  2:97431/1‑143 (MQ=255)

                                                                                                        acggaacggaaAAAATATGCGGGATGCTCGAAAGTGCCGAGCGCGCCAGACGGCGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGt                               <  1:452310/140‑1 (MQ=255)

                                                                                                          tgaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAg                             <  1:655081/143‑1 (MQ=255)

                                                                                                          tgaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAg                             >  1:5480/1‑143 (MQ=255)

                                                                                                           gaactgaaAAAATATGCGGTATTCTCTAAAGTGCCGATCGCGACAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTCACAGGAGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGt                            <  2:94880/143‑1 (MQ=255)

                                                                                                            aactgaaGAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTc                           <  2:152401/143‑1 (MQ=255)

                                                                                                                           cGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGCGTTGCCGGTTTTCAGGCGCGCGCCGCGCTCGCAAACATCTTTAGCGGACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGAcc            >  2:41071/1‑143 (MQ=255)

                                                                                                                                ttCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGCGTTGCCGGTTTTCACGCGCGCGCCGCGCCGGCAAATCTCTTTAGCGANCTGCCTTCGAAAGAACAAAAGGTAGTGACGGCAGCCGCGACCAGttt       >  2:142786/1‑143 (MQ=255)

                                                                                                                                   tctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGACAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCt    >  2:492678/1‑143 (MQ=255)

                                                                                                                                   tctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGGTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGGGCGACCGCGTTGcg    >  1:43445/1‑142 (MQ=255)

                                                                                                                                    ctAAAGTGCCCATCGCGCCAGACGACGAGCGTGGGCTGCTTGGTTTTTCCGGTTTTGAGGCGCGCGGCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTGAAAGCAGGCGCGACGACgttccga   >  1:589490/1‑139 (MQ=255)

                                                                                                                                     tAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTCTAAGGAGCGAGCAGCGCTGGCAAATCTATTTAGGGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGCAGGCGCGAGCACGTTGCTac  >  2:106708/1‑142 (MQ=255)

                                                                                                                                     gaaaGTCCCCAGCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAt  <  1:416755/142‑1 (MQ=255)

CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAT  >  NC_000913/3041494‑3041769

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