BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	459,694	G→A	E269E (GAG→GAA)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	459,694	0	G	A	96.0%	70.5 / NA	25	E269E (GAG→GAA)	lon	Lon protease
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (10/14); minor base C (1/0); total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTCGCAACCGCGTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATC  >  NC_000913/459522‑459899

gCATTCGCAACCGCGTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATc                                                                                                                                 <  1:115510/251‑1 (MQ=255)

    tCGCAACCGCGTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGGAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGan                                                                                                                              <  1:302189/250‑2 (MQ=255)

       cAACCGCGTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg                                                                                                                           >  1:290997/1‑250 (MQ=255)

              gTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAACGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGCAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCa                                                                                                                                                                                           >  1:263878/1‑179 (MQ=255)

                tAAAAAGCAGATGGATAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGGCCTGAAGCTCAAAAGCGACTCGGTGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAgn                                                                                                                                                                                      <  1:94221/182‑2 (MQ=255)

                         naTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGCACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCg                                                                                                            >  1:405436/2‑247 (MQ=255)

                                        gCGTGAGTACTATCTGAACGAGAAAATGAAAGCTATTCAGAAAGAACGCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCTACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAgg                                                                                         <  1:71131/251‑1 (MQ=255)

                                                       gAACGAGAAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAAGCGACGCGGCGAAGACGCCGAAAGAGGCAAAAGAAAAAGCGGAGGCAGAGTTGCAGAAGGCGAAAATGATGTCTCCGATGTCGGCAGAAGCTGCCGGAGTGTGTGTTTATACCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGAcc                                                                             <  1:293159/248‑1 (MQ=255)

                                                               aaaTGAAAGCTATGCAGAAAGACCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGGCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCg                                                                         <  1:459/244‑1 (MQ=255)

                                                                      aGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAt                                                                                                                                                                   <  1:461718/147‑1 (MQ=255)

                                                                             cAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGCAGCGCAAAATTGACGCGGCGAAGATGCCGAAAGAGGCATAAGCAAACGAGGAAGCAGAGTTGCAGAAGCTGAAAatg                                                                                                                                                                              >  1:357908/1‑129 (MQ=255)

                                                                                 tccgagctggTGAAAAGGACGCCGCGCTGGAAGAAAACGAAGCACTGAAGCGCAAAATCGACGCGGCGAAGATGCAGAAAGAGTCAAAAGAAAAAGCGGGAGCAGAGTTGCAGAAGCTGAAAATGATGTCGTCGATGGCGGCAGAAGCGTCTGTAGTGCGTTGTTAGATCGACTGGATGGTTCAGTGGCCGGGGGATTCGCGTAGCAAGATTAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAc                                                 <  1:352376/242‑1 (MQ=255)

                                                                                             aaaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGCAAGAGGCAAACGAAAACGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCACTGTCGGCAGAAGCGACCGTAGTGCGCGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGCTACCGACCATTATGGt                                    >  1:241982/1‑251 (MQ=255)

                                                                                                              gACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCTGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCa                                                                                            <  1:33088/178‑1 (MQ=255)

                                                                                                                          gCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAACGCGGCAAAAGACAAAGCGGTTGCAGAGTTGCAGAAGATGAAAATGATGTCTCCGATGCCGGCAGAAGAGACCGTAGTGCGTGGT‑ATATCGACTGGATGGTACAGGTGACGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTgg                                >  1:344478/1‑225 (MQ=255)

                                                                                                                           cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAACGCGGAAGCAGAGTTGCAGAAGCTGAGAATGAAGTCTACGATGTCGGCAGaa                                                                                                                                                         >  1:481483/1‑104 (MQ=255)

                                                                                                                            ccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGNTGGNNNNngg                                                                                                                >  2:311164/1‑143 (MQ=255)

                                                                                                                            ccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGNTGGNNNNNGtnnnnnn                                                                                                          >  2:308839/1‑144 (MQ=255)

                                                                                                                            ccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCgn                                                                                                           <  1:311164/149‑2 (MQ=255)

                                                                                                                             cTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGGTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGAGGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGTCCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCAGTATGGTCTGGAGCGCGTGaa                      <  1:269620/233‑1 (MQ=255)

                                                                                                                               gAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGCGTTGCAGAAGCTGAAAATGCTGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACCGGTGCCGTGGAATGCGCGTAGCAAGGTCCAAAAAGACCTGCGTCCGGCGCAGGCAATCCTTGATACCGACCCTTATGGTCTGGAGCGCGTGAAAGCTCGAATCCTTGAGTATc  >  1:330507/1‑251 (MQ=255)

                                                                                                                                 aGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAAAAAGCGGAAGCAGAGTTGCAGGAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCTTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATgg                                     <  1:498180/214‑1 (MQ=255)

                                                                                                                                                         tGCCGAAAGAGGAAAAAGAAAAAGCGGACGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATcc                                                        >  1:410443/1‑171 (MQ=255)

                                                                                                                                                                    gCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTTGTTATATCGACTGGATng                                                                                                                       <  1:394065/97‑1 (MQ=255)

                                                                                                                                                                    gCAAAAGAAAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGAc                                             <  1:210674/171‑1 (MQ=255)

GCATTCGCAACCGCGTTAAAAAGCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATC  >  NC_000913/459522‑459899

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