BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,957	G→T	A194S (GCG→TCG)	yhhI →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,957	0	G	T	97.4%	87.9 / NA	40	A194S (GCG→TCG)	yhhI	putative transposase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (20/18); total (20/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATT  >  NC_000913/3624841‑3625094

tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACAGGGGCGGCta                                                                                                                 >  2:311056/1‑143 (MQ=17)

             tAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAACAGGAACCCAGGGGCGGCTAAATAAAGCCTTTg                                                                                                    >  1:253762/1‑143 (MQ=17)

                    cTGGATATTAAAGGAAAAATCATCACAACTGATGCGCTGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaa                                                                                             >  2:355706/1‑143 (MQ=17)

                       gATAATAAAGGAAAAGTCATCAAAACAGATGAGAAGGATTGCGAAAAAGTTATTGCAGAGAAAAAACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAtt                                                                                          <  2:338317/143‑1 (MQ=11)

                            tAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAACAGGAACCCAGGGGCGGCCAAATCAAGCCTTTGAGGAAAAATTTCCGc                                                                                     >  1:930478/1‑143 (MQ=255)

                            tAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGc                                                                                     >  2:427775/1‑143 (MQ=17)

                                ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaa                                                                                 >  1:499912/1‑143 (MQ=17)

                                 gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                >  1:15401/1‑143 (MQ=17)

                                 gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                <  1:243903/143‑1 (MQ=17)

                                 gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACCGGGCCACCCAAAGAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                >  2:208692/1‑143 (MQ=11)

                                   aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGa                                                                              <  1:668545/143‑1 (MQ=17)

                                      atcctaacaacaGATGCGATAGGATGAGAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAtt                                                                           <  2:90573/137‑1 (MQ=11)

                                        catcaCACCTGATGCGATGGGATGTCAGAATAATATTGCAGAAAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGTAACCCAGTGGCGGCTAAATAAAGTCTTTGAGGAAAAATTTCCGCTGAAAGAATTaa                                                                         <  2:141974/143‑1 (MQ=11)

                                          tcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaat                                                                       >  1:241262/1‑143 (MQ=17)

                                             caACTGAAGCGAGGGGGTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataat                                                                    <  1:335290/143‑1 (MQ=255)

                                               aCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                                                                  >  1:707668/1‑143 (MQ=17)

                                                  gATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCaga                                                               >  2:421944/1‑143 (MQ=18)

                                                   aagCGAGGGGATGCCAGAAAGATATGGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCagag                                                              <  2:730156/141‑1 (MQ=255)

                                                      cGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGTGGAAAAATTTCCGCTTAAAGAATTAAATAATCCAGAGCAt                                                           <  1:249868/143‑1 (MQ=255)

                                                      cGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCAt                                                           >  1:88632/1‑143 (MQ=25)

                                                          gggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACa                                                       <  1:100798/143‑1 (MQ=25)

                                                                 cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAGATAATCCAGAGCATGACAGTTACGc                                                >  1:371754/1‑143 (MQ=25)

                                                                  aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCa                                               >  2:470830/1‑143 (MQ=25)

                                                                      aGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCa                                           <  1:355706/143‑1 (MQ=38)

                                                                         tatTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAACTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTg                                        <  1:811876/143‑1 (MQ=255)

                                                                           ttGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGaa                                      <  1:848917/143‑1 (MQ=255)

                                                                            tGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAATATGATCCTAGGTGCGTCTATATTAAGGCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGaaa                                     <  1:831094/143‑1 (MQ=11)

                                                                              ttgagaAGATACAACAAGAGGGAGGTGATGCTTTATGCGAGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCTCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAag                                   <  1:673606/141‑1 (MQ=11)

                                                                                   aagaTACAATAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCCGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              >  2:643948/1‑143 (MQ=255)

                                                                                   aagaTACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              <  1:606434/143‑1 (MQ=38)

                                                                                   aagaTACAAAAACAGGGAGGCGACTATTTATTCTCGGTAAAAGGAACGCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAACTTTCCGCTGACAGAATGAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              <  1:482204/143‑1 (MQ=11)

                                                                                          aaaaaCAGGGAGGTGATTATTTATTATCGGCAAAAGGGACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCaga                       >  1:697488/1‑143 (MQ=255)

                                                                                           aaaaCAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCagag                      <  2:460151/143‑1 (MQ=38)

                                                                                            aaaCAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAgaga                     >  2:628699/1‑143 (MQ=255)

                                                                                                   gAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTGACGGCAGAGAAGaaacc              >  2:573811/1‑141 (MQ=255)

                                                                                                      gTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGGTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGt           <  1:531879/143‑1 (MQ=255)

                                                                                                       tGATAATTTATTTACGACAAAAGAAACCCAAGGGCAGCTAAATAAAGTCTTTGAGGAAAAATTTCCGATGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTc          <  2:462598/143‑1 (MQ=11)

                                                                                                        gATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGGAGTATATAATCCAGAGCATGACAGGTACGCAATCAGTGAAAAGAGTCACGGCAGAGGAGAAATCGGtgt         >  1:38168/1‑141 (MQ=11)

                                                                                                              ttattCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCatat   >  2:286921/1‑143 (MQ=255)

                                                                                                               tattCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATAtt  <  1:39561/143‑1 (MQ=38)

TCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATT  >  NC_000913/3624841‑3625094

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