Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,263	G→C	25.7%	intergenic (+63/+79)	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,263	0	G	C	25.7%	67.6 / 23.3	39	intergenic (+63/+79)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base C (4/6);  total (17/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 = 4.10e-01

CTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACAT  >  NC_002947/1070119‑1070409                                                                                                                                                      |                                                                                                                                                   ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTGTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                                                   >  1:85633/1‑149 (MQ=255) ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTA‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcgggc                                                                                                                                                  >  5:269779/1‑149 (MQ=255)    ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                                                   <  7:52790/145‑1 (MQ=255)       aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGACGCGCGGGCGGCGCt                                                                                                                                            >  4:341747/1‑149 (MQ=255)        aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                            <  3:289936/148‑1 (MQ=37)         tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                          <  6:35883/149‑1 (MQ=38)           cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGTGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                        <  6:225152/149‑1 (MQ=38)           cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                         <  2:240014/148‑1 (MQ=25)           cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                        >  5:225152/1‑149 (MQ=38)              aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                      <  5:118282/148‑1 (MQ=255)                   ggCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                <  3:341747/149‑1 (MQ=255)                       tGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                                                            <  8:359441/149‑1 (MQ=255)                        ggATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                                                            <  4:342999/148‑1 (MQ=255)                             ggtgtgctgtgGGTGGCTT‑TCATGCTTAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                      <  2:382338/147‑1 (MQ=38)                                 tgctgtgGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTGGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                                                    <  1:120837/147‑1 (MQ=255)                                     gtgggtggCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                                                   >  8:145075/1‑112 (MQ=255)                                        ggtggCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                                                           <  1:241300/149‑1 (MQ=255)                                            gCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGACGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                                                       <  4:246354/149‑1 (MQ=38)                                                 tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                                                   >  4:212940/1‑149 (MQ=255)                                                 tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                                                   >  4:140349/1‑149 (MQ=255)                                                       tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                                             >  5:349853/1‑149 (MQ=255)                                                       tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                            <  2:262359/102‑1 (MQ=18)                                                       tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                            >  1:262359/1‑102 (MQ=18)                                                          agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                          <  3:140349/149‑1 (MQ=255)                                                               gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                                                     >  5:171662/1‑149 (MQ=255)                                                                 tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                           <  1:173075/93‑1 (MQ=25)                                                                 tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                           >  2:173075/1‑93 (MQ=25)                                                                     cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                                               >  6:85975/1‑149 (MQ=255)                                                                         ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                      <  4:343187/106‑1 (MQ=255)                                                                         ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                      >  3:343187/1‑106 (MQ=255)                                                                          cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                                                          <  2:85633/149‑1 (MQ=255)                                                                                                  aaGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                                                  >  1:56303/1‑149 (MQ=255)                                                                                                   aGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                                                  <  6:269779/148‑1 (MQ=255)                                                                                                              gTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg                                      >  6:325707/1‑149 (MQ=255)                                                                                                                  ggTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                  <  2:56303/149‑1 (MQ=255)                                                                                                                         ggCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc                                 <  8:43961/143‑1 (MQ=255)                                                                                                                         ggCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc                                 >  7:43961/1‑143 (MQ=255)                                                                                                                                   tgagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                  <  3:157690/131‑1 (MQ=255)                                                                                                                                   tgagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                  >  4:157690/4‑134 (MQ=255)                                                                                                                                   tgagATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                   <  3:212940/146‑1 (MQ=255)                                                                                                                                       aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAg                >  7:347721/1‑148 (MQ=255)                                                                                                                                         cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAc             >  5:361665/1‑149 (MQ=38)                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCg           >  7:181424/1‑149 (MQ=255)                                                                                                                                                    cgggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACAt  <  5:85975/149‑1 (MQ=255)                                                                                                                                                     gggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa   >  1:45381/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACAT  >  NC_002947/1070119‑1070409

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