Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,361,719	C→G	14.6%	G720R (GGC→CGC)	PP_2076 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,361,719	0	C	G	14.6%	105.8 / 9.4	41	G720R (GGC→CGC)	PP_2076	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/15);  new base G (3/3);  total (23/18)
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

ACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCATCC  >  NC_002947/2361573‑2361862                                                                                                                                                   |                                                                                                                                                 aCTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGt                                                                                                                                                <  4:69780/149‑1 (MQ=255)           ctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCGGTACGGCGCGaa                                                                                                                                      <  1:260780/149‑2 (MQ=255)                         aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcaca                                                                                                                         >  4:182590/1‑148 (MQ=255)                                             accacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAg                                                                                                    >  3:3190/1‑149 (MQ=255)                                                          cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                       >  8:201815/1‑149 (MQ=255)                                                               ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGac                                                                                                                                      <  2:246709/97‑1 (MQ=255)                                                               ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGac                                                                                                                                      >  1:246709/1‑97 (MQ=255)                                                                             aaaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                          >  3:178807/1‑127 (MQ=38)                                                                             aaaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                          <  4:178807/127‑1 (MQ=38)                                                                              aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCa                                                                     <  3:182590/147‑1 (MQ=37)                                                                               aCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  <  7:201815/149‑1 (MQ=38)                                                                                         cATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCtt                                                        >  5:19722/1‑149 (MQ=39)                                                                                                   aaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  <  7:185773/129‑1 (MQ=40)                                                                                                   aaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  >  8:185773/1‑129 (MQ=40)                                                                                                    aCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCAt                                             <  6:19722/149‑1 (MQ=39)                                                                                                     cGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATc                                            >  5:9833/1‑149 (MQ=39)                                                                                                       cccTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCAc                                          >  3:262694/1‑149 (MQ=39)                                                                                                            gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                     <  4:3190/149‑1 (MQ=38)                                                                                                                     ctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCAGCGCATCACTGTCAGG‑CAGTTcg                           >  6:188104/1‑149 (MQ=255)                                                                                                                      tgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTcgc                          >  8:36496/1‑149 (MQ=39)                                                                                                                       gctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGTTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGtt                                                               <  6:17852/112‑1 (MQ=25)                                                                                                                        ctcctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTcgc                          <  6:9833/147‑1 (MQ=39)                                                                                                                            tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                        <  2:158059/82‑1 (MQ=38)                                                                                                                            tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                        >  1:158059/1‑82 (MQ=38)                                                                                                                            tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCtt                    >  6:219631/1‑149 (MQ=39)                                                                                                                                 cGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGaa                <  7:36496/148‑1 (MQ=39)                                                                                                                                  ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAg                                    <  1:247103/128‑1 (MQ=38)                                                                                                                                  ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAg                                    >  2:247103/1‑128 (MQ=38)                                                                                                                                     cacaCTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCa           >  3:66618/1‑149 (MQ=39)                                                                                                                                        ggtggtCGCGCGGTACGGCGCGAACACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCAt          >  8:66807/6‑147 (MQ=255)                                                                                                                                          tggtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAg       >  3:257129/4‑148 (MQ=255)                                                                                                                                          tGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAg       >  6:77197/1‑148 (MQ=39)                                                                                                                                            gtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCAt    >  2:160787/2‑149 (MQ=255)                                                                                                                                             tCGCGCGGTACGGCGCGAACACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCATc   <  7:66807/149‑1 (MQ=255)                                                                                                                                              cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCt                     >  6:88885/1‑130 (MQ=40)                                                                                                                                              cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCt                     <  5:88885/130‑1 (MQ=40)                                                                                                                                              cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCATcc  >  5:255030/1‑149 (MQ=39)                                                                                                                                              cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCATcc  >  2:24907/1‑149 (MQ=39)                                                                                                                                               gcgcGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCCAGTTCGCGCGCTTTGAAGGCCATCAGCATc   <  1:160787/148‑1 (MQ=255)                                                                                                                                               gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAccc                                                   >  7:78088/1‑100 (MQ=40)                                                                                                                                               gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAccc                                                   <  8:78088/100‑1 (MQ=40)                                                                                                                                                   |                                                                                                                                                 ACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGG‑CAGTTCGCGCGCTTTGAAGGCCATCAGCATCC  >  NC_002947/2361573‑2361862

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