BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	445,241	A→G	9.7%	L207P (CTT→CCT)	qorB ←	Quinone oxidoreductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	445,241	0	A	G	9.7%	109.1 / 3.1	41	L207P (CTT→CCT)	qorB	Quinone oxidoreductase 2
Reads supporting (aligned to +/- strand): ref base A (16/21); new base G (2/2); total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

GTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGA  >  USA300TCH1516_ALE/445101‑445375

gTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCaa                                                                                                                                        <  2:240723/141‑1 (MQ=255)

      ttGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCtt                                                                                                                                  >  1:98412/1‑141 (MQ=255)

         cTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTtca                                                                                                                               >  1:331355/1‑141 (MQ=255)

                     aCATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCGGTTAAAATTGCAGCCAGTTCTTTCATATCATAACTGt                                                                                                                   >  1:24729/1‑141 (MQ=255)

                       aTTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAg                                                                                                                 >  1:145917/1‑141 (MQ=255)

                          gATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATACCTGTAGCCt                                                                                                              >  2:270848/1‑141 (MQ=255)

                                   acggtaCCAAAGCCTTGAGGTTCATAATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAg                                                                                                     <  2:299333/136‑1 (MQ=255)

                                    atGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGt                                                                                                    <  1:255142/141‑1 (MQ=255)

                                    atGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGt                                                                                                    <  1:255195/141‑1 (MQ=255)

                                        aCCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTTCAGGAAGTTCTTTCATATCATTACTGTAGCCTGATAATAAGTAgag                                                                                                >  2:94767/1‑139 (MQ=255)

                                         ccAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGt                                                                                               <  2:58506/141‑1 (MQ=255)

                                                   tAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGTCTGATAATAAGGAGCGTTTGCCCCAAg                                                                                     <  2:62647/141‑1 (MQ=255)

                                                   tAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAAGTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAg                                                                                     <  2:62583/141‑1 (MQ=255)

                                                   tAGGTTCATCATACATTTCGTCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAg                                                                                     >  2:187430/1‑141 (MQ=255)

                                                        tcatcaTACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCt                                                                                <  1:56510/141‑1 (MQ=255)

                                                         catcatACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTg                                                                               >  1:26442/1‑141 (MQ=255)

                                                             atACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGAtt                                                                           <  2:317017/141‑1 (MQ=255)

                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  1:10444/1‑141 (MQ=255)

                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  1:10434/1‑141 (MQ=255)

                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  2:122396/1‑141 (MQ=255)

                                                                              tctcTAATGAAACGGGCTCATATTTAATCTCTGTGCCGGATGCATAAAGGGATATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGa                                                          <  1:269761/141‑1 (MQ=255)

                                                                               ctctAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGac                                                         >  1:361096/1‑141 (MQ=255)

                                                                               ctctAATGAAACGGGCTCAGATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTACTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGac                                                         >  1:278319/1‑141 (MQ=255)

                                                                                         aCGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCAAAAGTGTCTGGATTTTTAATAATAGCAATGACACCTCTAGCa                                               >  1:68563/1‑141 (MQ=255)

                                                                                          cGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCaa                                              >  1:362148/1‑141 (MQ=255)

                                                                                          cGGGCTCATATTTAATTCCCATGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCaa                                              >  2:361436/1‑141 (MQ=255)

                                                                                             gCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAAAGACACCTCTAGCAatat                                           >  1:259817/1‑141 (MQ=255)

                                                                                               tCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGTTTTTTAATAATAGCAATGACACCTCTAGCAATATCa                                         <  2:186474/141‑1 (MQ=255)

                                                                                                  tatTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCAttt                                      <  2:336837/141‑1 (MQ=255)

                                                                                                         tttCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTaa                               <  2:230137/141‑1 (MQ=255)

                                                                                                           tCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAata                             <  2:10444/141‑1 (MQ=255)

                                                                                                           tCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAata                             <  2:10434/141‑1 (MQ=255)

                                                                                                            cTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAatat                            <  1:343647/141‑1 (MQ=255)

                                                                                                               tgCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATAtaat                         <  1:187430/141‑1 (MQ=255)

                                                                                                                    gATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCCTAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATa                    <  2:280837/141‑1 (MQ=255)

                                                                                                                      tGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACg                  <  1:140077/141‑1 (MQ=255)

                                                                                                                             gATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCg           <  2:78188/141‑1 (MQ=255)

                                                                                                                              aTAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGc          <  1:122396/141‑1 (MQ=255)

                                                                                                                                   actGCATCAAGTTCTTGCATATCATAACTGTAGCCTGATCATAGTTAGCGTGTGCCTCAAGTATCTGCATTTTTAATAATAGCAATGTCACCTCTAGCTTTATCATTTCTAGTAGTATGATTAATACGACCATTGCCAggt     <  2:141141/139‑3 (MQ=255)

                                                                                                                                    ttGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTg    <  1:287259/141‑1 (MQ=255)

                                                                                                                                     tGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTgg   <  1:360592/141‑1 (MQ=255)

                                                                                                                                      gcagcaAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGGTTGCCCCAAGTATCTGGATTTTTAATCATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGa  >  1:312154/1‑141 (MQ=255)

GTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGA  >  USA300TCH1516_ALE/445101‑445375

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