BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,744,285	C→T	10.2%	H203Y (CAC→TAC)	PP_2398 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,744,285	0	C	T	10.2%	99.8 / 3.5	39	H203Y (CAC→TAC)	PP_2398	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (19/16); new base T (2/2); total (21/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 = 4.92e-01

CTTCATCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAA  >  NC_002947/2744146‑2744427

cTTCATCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCgg                                                                                                                                      <  1:76617/150‑1 (MQ=255)

   cATCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCa                                                                                                                                   <  2:970918/150‑1 (MQ=255)

    aTCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCAt                                                                                                                                  >  1:673233/1‑150 (MQ=255)

     tCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCAtt                                                                                                                                 >  2:277817/1‑150 (MQ=255)

             gTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGaa                                                                                                                         <  2:258045/150‑1 (MQ=255)

                cTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGGTGCAGGTTGCCGGGCATTACGATGAAgca                                                                                                                      >  1:134363/1‑149 (MQ=255)

                    aagaagCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCgg                                                                                                                  <  2:673233/150‑1 (MQ=255)

                                      tgctgctgGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATAc                                                                                                  <  2:534183/148‑1 (MQ=255)

                                      tgctgctgGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATAc                                                                                                  >  1:534183/1‑148 (MQ=255)

                                           ctgGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATggg                                                                                           <  1:244580/150‑1 (MQ=255)

                                                  tCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATg                                                                                             <  1:222012/141‑1 (MQ=255)

                                                  tCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATg                                                                                             >  2:222012/1‑141 (MQ=255)

                                                         tGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCATGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAagga                                                                             >  1:60478/1‑150 (MQ=255)

                                                             tATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAgcgc                                                                                                                     <  1:80631/106‑1 (MQ=255)

                                                             tATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAgcgc                                                                                                                     >  2:80631/1‑106 (MQ=255)

                                                                      gCCTGCAACCATGGCTATGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGAGCGAGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGAATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGgcg                                                                >  2:175522/1‑150 (MQ=255)

                                                                       ccTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATg                                                                                             >  1:563914/1‑120 (MQ=255)

                                                                       ccTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATg                                                                                             <  2:563914/120‑1 (MQ=255)

                                                                       ccTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGgcgc                                                               <  2:735338/150‑1 (MQ=255)

                                                                          gCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGt                                                            >  2:493094/1‑150 (MQ=255)

                                                                            aaCCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTAc                                                          >  1:684507/1‑150 (MQ=255)

                                                                               cATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGcc                                                       >  2:176378/1‑150 (MQ=255)

                                                                                             ccAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGTACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGAt                                                                         <  1:616257/118‑1 (MQ=255)

                                                                                             ccAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGTACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGAt                                                                         >  2:616257/1‑118 (MQ=255)

                                                                                             ccAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTgcgcg                                         >  1:1201584/1‑150 (MQ=255)

                                                                                               aGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGTACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATgggg                                                                                          <  2:130053/99‑1 (MQ=255)

                                                                                               aGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGTACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATgggg                                                                                          >  1:130053/1‑99 (MQ=255)

                                                                                                  aaTTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTg                                    <  1:176378/150‑1 (MQ=255)

                                                                                                      ccTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGt                                <  1:277817/150‑1 (MQ=255)

                                                                                                       cTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTa                               >  1:322009/1‑150 (MQ=255)

                                                                                                            aGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGcc                          <  1:493094/150‑1 (MQ=255)

                                                                                                              gCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGa                        >  1:974707/1‑150 (MQ=255)

                                                                                                              gCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGa                        >  1:1002369/1‑150 (MQ=255)

                                                                                                                ttGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGCTGACCTGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACg                      <  1:921390/150‑1 (MQ=255)

                                                                                                                ttGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGt                                                            <  2:794388/112‑1 (MQ=255)

                                                                                                                ttGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGt                                                            >  1:794388/1‑112 (MQ=255)

                                                                                                                ttGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACg                      >  1:50229/1‑150 (MQ=255)

                                                                                                                           gcgcgTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACg           >  2:1080365/1‑150 (MQ=255)

                                                                                                                                    cGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCaa  <  2:684507/150‑1 (MQ=255)

CTTCATCCAGGCGGTGCTTGAAGAAGCAGATTGTGACCTGCTGCTGGACGTCAACAATGTCTATGTCAACGCCTGCAACCATGGCTACGACGCCCAGCAATTCCTTGCAGGCTTGCCTCAGGCGCGCGTGGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAA  >  NC_002947/2744146‑2744427

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