BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,028,269	G→C	100%	P126A (CCA→GCA)	USA300HOU_RS10230 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,028,269	0	G	C	100.0%	122.9 / NA	35	P126A (CCA→GCA)	USA300HOU_RS10230	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (12/23); total (12/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACAA  >  USA300TCH1516_ALE/2028134‑2028404

ccTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTAc                                                                                                                                    >  1:848068/1‑141 (MQ=255)

       tttACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTc                                                                                                                             <  1:835660/141‑1 (MQ=255)

        ttACCTCAACTGTGCCTAACAACAATGGTATTAATTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCa                                                                                                                            <  1:180027/141‑1 (MQ=255)

          aCCTAAAATGTGTCTAACAACAATTGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCAcc                                                                                                                          <  1:152001/141‑1 (MQ=255)

           ccTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCa                                                                                                                         <  1:828261/141‑1 (MQ=255)

                                aTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATaaaa                                                                                                    <  1:203201/141‑1 (MQ=255)

                                  ggTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGa                                                                                                  >  1:143298/1‑141 (MQ=255)

                                   gTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGAt                                                                                                 <  1:153708/141‑1 (MQ=255)

                                     attattTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATgg                                                                                               >  1:74950/1‑141 (MQ=255)

                                           ttCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGt                                                                                         <  1:519163/141‑1 (MQ=255)

                                                  atatTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAAcc                                                                                  >  1:526951/1‑141 (MQ=255)

                                                   tatTTGCAGCTTTTGTTAAGTCAAAATGCAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCt                                                                                 >  1:560451/1‑141 (MQ=255)

                                                         cAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCg                                                                           <  1:347570/141‑1 (MQ=255)

                                                                tGTTACGTCAATATGTACGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAac                                                                    <  1:143029/141‑1 (MQ=255)

                                                                 ggtAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAact                                                                   <  1:38467/139‑1 (MQ=255)

                                                                  ttAAGTCAATATGTAGGTCGCCAATTCAATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTAAGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGGTAAGGCATAACCTGATTCGTCGCacaga                                                                  >  1:81554/1‑137 (MQ=255)

                                                                   tAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactac                                                                 <  1:62990/141‑1 (MQ=255)

                                                                     aggcAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtt                                                               <  1:562187/138‑1 (MQ=255)

                                                                     aGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtt                                                               <  1:652574/141‑1 (MQ=255)

                                                                       tCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTAATCCGTCGCAACTACtttt                                                             >  1:398337/1‑141 (MQ=255)

                                                                         aaTATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtttttt                                                           <  1:468224/141‑1 (MQ=255)

                                                                           tatGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTg                                                         >  1:532253/1‑141 (MQ=255)

                                                                                  gTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTaaa                                                  >  1:703031/1‑141 (MQ=255)

                                                                                          ttCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTa                                          <  1:346104/141‑1 (MQ=255)

                                                                                          ttCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTa                                          <  1:266704/141‑1 (MQ=255)

                                                                                          ttCCACGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTa                                          <  1:319869/141‑1 (MQ=255)

                                                                                            ccATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTAtt                                        >  1:526790/1‑141 (MQ=255)

                                                                                             cATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTa                                       <  1:836914/141‑1 (MQ=255)

                                                                                               tGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTaaa                                     <  1:270793/141‑1 (MQ=255)

                                                                                                   gAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGGATTAAATAAc                                 >  1:803498/1‑141 (MQ=255)

                                                                                                    aaTTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAATCCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACg                                <  1:715181/141‑1 (MQ=255)

                                                                                                              tCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCTGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGt                      <  1:580732/141‑1 (MQ=255)

                                                                                                               cTTCCCACTTATATAAATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGta                     <  1:869644/141‑1 (MQ=255)

                                                                                                                                 aaaCTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACa   <  1:509066/141‑1 (MQ=255)

                                                                                                                                  aaCTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAACTAACGCACAATAAGTAATAATAAAAAGACGAACaa  >  1:811632/1‑141 (MQ=255)

CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACAA  >  USA300TCH1516_ALE/2028134‑2028404

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