Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,872	G→A	8.9%	C303Y (TGC→TAC)	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,872	0	G	A	8.9%	125.9 / 2.6	45	C303Y (TGC→TAC)	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (2/2);  total (21/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 = 8.97e-01

GGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACA  >  NC_002947/3593727‑3594013                                                                                                                                                    |                                                                                                                                              ggTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCg                                                                                                                                            >  4:34673/1‑149 (MQ=255)      ccTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATcc                                                                                                                                       <  6:232742/149‑1 (MQ=255)          gAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCAt                                                                                                                                     <  3:34673/147‑1 (MQ=255)           aaCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGAc                                                                                                                                  >  5:232664/1‑149 (MQ=255)               catGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACAGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGc                                                                                                                              >  4:222481/3‑149 (MQ=255)                ttGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTACCGCATCCATGACCCGCt                                                                                                                             >  2:86720/1‑149 (MQ=255)                ttGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTACCGCATCCATGACCCGCt                                                                                                                             >  2:91142/1‑149 (MQ=255)                    cGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATcc                                                                                                                                       >  8:330540/1‑135 (MQ=255)                    cGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATcc                                                                                                                                       <  7:330540/135‑1 (MQ=255)                      ctgatcGAGCACGGCGCGCTCGATCAGGTGAGTGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                        >  1:212353/6‑149 (MQ=255)                               cgcgcCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTc                                                                                                              <  1:149698/149‑1 (MQ=255)                                   ccGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                          <  5:14364/149‑1 (MQ=255)                                    cGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTcaca                                                                                                           <  3:239503/147‑1 (MQ=255)                                    cGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTcaca                                                                                                           >  4:239503/1‑147 (MQ=255)                                     gTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGc                                                                                                        <  7:130703/149‑1 (MQ=255)                                        cTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTa                                                                                                     >  3:242182/1‑149 (MQ=255)                                         tCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTAt                                                                                                    <  6:232664/149‑1 (MQ=255)                                           gATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGc                                                                                                  >  7:246279/1‑149 (MQ=255)                                                  tGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCa                                                                                                                                      <  6:244351/106‑1 (MQ=255)                                                        gCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgct                                                                                      >  7:231621/1‑149 (MQ=255)                                                          cTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctgg                                                                                    <  4:242182/149‑1 (MQ=255)                                                           tGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctgg                                                                                    >  1:217496/1‑149 (MQ=255)                                                                cGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAAc                                                                              <  8:231621/149‑1 (MQ=255)                                                                    gCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTggg                                                                          >  3:281285/1‑149 (MQ=255)                                                                    gCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTgg                                                                           <  8:246279/148‑1 (MQ=255)                                                                     ccGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGc                                                                         >  3:90986/1‑149 (MQ=255)                                                                          gggCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTACCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGaca                                                                    <  1:91142/149‑1 (MQ=255)                                                                           ggCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTACCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGaca                                                                    <  1:86720/148‑1 (MQ=255)                                                                              ccGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGAt                                                                <  2:217496/149‑1 (MQ=255)                                                                                    gggTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTa                                                          <  3:222481/149‑1 (MQ=255)                                                                                    gggTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTa                                                          <  2:212353/149‑1 (MQ=255)                                                                                            ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCg                                                   >  6:361951/1‑149 (MQ=255)                                                                                                 ggCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTggg                                                                          >  2:281212/1‑121 (MQ=255)                                                                                                 ggCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTggg                                                                          <  1:281212/121‑1 (MQ=255)                                                                                                 ggCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                             <  5:24635/134‑1 (MQ=255)                                                                                                 ggCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                             >  6:24635/1‑134 (MQ=255)                                                                                                               aTCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCct                                <  8:99957/149‑1 (MQ=255)                                                                                                                tCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTc                                  >  5:187406/1‑146 (MQ=255)                                                                                                                tCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTc                                  <  6:187406/146‑1 (MQ=255)                                                                                                                                      gCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTca                            >  7:113287/1‑130 (MQ=255)                                                                                                                                      gCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTca                            <  8:113287/130‑1 (MQ=255)                                                                                                                                          atgatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTcc       <  4:90986/147‑1 (MQ=255)                                                                                                                                          atgatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCg      >  7:65167/1‑148 (MQ=255)                                                                                                                                          atgatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGt     <  5:361951/149‑1 (MQ=255)                                                                                                                                             atCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTaca  >  6:234279/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                              GGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGA‑TGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGA‑GGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACA  >  NC_002947/3593727‑3594013

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