BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,390,147	G→A	A151V (GCC→GTC)	mpaA ←	murein peptide amidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,390,147	0	G	A	97.9%	156.4 / ‑4.0	48	A151V (GCC→GTC)	mpaA	murein peptide amidase A
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (24/23); total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CAGAGGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTC  >  NC_000913/1390050‑1390241

cAGAGGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcagaca                                                                                              <  2:179940/100‑1 (MQ=255)

    ggCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGtt                                                                                          >  2:71167/1‑100 (MQ=255)

    ggCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGtt                                                                                          >  1:943764/1‑100 (MQ=255)

     gCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTg                                                                                         <  1:1121844/100‑1 (MQ=255)

         ttCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATg                                                                                     <  1:284836/100‑1 (MQ=255)

             aaCGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGgt                                                                                 <  2:339018/100‑1 (MQ=255)

              aCGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACATGTTGAATGCGgtg                                                                                <  2:803669/100‑1 (MQ=255)

                gCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGta                                                                              <  2:47195/100‑1 (MQ=255)

                gCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGCGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGta                                                                              >  2:835678/1‑100 (MQ=255)

                 ccTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGTATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGtat                                                                             >  2:917847/1‑100 (MQ=255)

                 ccTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCATACAGGTTGAATGCGGTGtat                                                                             <  2:796125/100‑1 (MQ=255)

                        cagccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTgg                                                                      >  2:74706/1‑100 (MQ=255)

                         agccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGc                                                                     <  1:508046/100‑1 (MQ=255)

                           ccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGcac                                                                   <  1:436373/100‑1 (MQ=255)

                               tCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAgcg                                                               >  1:413514/1‑100 (MQ=255)

                                  ccTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCt                                                            >  1:869391/1‑100 (MQ=255)

                                   cTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTg                                                           <  2:1278445/100‑1 (MQ=255)

                                        tCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCCCAGAGCCTGGGttt                                                      >  1:912022/1‑100 (MQ=255)

                                             gtgtCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGGGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAgg                                                  <  1:138966/99‑1 (MQ=255)

                                             gtgtCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGGCGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCa                                                    >  2:639661/1‑97 (MQ=255)

                                              tgtCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACACAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGCTTCAGGtt                                                >  2:205172/1‑100 (MQ=255)

                                                  tGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACGGGTTGAATGCGGTGTATAAGCGGGCACAGCGCCTGGGTTTCAGGTGCg                                              >  2:777655/1‑98 (MQ=255)

                                                  tGGGATCTTCAATACAGGCCAGTGGACCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGGTTCAGGTTCGGa                                            >  2:728566/1‑100 (MQ=255)

                                                   gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAg                                           <  2:978546/100‑1 (MQ=255)

                                                      aTCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCa                                        >  1:932569/1‑100 (MQ=255)

                                                      aTCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCa                                        <  2:614327/100‑1 (MQ=255)

                                                      aTCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGGTCGGAGcc                                         >  2:663622/1‑99 (MQ=255)

                                                         ttCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGt                                     >  1:479107/1‑100 (MQ=255)

                                                           cAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGttt                                   >  2:998412/1‑100 (MQ=255)

                                                            aaTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTg                                  <  1:728907/100‑1 (MQ=255)

                                                                cAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc                              >  1:12128/1‑100 (MQ=255)

                                                                cAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGACTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc                              >  2:481753/1‑100 (MQ=255)

                                                                 aGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAcc                             <  1:697837/100‑1 (MQ=255)

                                                                 aGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAcc                             <  2:415247/100‑1 (MQ=255)

                                                                  ggCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCt                            >  2:723016/1‑100 (MQ=255)

                                                                     cAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTc                         <  2:569333/100‑1 (MQ=255)

                                                                       gTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAg                       <  2:978762/100‑1 (MQ=255)

                                                                       gTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAg                       <  1:498609/100‑1 (MQ=255)

                                                                       gTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAg                       <  1:71168/100‑1 (MQ=255)

                                                                         ggATCGTGGAATGAGACCACCCAGACAGGTGGAATGCGGTGTATAAGCAGGCACAGCGCCTGCGGTTAAGGTTCGAAGCCAGGTTTGTCACCTGTCAg                       >  2:227830/1‑98 (MQ=255)

                                                                              ctGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACg                <  1:829567/99‑1 (MQ=255)

                                                                                ggAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGag              >  2:829561/1‑99 (MQ=255)

                                                                                ggAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGa               <  2:1107712/99‑1 (MQ=255)

                                                                                    tGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATcg          >  2:408336/1‑100 (MQ=255)

                                                                                     gagaCCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATcgc         <  2:355867/100‑1 (MQ=255)

                                                                                        accaccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCGTGGGTTTCAGGTTCGGAGCCAGGTTCGGCACCTGTCAGTAATACGACATcgcg        >  2:1153619/1‑98 (MQ=255)

                                                                                          caccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGttct    >  1:465235/1‑100 (MQ=255)

                                                                                            cccAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGttcttc  >  1:544664/1‑100 (MQ=255)

                                                                                            cccAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGttcttc  <  2:136186/100‑1 (MQ=255)

CAGAGGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTC  >  NC_000913/1390050‑1390241

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