BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,911,319	C→A	Q69H (CAG→CAT)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,911,319	0	C	A	100.0%	124.0 / NA	36	Q69H (CAG→CAT)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (16/20); total (16/20)

CAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTG  >  NC_000913/3911187‑3911452

cAGCCCTGAAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTa                                                                                                                                <  2:565240/140‑1 (MQ=255)

  gCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACg                                                                                                                              <  1:432855/140‑1 (MQ=255)

        gAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTa                                                                                                                                >  2:80996/1‑132 (MQ=255)

        gAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTa                                                                                                                                <  1:80996/132‑1 (MQ=255)

                aCCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAAcc                                                                                                                <  2:187257/140‑1 (MQ=255)

                  cGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGa                                                                                                              <  2:718477/140‑1 (MQ=255)

                  cGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGa                                                                                                              >  2:665491/1‑140 (MQ=255)

                    ccAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATa                                                                                                            <  2:487779/140‑1 (MQ=255)

                    ccAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATa                                                                                                            <  1:56500/140‑1 (MQ=255)

                        caccacGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCt                                                                                                        <  2:78554/140‑1 (MQ=255)

                          ccacGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTgg                                                                                                      <  2:175278/140‑1 (MQ=255)

                             cgccgcCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAg                                                                                                   >  1:593305/1‑140 (MQ=255)

                                  ccAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAAc                                                                                              >  1:687531/1‑140 (MQ=255)

                                          ggTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTAc                                                                                      >  2:42950/1‑140 (MQ=255)

                                          ggTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTAc                                                                                      <  1:688537/140‑1 (MQ=255)

                                          ggTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTAc                                                                                      <  1:634262/140‑1 (MQ=255)

                                            tCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCg                                                                                    >  1:273071/1‑140 (MQ=255)

                                                      gAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTttctttct                                                                          >  2:293844/1‑140 (MQ=255)

                                                              cTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAgtg                                                                  >  2:466575/1‑140 (MQ=255)

                                                                 ccTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCa                                                               >  2:244412/1‑140 (MQ=255)

                                                                    gAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGcc                                                            <  1:645900/140‑1 (MQ=255)

                                                                       ccAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCAt                                                         <  1:571040/140‑1 (MQ=255)

                                                                           tttttCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTgg                                                     >  2:715124/1‑140 (MQ=255)

                                                                           tttttCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGACAGCCCATTTgg                                                     <  2:601689/140‑1 (MQ=255)

                                                                                  tCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATacac                                              <  1:674467/140‑1 (MQ=255)

                                                                                        aGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGgcgc                                        <  1:395883/140‑1 (MQ=255)

                                                                                          cggcgCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAg                                      >  2:635814/1‑140 (MQ=255)

                                                                                            gcgcGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAggg                                    <  2:123822/140‑1 (MQ=255)

                                                                                            gcgcGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAggg                                    <  1:399686/140‑1 (MQ=255)

                                                                                               cgTCAGAGGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAg                                 >  2:31797/1‑140 (MQ=255)

                                                                                                    gagGCACCAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgc                            >  1:544197/1‑140 (MQ=255)

                                                                                                          ccAAAATCAACGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTtgcacctgc                      >  2:405809/1‑140 (MQ=255)

                                                                                                                    cGGTATTAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAgg            >  1:452857/1‑140 (MQ=255)

                                                                                                                          tAGCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGgcttgct      >  1:126406/1‑140 (MQ=255)

                                                                                                                            gCGATAATATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTc    <  2:273071/140‑1 (MQ=255)

                                                                                                                                 aaTATGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTg  <  1:218810/137‑1 (MQ=255)

CAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAACCGATACCCTGGTAGTTAACTTTATTACCGGTTTCTTTCTGGTAAGTGTCAGCCCATTTGGCATACACCGGCGCAGGGAAGGTTGCACCTGCACCTGTCAGGCTTGCTTCTG  >  NC_000913/3911187‑3911452

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