Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,839,682	G→C	9.2%	S71R (AGC→AGG)	PP_5529 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,839,682	0	G	C	9.2%	119.8 / 2.1	41	S71R (AGC→AGG)	PP_5529	epimerase
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base C (2/2);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTC  >  NC_002947/2839535‑2839825                                                                                                                                                     |                                                                                                                                                cAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGATGAGACGATGCTTTAGCGCCGCCTGGTAAGCCggg                                                                                                                                                 <  1:51858/149‑1 (MQ=255) cAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCggg                                                                                                                                                 <  6:277886/148‑1 (MQ=255)         tGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTAc                                                                                                                                        <  7:109334/149‑1 (MQ=255)           aCATGCACTAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTTGTAAGCCGGGCTGGCATAcc                                                                                                                                       <  5:38388/148‑1 (MQ=255)               gCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCt                                                                                                                                  >  4:194128/1‑149 (MQ=255)                             ataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCa                                                                                                                     >  6:365147/1‑148 (MQ=255)                             ataCGCTCATAACGACTAGTGCACGCCTGTCAGTTAGGTGATCTGTCGGTGCCTTCGAATAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGATTTAGCGCCGCCTGGTAATCCCGGCAGGCATACCAGGCTCGAGCGGCAGTCa                                                                                                                     >  2:314279/1‑148 (MQ=255)                                     aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                              <  5:224707/147‑1 (MQ=255)                                     aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                              <  7:80007/147‑1 (MQ=255)                                           aCTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                                                                                      <  8:210609/149‑1 (MQ=255)                                                   aCGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                              <  1:158059/149‑1 (MQ=255)                                                          tCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCt                                                                                       <  3:191597/149‑1 (MQ=255)                                                              ttAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                                                                                             >  6:25030/1‑91 (MQ=255)                                                              ttAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                                                                                             <  5:25030/91‑1 (MQ=255)                                                              ttAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                   >  5:180921/1‑149 (MQ=255)                                                                  gTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                               <  3:174370/149‑1 (MQ=255)                                                                   tGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATg                                                                              <  3:194128/149‑1 (MQ=255)                                                                    gATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGa                                                                                                    <  6:358587/126‑1 (MQ=255)                                                                    gATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGa                                                                                                    >  5:358587/1‑126 (MQ=255)                                                                         gtCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCtt                                                                        >  4:54625/1‑149 (MQ=255)                                                                          tCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGTCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCt                                                                                       <  1:331483/133‑1 (MQ=255)                                                                          tCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGGCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCt                                                                                       >  2:331483/1‑133 (MQ=255)                                                                            ggTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTcc                                                                      >  6:371543/1‑148 (MQ=255)                                                                                ccTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATg                                                                              <  3:15735/136‑1 (MQ=255)                                                                                ccTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATg                                                                              >  4:15735/1‑136 (MQ=255)                                                                                  ttCGATTACAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCt                                                               >  8:104563/1‑149 (MQ=255)                                                                                       ttAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGt                                                          <  5:365147/149‑1 (MQ=255)                                                                                                          cGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGAc                                       >  1:183302/1‑149 (MQ=255)                                                                                                           gCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                              <  7:309011/77‑1 (MQ=255)                                                                                                           gCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                              >  8:309011/1‑77 (MQ=255)                                                                                                              ccGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                   <  4:363455/149‑1 (MQ=255)                                                                                                              ccGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                   <  3:54625/149‑1 (MQ=255)                                                                                                                        cGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTg                          <  2:183302/149‑1 (MQ=255)                                                                                                                                  gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                        >  4:69190/1‑125 (MQ=255)                                                                                                                                  gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                        <  3:69190/125‑1 (MQ=255)                                                                                                                                    cgagccTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGc                                            <  2:328464/116‑1 (MQ=38)                                                                                                                                    cgagccTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGc                                            >  1:328464/4‑119 (MQ=38)                                                                                                                                            gTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGGGACGTCGTGCCCTTCTCGTGCGGCAGgc       <  7:104563/148‑1 (MQ=255)                                                                                                                                                gCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtc                             >  5:311394/1‑122 (MQ=255)                                                                                                                                                gCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtc                             <  6:311394/122‑1 (MQ=255)                                                                                                                                                gCCGGCCTGGCATACCAGGCTCGACCGGCTCTCAAATCTTCGAAACGTAGAATGACCGCTCCCTAGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACCTGTGACGACGTGCCCTTCTCGTGCGGCAGGCGcctc  >  2:297186/1‑149 (MQ=255)                                                                                                                                                 ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGa                                                       >  2:255267/1‑95 (MQ=255)                                                                                                                                                 ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGa                                                       <  1:255267/95‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGA‑GAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTC  >  NC_002947/2839535‑2839825

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