BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,418,824	G→A	11.8%	L99L (CTG→TTG)	PP_4758 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,418,824	0	G	A	11.8%	80.2 / 3.3	34	L99L (CTG→TTG)	PP_4758	MFS transporter
Reads supporting (aligned to +/- strand): ref base G (14/16); new base A (2/2); total (16/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 = 7.54e-01

CAACGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGA  >  NC_002947/5418679‑5418961

cAACGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAACcac                                                                                                                                        >  2:100561/1‑150 (MQ=255)

cAACGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAACcac                                                                                                                                        >  2:100536/1‑150 (MQ=255)

   cGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGc                                                                                                                                     <  2:210026/150‑1 (MQ=255)

   cGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAAAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGc                                                                                                                                     <  1:215312/150‑1 (MQ=255)

      gTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCgg                                                                                                                                  <  2:195064/150‑1 (MQ=255)

       tGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCggg                                                                                                                                 >  1:51224/1‑150 (MQ=255)

                   tCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCtttt                                                                                                                     >  2:223202/1‑150 (MQ=255)

                       gCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTgcag                                                                                                                 >  2:372867/1‑150 (MQ=255)

                                 aagCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAACCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTc                                                                                                       <  1:100536/150‑1 (MQ=255)

                                 aagCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTCCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAACCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTc                                                                                                       <  1:100561/150‑1 (MQ=255)

                                    ccTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGaa                                                                                                    <  1:223202/150‑1 (MQ=255)

                                            ggCCACCCCGAGCAGAACACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTc                                                                                            <  1:151202/150‑1 (MQ=255)

                                                   ccGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAgc                                                                                    <  2:370568/150‑1 (MQ=255)

                                                     gAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAgcgc                                                                                  <  2:334467/150‑1 (MQ=255)

                                                             aaCGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGa                                                                          <  2:51224/150‑1 (MQ=255)

                                                                  agcaCGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATgaag                                                                     <  2:323301/150‑1 (MQ=255)

                                                                    caCGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATgaagaa                                                                   >  2:144269/1‑150 (MQ=255)

                                                                        tAGAACTGGGTTTCGGTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGc                                                               >  2:297298/1‑150 (MQ=255)

                                                                              tGGGTTTCGTTCTGCACAAAGGCCATTGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTCCTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGcc                                                         >  2:55174/1‑149 (MQ=255)

                                                                                     cGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGc                                                  <  1:372867/150‑1 (MQ=255)

                                                                                                      aTGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATga                                                                       <  2:435984/112‑1 (MQ=255)

                                                                                                      aTGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATga                                                                       >  1:435984/1‑112 (MQ=255)

                                                                                                               gCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAgg                        <  2:171885/150‑1 (MQ=255)

                                                                                                                 caccaGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATCGGCCGCAGCGCTGATGCCCAGGtg                      >  1:6711/1‑149 (MQ=255)

                                                                                                                          ccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCa                                       <  1:444982/124‑1 (MQ=255)

                                                                                                                          ccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCa                                       >  2:444982/1‑124 (MQ=255)

                                                                                                                          ccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAg             >  1:270274/1‑150 (MQ=255)

                                                                                                                           cAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGACACTGTG‑GTAACCGATGAAAAAACACACATCTGAACTCCCGTAAGGCGCAGCGCTGATGCCCAGGTCGCTTTCCAgg            >  1:211393/1‑150 (MQ=255)

                                                                                                                              gTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgc         <  1:249888/150‑1 (MQ=255)

                                                                                                                                 aaCATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATgaag                                                                     >  1:418634/1‑87 (MQ=255)

                                                                                                                                 aaCATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATgaag                                                                     <  2:418634/87‑1 (MQ=255)

                                                                                                                                   cATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg    >  2:32354/1‑150 (MQ=255)

                                                                                                                                    aTGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcgg   >  2:364255/1‑150 (MQ=255)

                                                                                                                                     tGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGa  <  1:144269/150‑1 (MQ=255)

CAACGGGTGAAATAGTAGATCACGCCAGGAAAGAAGCCTGCTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGA‑ACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGA  >  NC_002947/5418679‑5418961

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