Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,298	A→C	9.1%	intergenic (+98/+44)	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,298	0	A	C	9.1%	134.3 / 5.6	44	intergenic (+98/+44)	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/20);  new base C (2/2);  total (22/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 = 6.92e-01

GTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCA  >  NC_002947/1070159‑1070438                                                                                                                                                  |                                                                                                                                             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)                   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                                                                                           <  6:276544/149‑1 (MQ=255)                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                           <  5:71546/149‑1 (MQ=255)                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                                           >  1:165450/1‑149 (MQ=255)                                                      cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                                                             >  3:29165/1‑106 (MQ=255)                                                      cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                                                             <  4:29165/106‑1 (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)                                                                                                                                          aaGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgcca  >  4:42885/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                             GTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCA  >  NC_002947/1070159‑1070438

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