BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,290,355	G→A	16.9%	A50V (GCA→GTA)	ccmB ←	heme exporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,290,355	0	G	A	16.9%	1222.6 / 232.2	493	A50V (GCA→GTA)	ccmB	heme exporter subunit
Reads supporting (aligned to +/- strand): ref base G (189/219); new base A (35/48); total (225/268)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 510 total aligned reads displayed.

GTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAGGAA  >  CP009273/2290288‑2290422

gTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTg                                                                     <  1:23491145/68‑1 (MQ=255)

  cACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACa                                                                   <  1:36851663/68‑1 (MQ=255)

  cACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACa                                                                   <  1:7168035/68‑1 (MQ=255)

       aaCAGGCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATAc                                                               <  1:8218644/67‑1 (MQ=255)

        aCAGTCGTTCCAGCGCCAGCAAGGATGATAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgc                                                             >  1:26798859/1‑68 (MQ=255)

        aCAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgc                                                             >  1:4532496/1‑68 (MQ=255)

        aCAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgc                                                             >  1:16477349/1‑68 (MQ=255)

        aCAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgc                                                             >  1:346315/1‑68 (MQ=255)

           gTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGcc                                                          >  1:8816014/1‑68 (MQ=255)

             cGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgcg                                                            <  1:3089647/64‑1 (MQ=255)

             cGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgcg                                                            <  1:21520561/64‑1 (MQ=255)

             cGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATAcgcg                                                            <  1:35000157/64‑1 (MQ=255)

             cGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCca                                                         >  1:7293215/1‑67 (MQ=255)

             cGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCca                                                         >  1:3529207/1‑67 (MQ=255)

                tCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCcagcag                                                     <  1:5756766/68‑1 (MQ=255)

                tCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCcagca                                                      <  1:22870435/67‑1 (MQ=255)

                tCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCcagca                                                      <  1:19720001/67‑1 (MQ=255)

                 ccAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGt                                                    >  1:5163390/1‑68 (MQ=255)

                 ccAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGt                                                    >  1:10257199/1‑68 (MQ=255)

                  cAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGtt                                                   <  1:12706525/68‑1 (MQ=255)

                      gcCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCg                                                >  1:32634984/1‑67 (MQ=255)

                      gcCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCg                                                >  1:30412260/1‑67 (MQ=255)

                      gcCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCg                                                >  1:672380/1‑67 (MQ=255)

                        cAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGc                                              <  1:16747333/67‑1 (MQ=255)

                        cAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGc                                              <  1:17080999/67‑1 (MQ=255)

                           cAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCg                                          >  1:3145317/1‑68 (MQ=255)

                            aaGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTcc                                           <  1:22546941/66‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:14067283/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:26947181/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:27480843/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:17647193/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:6245056/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:33694931/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:29958030/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:17436156/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:30754204/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:3546886/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:21247206/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:1469711/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:32424671/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:18740551/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCGCCGGa                                        <  1:20351689/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:14908469/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:27960516/68‑1 (MQ=255)

                             aGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGa                                        <  1:10034276/68‑1 (MQ=255)

                                aTGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGAcc                                      <  1:34355752/67‑1 (MQ=255)

                                 tGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGa                                    <  1:11734304/68‑1 (MQ=255)

                                 tGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGa                                    <  1:33020016/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:19332649/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:28612746/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:11251580/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:15809669/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:20034229/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:10372855/68‑1 (MQ=255)

                                  gAAAGCAGCGCAGCAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   <  1:11017789/68‑1 (MQ=255)

                                    aaGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGAt                                   >  1:24913905/1‑66 (MQ=255)

                                      gcagcgcagcAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTg                               >  1:32159555/1‑68 (MQ=255)

                                      gcagcgcagcAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTg                               >  1:2912811/1‑68 (MQ=255)

                                        agcgcagcAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGa                              <  1:13792402/67‑1 (MQ=255)

                                        agcgcagcAACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGa                              <  1:23917461/67‑1 (MQ=255)

                                           gcagcaACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGc                            <  1:31745533/66‑1 (MQ=255)

                                           gcagcaACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGc                            <  1:2529326/66‑1 (MQ=255)

                                           gcagcaACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGc                            <  1:35981904/66‑1 (MQ=255)

                                             agcaACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGa                         <  1:29580483/67‑1 (MQ=255)

                                             agcaACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGa                         <  1:22253480/67‑1 (MQ=255)

                                             agcaACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGa                         <  1:28285346/67‑1 (MQ=255)

                                             agcaACCCAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGa                         <  1:6667234/67‑1 (MQ=255)

                                                  cccAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAgg                   <  1:36347415/68‑1 (MQ=255)

                                                  cccAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAgg                   <  1:27607320/68‑1 (MQ=255)

                                                  cccAGATAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAgg                   <  1:5940351/68‑1 (MQ=255)

                                                   ccAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAggg                  <  1:29369527/68‑1 (MQ=255)

                                                   ccAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAggg                  <  1:10629444/68‑1 (MQ=255)

                                                   ccAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAggg                  <  1:769957/68‑1 (MQ=255)

                                                       aTAATGCCCGGTGCTATACGCCCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAACGGTAAt              >  1:14503832/1‑68 (MQ=255)

                                                       aTAATGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAt              >  1:37016968/1‑68 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAtt             >  1:26311075/1‑68 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAt              <  1:31319152/67‑1 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAt              <  1:33559035/67‑1 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAt              <  1:34123193/67‑1 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAt              <  1:12434025/67‑1 (MQ=255)

                                                        tAATGCCCGGTGCAATACGCGCCAGCACTTGCGGCTCCGGACCGATACTCAGCCGAAAAAGGGTAAtt             >  1:4805345/1‑68 (MQ=255)

                                                         aaTGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAtt             >  1:18319496/1‑67 (MQ=255)

                                                         aaTGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAAtt             >  1:35744216/1‑67 (MQ=255)

                                                           tGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACa          <  1:851043/68‑1 (MQ=255)

                                                           tGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACa          <  1:9107300/68‑1 (MQ=255)

                                                           tGCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACa          <  1:15413683/68‑1 (MQ=255)

                                                            gCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:24392493/1‑68 (MQ=255)

                                                            gCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:7631790/1‑68 (MQ=255)

                                                            gCCCGGTACAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:20468392/1‑68 (MQ=255)

                                                             cccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTaaaa         >  1:23279225/1‑67 (MQ=255)

                                                             cccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:33276083/1‑67 (MQ=255)

                                                             cccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:25530662/1‑67 (MQ=255)

                                                             cccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACaa         >  1:16080841/1‑67 (MQ=255)

                                                              ccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAAt        <  1:27295540/67‑1 (MQ=255)

                                                              ccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATc       >  1:20987964/1‑68 (MQ=255)

                                                              ccGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATc       >  1:13125832/1‑68 (MQ=255)

                                                                  tGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAgg    <  1:20612476/67‑1 (MQ=255)

                                                                   gCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAGgaa  >  1:25985746/1‑68 (MQ=255)

                                                                   gCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAGgaa  >  1:25822768/1‑68 (MQ=255)

                                                                   gCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAGgaa  >  1:18057929/1‑68 (MQ=255)

GTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGCAGCAACCCAGATAATGCCCGGTGCAATACGCGCCAGCAGTTGCGGCTCCGGACCGATACTGAGCGGAAAAAGGGTAATTACAATCAGGAA  >  CP009273/2290288‑2290422

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