Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	26,837:1	+A	9.5%	intergenic (+33/+152)	PP_0021 → / ← PP_0022	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	26,837	1	.	A	9.5%	153.4 / 10.1	42	intergenic (+33/+152)	PP_0021/PP_0022	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (18/20);  new base A (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGT  >  NC_002947/26694‑26985                                                                                                                                                 |                                                                                                                                                     tGGGTCGATTTCCTTATCCGGATGTCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACcac                                                                                                                                                <  1:116806/150‑1 (MQ=255)     tCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTg                                                                                                                                            <  2:360472/150‑1 (MQ=255)         aTCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTGTCACTGT‑AACCACAGTGGttt                                                                                                                                        >  2:409966/1‑150 (MQ=255)                aTCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGttaca                                                                                                                                      >  2:125555/1‑143 (MQ=255)                aTCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGttaca                                                                                                                                      <  1:125555/146‑4 (MQ=255)                      aTGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                           <  2:224358/150‑1 (MQ=255)                          cGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTgg                                                                                                                       <  1:409966/150‑1 (MQ=255)                                    tcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGt                                                                                                             >  1:103483/1‑150 (MQ=255)                                       ttATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGc                                                                                                          >  2:199627/1‑150 (MQ=255)                                             cAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGac                                                                                                    >  1:544552/1‑150 (MQ=255)                                                 aTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACAcgc                                                                                                <  1:274750/150‑1 (MQ=255)                                                      tAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCt                                                                                           >  2:148758/1‑150 (MQ=255)                                                         agTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                        <  1:304005/150‑1 (MQ=255)                                                           tGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCa                                                                                      <  1:148758/150‑1 (MQ=255)                                                                     ccGGGAGGAACTGGTCGTTCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                            <  1:610974/150‑1 (MQ=255)                                                                            gaACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAgcgcg                                                                     <  1:199627/150‑1 (MQ=255)                                                                               cTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGaaaaaa                                                                                                                            <  2:48370/92‑1 (MQ=255)                                                                               cTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGaaaaaa                                                                                                                            >  1:48370/1‑92 (MQ=255)                                                                                             ggAGCTGCATCCGGTCTGGAGATGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                    <  1:124625/150‑1 (MQ=37)                                                                                                 cTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGt                                                >  1:44636/1‑150 (MQ=255)                                                                                                      tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                           >  1:787127/1‑150 (MQ=255)                                                                                                      tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                           >  1:787134/1‑150 (MQ=255)                                                                                                          gCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                        >  1:794082/1‑149 (MQ=255)                                                                                                          gCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                        <  2:794082/149‑1 (MQ=255)                                                                                                            cTGAGTAGTCCGGTTTTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGaa                                     <  2:44636/150‑1 (MQ=255)                                                                                                            cTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGaa                                     <  2:103483/150‑1 (MQ=255)                                                                                                                    gCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtg                             >  1:327269/1‑150 (MQ=255)                                                                                                                       tgttgtAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCGCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctg                          >  1:314305/1‑150 (MQ=255)                                                                                                                        gttgtAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctgc                         >  1:18512/1‑150 (MQ=255)                                                                                                                           gtAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                           >  1:748460/1‑49 (MQ=255)                                                                                                                           gtAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                           <  2:748460/49‑1 (MQ=255)                                                                                                                                    tttttttCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGTTCGGGGCCTTga             >  1:50375/1‑149 (MQ=255)                                                                                                                                       ttttCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAc          >  2:38730/1‑150 (MQ=255)                                                                                                                                         ttCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                  >  1:586187/1‑124 (MQ=255)                                                                                                                                         ttCACTGT‑AACCACAGTGGTTTAAAGTGAATAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                  <  2:586187/124‑1 (MQ=255)                                                                                                                                          tCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                        >  2:726563/1‑69 (MQ=255)                                                                                                                                          tCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                        <  1:726563/69‑1 (MQ=255)                                                                                                                                            aCTGTAAACCACTGTGG‑TTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                            <  2:335522/95‑1 (MQ=25)                                                                                                                                            aCTGTAAACCACTGTGG‑TTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                            >  1:335522/1‑95 (MQ=25)                                                                                                                                               gT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgt  <  2:787127/150‑1 (MQ=255)                                                                                                                                               gT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgt  <  2:787134/150‑1 (MQ=255)                                                                                                                                               gT‑AACCACAGTGGGTTAAAGTAAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgt  <  2:272959/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGT‑AACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGT  >  NC_002947/26694‑26985

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