Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,973,403	Δ1 bp	12.1%	coding (471/1713 nt)	PP_2600 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,973,403	0	G	.	12.1%	97.9 / 7.8	33	coding (471/1713 nt)	PP_2600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base . (2/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTA  >  NC_002947/2973262‑2973549                                                                                                                                              |                                                                                                                                                   aTGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTgc                                                                                                                                             <  3:321806/149‑1 (MQ=255)           aGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTgg                                                                                                                                   <  6:63962/149‑1 (MQ=255)                      gTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGgcgc                                                                                                                        >  7:14530/1‑149 (MQ=255)                              cAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTcg                                                                                                                >  5:110321/1‑149 (MQ=255)                              cAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCAGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTcg                                                                                                                >  3:6369/1‑149 (MQ=255)                                        gCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCg                                                                                                       <  8:14530/148‑1 (MQ=255)                                        gCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCAT‑GTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCg                                                                                                       >  6:263792/1‑147 (MQ=255)                                         cacaCGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGTGGGCATCGTa                                                                                                     <  6:110321/149‑1 (MQ=255)                                         cacaCGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTa                                                                                                     <  2:294530/149‑1 (MQ=255)                                          acacGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAg                                                                                                    >  8:113575/1‑149 (MQ=255)                                            acGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAggg                                                                                                  >  2:263254/1‑149 (MQ=255)                                                         cTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCAT‑GTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTGGCGGGCATGGTAGGGCGAGGCGATCACCt                                                                                    >  4:118954/1‑149 (MQ=255)                                                         cTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCAGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCAcc                                                                                     <  4:6369/149‑1 (MQ=255)                                                          tttGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATCAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCt                                                                                    >  6:53743/1‑149 (MQ=255)                                                            tGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCt                                                                                    <  7:113575/147‑1 (MQ=255)                                                                          aTGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAggg                                                                     <  1:41530/148‑1 (MQ=255)                                                                             ttGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGAGATCACCTGGTGGCCGCCAAGGGCTgc                                                                 >  7:245091/1‑149 (MQ=255)                                                                              tGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCt                                                                                    <  4:102250/129‑1 (MQ=255)                                                                              tGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCt                                                                                    >  3:102250/1‑129 (MQ=255)                                                                                 ggAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCt                                                             >  2:55368/1‑149 (MQ=255)                                                                                  gAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTg                                                            >  2:61676/1‑149 (MQ=255)                                                                                  gAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTg                                                            >  8:274289/1‑149 (MQ=255)                                                                                   aaGGCGGCAAGGTCGGCACCGGCCAGGATCGGCATGCCCATGAACGAGGTTTCGCGCAT‑GTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGc                                                          <  3:118954/149‑1 (MQ=255)                                                                                                 ggCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgc                                              <  8:245091/148‑1 (MQ=255)                                                                                                       ggCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgc                                              <  5:85149/142‑1 (MQ=255)                                                                                                       ggCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgc                                              >  6:85149/1‑142 (MQ=255)                                                                                                                         cATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCgg                     >  2:98622/1‑149 (MQ=255)                                                                                                                         cATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCgg                     >  7:189350/1‑149 (MQ=255)                                                                                                                              aCGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAAc                >  7:92458/1‑149 (MQ=255)                                                                                                                                gAGGTTTCGCGCAT‑GTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGc             <  5:263792/149‑1 (MQ=255)                                                                                                                                 aGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGc             <  2:33452/149‑1 (MQ=255)                                                                                                                                        gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCg      >  5:276552/1‑149 (MQ=255)                                                                                                                                            aTGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTa  >  8:66745/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                   ATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGGCGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTA  >  NC_002947/2973262‑2973549

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