BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,828	A→G	Y42H (TAC→CAC)	pepA ←	aminopeptidase A/I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,828	0	A	G	100.0%	124.4 / NA	35	Y42H (TAC→CAC)	pepA	aminopeptidase A/I
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (19/16); total (19/16)

CGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACG  >  NC_000913/4485692‑4485959

cGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGccc                                                                                                                                  >  2:72198/1‑140 (MQ=255)

    cATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCg                                                                                                                              >  2:206733/1‑140 (MQ=255)

           cTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGAttt                                                                                                                       >  1:333643/1‑140 (MQ=255)

           cTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGAttt                                                                                                                       >  1:374388/1‑140 (MQ=255)

               cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATc                                                                                                                   >  1:3707/1‑140 (MQ=255)

                gTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCg                                                                                                                  >  2:21219/1‑140 (MQ=255)

                gTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCg                                                                                                                  >  1:17716/1‑140 (MQ=255)

                  tCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAg                                                                                                                <  1:21219/140‑1 (MQ=255)

                  tCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAg                                                                                                                <  2:336002/140‑1 (MQ=255)

                  tCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAg                                                                                                                <  1:235466/140‑1 (MQ=255)

                      tGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGt                                                                                                            >  1:339921/1‑140 (MQ=255)

                                        gaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                                          >  1:262737/1‑140 (MQ=255)

                                         aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                         >  2:388153/1‑140 (MQ=255)

                                         aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                         >  2:311430/1‑140 (MQ=255)

                                          aTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                        >  1:4883/1‑140 (MQ=255)

                                           ttCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                             <  1:376829/134‑1 (MQ=255)

                                           ttCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                             >  2:376829/1‑134 (MQ=255)

                                                       gTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAg                                                                           <  2:313823/140‑1 (MQ=255)

                                                            ttCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGcc                                                                      <  1:311430/140‑1 (MQ=255)

                                                                 aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacga                                                                 >  1:210796/1‑140 (MQ=255)

                                                                  aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                                <  2:254261/140‑1 (MQ=255)

                                                                  aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                                <  2:262737/140‑1 (MQ=255)

                                                                       gTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGc                                                           >  2:199567/1‑140 (MQ=255)

                                                                               acaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTc                                                   <  1:388153/140‑1 (MQ=255)

                                                                                  aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGc                                                >  1:284645/1‑140 (MQ=255)

                                                                                  aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGc                                                >  1:34953/1‑140 (MQ=255)

                                                                                  aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGc                                                >  2:45762/1‑140 (MQ=255)

                                                                                          cccGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                                        <  2:4883/140‑1 (MQ=255)

                                                                                                   ccttccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAcc                               <  2:210796/140‑1 (MQ=255)

                                                                                                               ccccGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                   >  1:276440/1‑140 (MQ=255)

                                                                                                               ccccGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                   <  2:3707/140‑1 (MQ=255)

                                                                                                                         gCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg         <  1:206733/140‑1 (MQ=255)

                                                                                                                              gcgcTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTa    <  1:234864/140‑1 (MQ=255)

                                                                                                                               cgcTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt              <  2:340060/129‑1 (MQ=255)

                                                                                                                               cgcTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt              >  1:340060/1‑129 (MQ=255)

                                                                                                                                gcTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACg  <  2:68310/140‑1 (MQ=255)

CGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACG  >  NC_000913/4485692‑4485959

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