BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	630,463	G→T	10.2%	T411T (ACC→ACA)	PP_0544 ←	ethanolamine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	630,463	0	G	T	10.2%	127.2 / 2.7	39	T411T (ACC→ACA)	PP_0544	ethanolamine transporter
Reads supporting (aligned to +/- strand): ref base G (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 = 3.17e-01

GTAGGCAATTAATACCGCATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCTGGCCCGTCAGC  >  NC_002947/630337‑630592

gTAGGCAATTAATACCGCATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAg                                                                                                                          <  2:259671/136‑1 (MQ=255)

    gCAATTAATACCGCATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGgc                                                                                                                      <  2:564986/136‑1 (MQ=255)

              ccGCATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCtt                                                                                                            >  2:597613/1‑136 (MQ=255)

                ccATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTg                                                                                                          <  2:612748/135‑1 (MQ=255)

                   tAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGcc                                                                                                       >  2:37587/1‑136 (MQ=255)

                    aGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcg                                                                                                      >  2:662061/1‑136 (MQ=255)

                    aGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcg                                                                                                      >  2:227043/1‑136 (MQ=255)

                    aGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcg                                                                                                      >  2:251563/1‑136 (MQ=255)

                    aGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcg                                                                                                      <  1:6185/136‑1 (MQ=255)

                    aGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcg                                                                                                      >  1:317050/1‑136 (MQ=255)

                     gATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcgg                                                                                                     <  2:27951/136‑1 (MQ=255)

                      aTCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCcggc                                                                                                    <  1:359496/136‑1 (MQ=255)

                         aCTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGa                                                                                                 >  2:263621/1‑136 (MQ=255)

                         aCTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGa                                                                                                 >  1:402369/1‑136 (MQ=255)

                            gCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATgc                                                                                              <  1:151326/136‑1 (MQ=255)

                              aGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATgcgc                                                                                            >  1:48470/1‑136 (MQ=255)

                                   cAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGt                                                                                       >  1:467960/1‑136 (MQ=255)

                                        acaaCCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGCATGTGAAGATGCCACCCGGGGGGCGCTAGGGGCGCGCCATCTTTGGCCGGCGGATGCGCAGGGAGCtgt                                                                                  >  2:236140/4‑136 (MQ=255)

                                        aCTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGAtgt                                                                                  <  1:618921/136‑1 (MQ=255)

                                             ccGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCa                                                                             >  2:475522/1‑136 (MQ=255)

                                                    tCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCcatca                                                                      <  2:521166/136‑1 (MQ=255)

                                                    tCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCcatca                                                                      <  2:318438/136‑1 (MQ=255)

                                                         cAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGc                                                                 <  1:15537/136‑1 (MQ=255)

                                                          aGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCa                                                                >  1:292556/1‑136 (MQ=255)

                                                              aTCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTa                                                            <  1:547817/136‑1 (MQ=255)

                                                                ccGGCGCCCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCCTCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAg                                                          >  1:9834/1‑136 (MQ=255)

                                                                    cGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAg                                                      <  2:340792/136‑1 (MQ=255)

                                                                       ccACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGt                                                   >  2:260468/1‑136 (MQ=255)

                                                                                     gCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGTGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAATACAGGGTAGCGCCAAACACTg                                     >  2:614862/1‑136 (MQ=255)

                                                                                     gCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGTGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAATACAGGGTAGCGCCAAACACTg                                     >  2:657723/1‑136 (MQ=255)

                                                                                      ccAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGc                                    <  1:187588/136‑1 (MQ=255)

                                                                                        aGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGcca                                  >  2:443211/1‑136 (MQ=255)

                                                                                                   cAGGCCGGATGTGAAGATGCCACCCGGTGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAATACAGGGTAGCGCCAAACACTGCCACCAGGATCAGc                       <  1:614862/136‑1 (MQ=255)

                                                                                                   cAGGCCGGATGTGAAGATGCCACCCGGTGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAATACAGGGTAGCGCCAAACACTGCCACCAGGATCAGc                       <  1:657723/136‑1 (MQ=255)

                                                                                                             gtgAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCt             >  2:389221/1‑136 (MQ=255)

                                                                                                             gtgAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCt             >  1:269624/1‑136 (MQ=255)

                                                                                                             gtgAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCt             <  1:263621/136‑1 (MQ=255)

                                                                                                                       cACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCTGGCCCGTCAg   <  2:240761/136‑1 (MQ=255)

                                                                                                                        aCCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCTGGCCCGTCAGc  >  1:59055/1‑136 (MQ=255)

GTAGGCAATTAATACCGCATAGATCACTGCAGCGCCAATGACTACCCGCGGGTCAACCAGGAATCCGGCGACCACGGCGATGCACGCCAGCACCAGCGCCAGGCCGGATGTGAAGATGCCACCCGGGGTGCGGTAGGGGCGCTCCATCTTTGGCCGGCGGATGCGCAGGGTGATGTGCGCAGCCATCATCAGCACGTAAGACAGGGTAGCGCCAAACACTGCCACCAGGATCAGCAAGTCACCCTGGCCCGTCAGC  >  NC_002947/630337‑630592

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