Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,286	A→C	8.7%	intergenic (+86/+56)	aroP‑I → / ← PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,070,286	0	A	C	8.7%	140.3 / 5.3	46	intergenic (+86/+56)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base C (2/2);  total (22/24)
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.84e-01

GCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTC  >  NC_002947/1070147‑1070428                                                                                                                                                  |                                                                                                                                               gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                                                                          <  6:122221/148‑1 (MQ=255)             gtggCTTGCAAGCTGACGAAGGCTCGCCAGGCCAAGGGCGGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                                                             <  2:120993/149‑1 (MQ=38)                gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGt                                                                                                                           <  7:84656/148‑1 (MQ=255)                         cTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                                                                                                 >  5:202334/1‑149 (MQ=255)                          tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                                                                >  6:228928/1‑149 (MQ=255)                           gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                                                                          <  2:150705/122‑1 (MQ=255)                           gaagaaGGCTCGCCAGGCCAAGGGCTGAACGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                                                                          >  1:150705/1‑122 (MQ=37)                               aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                                                                           <  8:251483/149‑1 (MQ=255)                                    tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                                                                      >  8:136052/1‑149 (MQ=255)                                       ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                                                                                   >  1:114390/1‑149 (MQ=255)                                              aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                                                                            <  1:225838/149‑1 (MQ=255)                                               aGGGCGGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGGTTTCGGCGCCGGTGA‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                                                                                           <  2:223844/149‑1 (MQ=38)                                                 ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                                                                                         >  2:182002/1‑149 (MQ=255)                                                  gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                                                                        >  5:276544/1‑149 (MQ=255)                                                  gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTACCCTTCCATGCCCCt                                                                                        >  7:70431/1‑149 (MQ=255)                                                      aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                                                                                      <  1:182002/147‑1 (MQ=255)                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                 <  5:71546/149‑1 (MQ=255)                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                 <  6:276544/149‑1 (MQ=255)                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                 >  1:165450/1‑149 (MQ=255)                                                                  cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                                                   <  4:29165/106‑1 (MQ=255)                                                                  cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                                                   >  3:29165/1‑106 (MQ=255)                                                                      aGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                                                                     >  3:235375/1‑148 (MQ=255)                                                                       ggggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgca                                                                   <  2:114390/149‑1 (MQ=255)                                                                        gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                                                                  >  6:130887/1‑149 (MQ=255)                                                                          gcaTGCCTTTAGGCTTTCATCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg                                                                 >  4:246069/3‑148 (MQ=17)                                                                             tGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgcac                                                              <  4:235375/148‑1 (MQ=255)                                                                                ctttAGGCTTTCATCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                                          <  3:246069/147‑1 (MQ=18)                                                                                cGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                                          <  5:228928/149‑1 (MQ=255)                                                                                 gTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca                                                             >  1:30423/1‑145 (MQ=255)                                                                                 gTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca                                                             <  2:30423/145‑1 (MQ=255)                                                                                    aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgttgtg                                                       <  6:202334/148‑1 (MQ=255)                                                                                            ggCGCC‑‑TGT‑‑GCTCTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa                                               <  4:10467/148‑1 (MQ=255)                                                                                             gcgcC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAaca                                             >  1:219831/1‑149 (MQ=255)                                                                                                      tgagA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAg                                       <  7:136052/145‑1 (MQ=255)                                                                                                                aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                 >  5:129558/1‑98 (MQ=25)                                                                                                                aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                 <  6:129558/98‑1 (MQ=25)                                                                                                                     cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAg                          >  1:108886/1‑148 (MQ=255)                                                                                                                       cgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTAc                       <  1:186460/149‑1 (MQ=255)                                                                                                                               gcgcTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGTACa               <  3:123952/149‑1 (MQ=255)                                                                                                                                       tctcACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGaa                 >  4:193626/1‑139 (MQ=255)                                                                                                                                       tctcACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGaa                 <  3:193626/139‑1 (MQ=255)                                                                                                                                       tctcACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCaa        <  2:165450/148‑1 (MQ=255)                                                                                                                                           acaGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCt   >  4:56613/1‑149 (MQ=255)                                                                                                                                           acaGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCt   >  8:55230/1‑149 (MQ=255)                                                                                                                                             aGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACCGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTc  >  3:53209/1‑148 (MQ=255)                                                                                                                                             aGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAAAGCAACACCGTAGCGAGCACCGGCAGTACATAGAAAAGCATCAATCGCTc  >  7:167600/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                               GCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTC  >  NC_002947/1070147‑1070428

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