BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,125,097	C→T	E572K (GAG→AAG)	priA ←	primosome factor N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,125,097	0	C	T	95.5%	123.0 / NA	45	E572K (GAG→AAG)	priA	primosome factor N'
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (23/19); minor base G (1/1); total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124957‑4125235

ggCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACAACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTg                                                                                                                                          <  1:496192/143‑1 (MQ=255)

     aGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAg                                                                                                                                     <  2:643281/143‑1 (MQ=255)

     aGGCGCTGTGGCCTTTATAGAGCAACGCTGGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTCCCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAg                                                                                                                                     <  2:794983/143‑1 (MQ=255)

        cgtcgtAGCGTTTATGTCGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTGCGCCCTGTTTACCCGCAGGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGa                                                                                                                                  <  2:208626/143‑1 (MQ=255)

          tcgtAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGaag                                                                                                                                >  1:773187/1‑142 (MQ=255)

             tAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATCGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGCCCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCCAAAATc                                                                                                                             >  1:401460/1‑143 (MQ=255)

                ccTTTATAGAGCAACGTTTGCAACCGACGATGTTCCGGATGGTGCGTATGCAGCACCACTTCGCCCCCTTTACGCGCACGCCCGTCACGACACGCGCCACGGGTGTAAAGCTGCGCACAGCGCTCGGCCGCTCGaaccgcggg                                                                                                                          >  2:351151/1‑136 (MQ=255)

                 cTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTCAAGCTGAGCGCAACGCTTTGCCGAGCGAAACTCCGcc                                                                                                                         >  2:439632/1‑142 (MQ=255)

                  tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCGTGTTGACCCGCACGCCGGGCGCGACCGGCCACCTGGGGGTAAAGCGGAGGGAAACGCTGCGCCGAGCAAGAATGCGCAg                                                                                                                        >  2:731777/1‑143 (MQ=255)

                  tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGGCCTGGGTGTAAAGCTGAGCGAAACGCATTGCCGAGCGAAAATCGGccg                                                                                                                        >  1:291274/1‑141 (MQ=255)

                     ataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaaa                                                                                                                     >  1:223830/1‑143 (MQ=255)

                     ataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaaa                                                                                                                     >  1:841452/1‑143 (MQ=255)

                      taGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGaaaa                                                                                                                    >  2:361713/1‑143 (MQ=255)

                       agagCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGACAATCGGCAGCAAAc                                                                                                                   >  2:540874/1‑143 (MQ=255)

                            aaCGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGCGTAAAGCTGAGCGCAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAgcgc                                                                                                              >  2:482808/1‑143 (MQ=255)

                                     aaCAGAGTATGGTCCGGATGGGGCGTTTGCAGCACCACTTCGCCCGGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccac                                                                                                     <  1:751860/143‑1 (MQ=255)

                                      aCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGCTTACCCGCCCGCCCGGCACGACCGCCGCCCCGGGTGTACCGCTGGGCGAGACGCCTCGCCGCGCGCAACTCGGCAGACACCACCCCGCcccccacg                                                                                                    >  2:102873/1‑143 (MQ=255)

                                           ggATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAg                                                                                               <  1:122485/143‑1 (MQ=255)

                                                 tCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAGCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGc                                                                                         >  2:24085/1‑143 (MQ=255)

                                                   cGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGCAAATCGGCCGAAAACAGCGCGCCGGCCACGTCCCGTAATGCaa                                                                                       >  2:233076/1‑143 (MQ=255)

                                                   cGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCaa                                                                                       <  1:361713/143‑1 (MQ=255)

                                                   cGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAAGCGGCAGAAAACAGCGCGCCGTCCACGTCCGGTAGTGcga                                                                                       >  2:210284/1‑141 (MQ=255)

                                                           gCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTc                                                                               <  2:86462/143‑1 (MQ=255)

                                                                  cagcaCCACTTCGCCCGGTTTACCCGCACGCCCGGCACGACCCCCGACCTGCGTGTAAAGCTGCGCGAAACGCTTTGCCGACCGCAAATCGGCCGAAAACACCCCGCCGGACCCATCACGTAATCCCACCAGCGTccgacccc                                                                        >  2:238704/1‑136 (MQ=255)

                                                                      accacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGaa                                                                    <  2:829379/143‑1 (MQ=255)

                                                                          cTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATgg                                                                <  1:592623/143‑1 (MQ=255)

                                                                              gCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAGTGCAACCAGCGTCACATCCGGGAAATGGTGAc                                                            >  2:25925/1‑143 (MQ=255)

                                                                              gCCCGGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCGGCGTCACATCCGGGAAATGGTGAc                                                            <  2:436212/143‑1 (MQ=255)

                                                                               cgcTGGTTACCCGCACGGGGGGCACGACCGGCGAGTTGGGTGTAAGGGTGAGCGAAGCGCTGTGTGGAGCGAAAATCGGCAGAAAACAGCGCGCCGTGCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAcc                                                           <  2:12653/141‑1 (MQ=255)

                                                                                ccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCt                                                          >  2:846671/1‑143 (MQ=255)

                                                                                      tACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGAAGAAAACAGCGCGCCGTCCACGTGCCGTAATGCAACCAGCGTGACATCCGGGAAATGGGGACCtctgggc                                                    >  2:848322/1‑137 (MQ=255)

                                                                                       aCCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCa                                                   <  2:661174/143‑1 (MQ=255)

                                                                                         ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGc                                                 <  1:434308/143‑1 (MQ=255)

                                                                                          cGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCa                                                >  1:408723/1‑143 (MQ=255)

                                                                                              cgtccGGCTCGCCCGGCGGCCTGTGTGTAAAGCTGAGCGGAACGCTTTGCCGAGCGAAACTCGGCGGAAAACAGCGCGCCGTCCACGTCGAGTAATGCTACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtg                                            <  2:75174/140‑1 (MQ=255)

                                                                                                cccGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtgtg                                          <  1:600991/143‑1 (MQ=255)

                                                                                                      aCGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGGAccga                                    >  1:432156/1‑141 (MQ=255)

                                                                                                                aCCTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATcgg                          >  1:408716/1‑141 (MQ=255)

                                                                                                                 ccTGGGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAGTCAAAATCcgc                         >  1:152233/1‑143 (MQ=255)

                                                                                                                      gcgtAACGCTGAGCGAAACGCCTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGcgcc                    <  1:210284/141‑1 (MQ=255)

                                                                                                                             gCTGAGCGATGCGCTTGGCCGGGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACAGCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGa             <  2:95049/143‑1 (MQ=255)

                                                                                                                             gCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGa             <  1:512392/143‑1 (MQ=255)

                                                                                                                              cTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGAt            >  1:307678/1‑143 (MQ=255)

                                                                                                                                   cGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTagt       >  1:531824/1‑141 (MQ=255)

                                                                                                                                        cGCTTTGCCGAGCGCACACCGGCAGCAACCAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGc  <  2:576041/143‑1 (MQ=255)

GGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124957‑4125235

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: `‑`