BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,446,968	A→T	9.5%	N89K (AAT→AAA)	yojO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,446,968	0	A	T	9.5%	120.0 / 4.5	42	N89K (AAT→AAA)	yojO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (22/16); new base T (2/2); total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TCTAAGCGAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAG  >  CP000730/1446831‑1447098

tCTAAGCGAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCAttt                                                                                                                                 >  1:344573/1‑141 (MQ=255)

tCTAAGCGAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCAttt                                                                                                                                 <  1:174282/141‑1 (MQ=255)

    aGCGAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACaa                                                                                                                             >  1:430942/1‑141 (MQ=255)

       gAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAAttt                                                                                                                          >  1:275452/1‑141 (MQ=255)

          gATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTa                                                                                                                       >  2:295821/1‑141 (MQ=255)

            tCAATTAATTTAGCAGTACTTTGGCGCTGAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTAca                                                                                                                     <  2:339693/141‑1 (MQ=255)

            tCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTAca                                                                                                                     >  1:314806/1‑141 (MQ=255)

              aattaattTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTAcatc                                                                                                                   <  2:357410/141‑1 (MQ=255)

                ttaattTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTAcatcca                                                                                                                 >  1:401581/1‑141 (MQ=255)

                    tttAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCcatcat                                                                                                             <  1:155737/141‑1 (MQ=255)

                         cAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTATTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGAt                                                                                                        >  2:23707/1‑141 (MQ=255)

                                tGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTa                                                                                                 <  1:299735/141‑1 (MQ=255)

                                    cgcTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGGAagca                                                                                             >  2:369962/1‑141 (MQ=255)

                                     gcTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAagcag                                                                                            <  2:55151/141‑1 (MQ=255)

                                        cAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCta                                                                                         >  2:450177/1‑141 (MQ=255)

                                          aCTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACTATTTGTACATCCATCATATGATAACCGTAAGCAGCtaac                                                                                       >  1:231013/1‑141 (MQ=255)

                                          aCTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCtaac                                                                                       >  1:231009/1‑141 (MQ=255)

                                               aaTAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTACAAGATTTATATAGAACTAACTCAATACCAATTTGTACATCCATAATATGATAACCGTAAGCAGCtaacataac                                                                                  <  1:65258/141‑1 (MQ=255)

                                               aaTAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCtaacataac                                                                                  <  1:136267/141‑1 (MQ=255)

                                                aTAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCtaacataaca                                                                                 >  2:316812/1‑141 (MQ=255)

                                                          aTACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTa                                                                       >  1:235200/1‑141 (MQ=255)

                                                                         ctAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATATTATTTGTTTTTAAACAAGCCATTTCGa                                                        >  2:362043/1‑141 (MQ=255)

                                                                          taataaTTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATATCGaa                                                       <  2:141000/141‑2 (MQ=255)

                                                                               atTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGc                                                  <  1:170688/141‑1 (MQ=255)

                                                                                   aaaaaGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCTTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTgg                                              >  2:442843/1‑141 (MQ=255)

                                                                                   aaaaaGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCTTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTgg                                              >  2:443104/1‑141 (MQ=255)

                                                                                    aaaaGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGt                                             >  2:445098/1‑141 (MQ=255)

                                                                                        gttgttGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAt                                         >  1:383703/1‑141 (MQ=255)

                                                                                           gttgATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGAtgtg                                      >  1:267123/1‑141 (MQ=255)

                                                                                              gATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGAcc                                   >  1:85594/1‑141 (MQ=255)

                                                                                                tAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCaa                                 <  2:299856/141‑1 (MQ=255)

                                                                                                    aTTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAAcc                             >  2:180405/1‑141 (MQ=255)

                                                                                                     tttAGGGTATTTAAAAGTTTTATCTTGGACTACTGCATTTGCAAGTTGTACATCCATCATATGATAACCGTAAGCTGATAACATAACAGCTGTTTTTAAACCACCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCa                            <  2:240067/141‑1 (MQ=255)

                                                                                                       tAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATGAGAGGATGACCGAATTCATCTTACATGACATCGGTTTTTAAGCTAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAAc                          >  1:154838/1‑141 (MQ=255)

                                                                                                        aGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACt                         >  1:315773/1‑141 (MQ=255)

                                                                                                                tAAAAGTTTTATCTTGTACTACTTCTTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                 <  1:443104/141‑1 (MQ=255)

                                                                                                                tAAAAGTTTTATCTTGTACTACTTCTTTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                 <  1:442843/141‑1 (MQ=255)

                                                                                                                tAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCa                 <  2:235200/141‑1 (MQ=255)

                                                                                                                 aaaaGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCaa                <  1:14073/141‑1 (MQ=255)

                                                                                                                     gTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACa            <  1:316812/141‑1 (MQ=255)

                                                                                                                        ttATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACAttt         <  2:20639/141‑1 (MQ=255)

                                                                                                                               gTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAg  >  1:176534/1‑141 (MQ=255)

TCTAAGCGAAGATCAATTAATTTAGCAGTACTTGGGCGCTCAACTTTAATAGAATTTAATACGCGCATATCTTCTAATAATTTAAAAAGTTGTTGATAAAATTTAGGGTGTTTAAAAGTTTTATCTTGTACTACTTCATTTACAATTTGTACATCCATCATATGATAACCGTAAGCAGCTAACATAACATCTGTTTTTAAACCAGCCATTTCGATATGGCTTGGTCGATGTGACCAAAACCAACTCGTAATCAAAACATTTTGAACAG  >  CP000730/1446831‑1447098

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