BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	594,985	G→A	11.7%	A391V (GCT→GTT)	PP_0506 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	594,985	0	G	A	11.7%	87.7 / 3.5	34	A391V (GCT→GTT)	PP_0506	ABC transporter permease
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 = 1.00e+00

CGTGCCCGCGAAGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGCCACGATCCCGACCGAGGCCAGTGAC  >  NC_002947/594846‑595125

cGTGCCCGCGAAGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGt                                                                                                                                    <  2:446080/150‑1 (MQ=255)

  tGCCCGCGAAGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTcc                                                                                                                                  <  1:447013/150‑1 (MQ=255)

          aaGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAg                                                                                                                                      >  2:52350/1‑138 (MQ=255)

          aaGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAg                                                                                                                                      <  1:52350/138‑1 (MQ=255)

          aaGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATgg                                                                                                                          <  2:212234/150‑1 (MQ=255)

          aaGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATgg                                                                                                                          <  1:471607/150‑1 (MQ=255)

                  ggTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCgg                                                                                                                  <  2:33013/150‑1 (MQ=255)

                        gtgtTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGc                                                                                                            >  2:533510/1‑150 (MQ=255)

                        gtgtTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGc                                                                                                            >  1:416939/1‑150 (MQ=255)

                        gggtTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGCTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGc                                                                                                            >  2:243859/3‑150 (MQ=255)

                                  aTTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAg                                                                                                  >  2:467206/1‑150 (MQ=255)

                                    tACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTa                                                                                                <  2:426108/150‑1 (MQ=255)

                                       gcgcACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACaa                                                                                             >  1:20123/1‑150 (MQ=255)

                                       gcgcACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACaa                                                                                             >  1:531441/1‑150 (MQ=255)

                                             tCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCa                                                                                                                                           <  2:140193/98‑1 (MQ=255)

                                             tCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCa                                                                                                                                           >  1:140193/1‑98 (MQ=255)

                                                        tGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCa                                                                                                                             <  1:175776/101‑1 (MQ=255)

                                                        tGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCa                                                                                                                             >  2:175776/1‑101 (MQ=255)

                                                           aCAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTg                                                                                   <  2:272042/140‑1 (MQ=255)

                                                           aCAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTg                                                                                   >  1:272042/1‑140 (MQ=255)

                                                                 cGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACACACCATAGTTGGCCTGCACATATCCCt                                                                   <  1:135104/150‑1 (MQ=255)

                                                                                   ggTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCg                                                 >  2:355913/1‑150 (MQ=255)

                                                                                    gTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAACCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGc                                                >  2:20211/1‑150 (MQ=255)

                                                                                      aaGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCat                                              <  1:4310/150‑1 (MQ=255)

                                                                                         cccGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAACCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCAtata                                           >  2:320445/1‑150 (MQ=255)

                                                                                             ccACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAACCCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGcaa                                       <  2:362573/150‑1 (MQ=255)

                                                                                             ccACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAACCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGcaa                                       <  1:20211/150‑1 (MQ=255)

                                                                                              cACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGcaag                                      >  1:498774/1‑150 (MQ=255)

                                                                                                         gCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGc                           <  1:168660/150‑1 (MQ=255)

                                                                                                         gCTGCCCATCAGCAACGCGGCCCCCAGGATGATAACCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGc                           <  1:320445/150‑1 (MQ=255)

                                                                                                            gCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGACAAGCCCGGCGAc                        >  1:212677/1‑150 (MQ=255)

                                                                                                              ccATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGCCACGa                      <  2:20123/150‑1 (MQ=255)

                                                                                                                    gCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGCCACGATCccga                <  1:467206/150‑1 (MQ=255)

                                                                                                                                  aGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGCCACGATCCCGACCGAGGCCAGTGAc  >  2:266667/1‑150 (MQ=255)

CGTGCCCGCGAAGCAGGCGGTGAGGTGTTTGATCATTACGCGCACTCAGAGGTGTATGGACAGGCCGTCAGCCAGCGATTGCCGGTAAGCCCGCCACGCCGGCACGCTGCCCATCAGCAACGCGGCCCCCAGGATGATAGCCAGCAAGGTCCATTCATGGGTACTCGGCAGAGCCAGTGGCAAGTACAAACCATAGTTGGCCTGCACATATCCCTGTGCCAGGGCGATACCCGCATATAGCAAGCCAAGCCCGGCCACGATCCCGACCGAGGCCAGTGAC  >  NC_002947/594846‑595125

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: `‑`