BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	380,797	G→A	10.3%	P16S (CCC→TCC)	yaiO ←	outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,797	0	G	A	10.3%	102.0 / 4.5	39	P16S (CCC→TCC)	yaiO	outer membrane protein
Reads supporting (aligned to +/- strand): ref base G (15/20); new base A (2/2); total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

CCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGG  >  NC_000913/380656‑380939

ccTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATgcg                                                                                                                                         <  3:571768/149‑1 (MQ=255)

 cTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATgcgc                                                                                                                                        >  4:89465/1‑149 (MQ=255)

 cTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATgcgc                                                                                                                                        >  1:512853/1‑149 (MQ=255)

   gggAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTa                                                                                                                                      >  1:248211/1‑149 (MQ=255)

        aGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAAt                                                                                                                                 <  1:302569/149‑1 (MQ=255)

        aGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAAt                                                                                                                                 <  2:512853/149‑1 (MQ=255)

           ttAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAAGAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGgc                                                                                                                              >  4:255462/1‑149 (MQ=255)

             aaaaaGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTCCACAATACGCT‑AATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGgccgc                                                                                                                           <  4:578612/149‑1 (MQ=255)

                  gcagcGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCt                                                                                                                       >  4:209352/1‑149 (MQ=255)

                                  aaCTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTa                                                                                                       <  2:190634/149‑1 (MQ=255)

                                  aaCTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTa                                                                                                       >  2:610449/1‑149 (MQ=255)

                                    cTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATc                                                                                                     >  1:501211/1‑149 (MQ=255)

                                     tttCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCa                                                                                                    <  2:641/149‑1 (MQ=255)

                                         gCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCAtttt                                                                                                >  1:233/1‑149 (MQ=255)

                                            accaGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAAt                                                                                             >  4:108045/1‑149 (MQ=255)

                                                         aCGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGAt                                                                                <  3:108045/149‑1 (MQ=255)

                                                           gCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATaa                                                                              <  3:255462/149‑1 (MQ=255)

                                                               aaTTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTc                                                                          <  2:248211/149‑1 (MQ=255)

                                                               aaTTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTc                                                                          <  4:65784/149‑1 (MQ=255)

                                                                          ccATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTTGGAGTg                                                               >  1:583552/1‑149 (MQ=255)

                                                                          ccATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTg                                                               >  4:139916/1‑149 (MQ=255)

                                                                            aTGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGaa                                                             <  1:593801/149‑1 (MQ=255)

                                                                            aTGATCGCCAGAATAATCGGTAAAATCGTAGACCGCGGTAATGGAAGTTAACCCGTCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGaa                                                             <  1:179394/149‑1 (MQ=255)

                                                                             tGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAAt                                                            >  3:279095/1‑149 (MQ=255)

                                                                                  gCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGtt                                                       <  1:454144/149‑1 (MQ=255)

                                                                                     aGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCg                                                    <  4:519199/149‑1 (MQ=255)

                                                                                        ataatCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGACAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCggta                                                 >  3:570882/1‑149 (MQ=255)

                                                                                        ataatCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGACAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCggta                                                 >  3:589938/1‑149 (MQ=255)

                                                                                                tAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAAtt                                         >  4:211544/1‑149 (MQ=255)

                                                                                                  aaaTCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTAc                                       <  3:461436/149‑1 (MQ=255)

                                                                                                    aTCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGt                                     <  3:89465/149‑1 (MQ=255)

                                                                                                                      ggAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCggg                   <  4:279095/149‑1 (MQ=255)

                                                                                                                       gAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGc                  >  1:498320/1‑149 (MQ=255)

                                                                                                                          gTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCaa               <  1:610449/149‑1 (MQ=255)

                                                                                                                               cccGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGc          <  3:367047/149‑1 (MQ=255)

                                                                                                                                   gCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCgtgt      >  4:447760/1‑149 (MQ=255)

                                                                                                                                       aaGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCgg  <  3:233482/149‑1 (MQ=255)

                                                                                                                                       aaGCGGACAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCgg  <  4:570882/149‑1 (MQ=255)

                                                                                                                                       aaGCGGACAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCgg  <  4:589938/149‑1 (MQ=255)

CCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGG  >  NC_000913/380656‑380939

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