Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,289,884	T→C	9.1%	W306R (TGG→CGG)	PP_2892 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,289,884	0	T	C	9.1%	119.1 / 2.5	45	W306R (TGG→CGG)	PP_2892	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base C (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GGCGAAGAAGTCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCACCGACACGCTCAGCAAG  >  NC_002947/3289737‑3290024                                                                                                                                                    |                                                                                                                                             ggCGAAGAAGTCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTg                                                                                                                                             <  6:34464/149‑1 (MQ=255) ggCGAAGAAGTCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTg                                                                                                                                             <  4:343067/149‑1 (MQ=255)           tCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGTAAAGCGaa                                                                                                                                   <  4:299255/149‑1 (MQ=255)           tCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGaa                                                                                                                                   >  6:102600/1‑149 (MQ=255)             ccGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACt                                                                                                                                 <  8:3365/149‑1 (MQ=255)             ccGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACt                                                                                                                                 >  5:13278/1‑149 (MQ=255)             ccGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCCCGCGGCAAGGCCAACTGGACCCGCCACAGGGTCAACCCACAAGGCCACAGCTATGGCCGGGAAACCGAACt                                                                                                                                 >  4:319860/1‑149 (MQ=255)               ggTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGt                                                                                                                               <  8:40704/149‑1 (MQ=255)                   gATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCgg                                                                                                                            <  5:190665/148‑1 (MQ=255)                    aTGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGGTGGGAACCCGAAGTGTCGGcc                                                                                                                          >  6:379641/1‑149 (MQ=255)                        tCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCcgc                                                                                                                      >  4:247570/1‑149 (MQ=255)                          aTCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCgcgc                                                                                                                    >  6:275711/1‑149 (MQ=255)                           tCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGAGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCAGGGAAACCGAACTGTCGGCCCCgcgc                                                                                                                    >  1:301782/1‑148 (MQ=255)                             tGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCgcgc                                                                                                                    >  7:198582/1‑146 (MQ=255)                              ggAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAg                                                                                                                 >  2:52016/1‑148 (MQ=255)                               gAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGa                                                                                                                <  2:367441/148‑1 (MQ=255)                               caaaCCGGTAAACCAGGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACTTCGGGCCGCGCGGCATGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCa                                                                                                                  <  7:234757/145‑1 (MQ=255)                                 aaCCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAAc                                                                                                              <  5:31260/148‑1 (MQ=255)                                       tAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGcc                                                                                                       <  5:275711/149‑1 (MQ=255)                                             aCGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGt                                                                                                 <  7:317054/149‑1 (MQ=255)                                                 aCCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCtgt                                                                                                    <  6:33065/142‑1 (MQ=255)                                                 aCCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCtgt                                                                                                    >  5:33065/1‑142 (MQ=255)                                                     tGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACAtt                                                                                         >  4:24369/1‑149 (MQ=255)                                                          cTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCt                                                                                    <  6:13278/149‑1 (MQ=255)                                                             ggCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGTGTCAACCCACAAGGCCACAGCTATGGCTGGGAATCCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTc                                                                                   <  2:26526/147‑1 (MQ=255)                                                                   ggCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACccg                                                                             <  5:379641/147‑1 (MQ=255)                                                                             ccgacGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGACCTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGAc                                                                 <  5:102600/149‑1 (MQ=255)                                                                              cgacGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGAc                                                                 >  6:352459/1‑148 (MQ=255)                                                                                   tCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTg                                                           <  8:198582/149‑1 (MQ=255)                                                                                    cGGGCCGCGCGGCAAGGCCAACTGGACCCGGCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTgc                                                          >  1:31556/1‑149 (MQ=255)                                                                                                          tGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGaa                                    >  2:341695/1‑149 (MQ=255)                                                                                                             aCCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAAc                                   >  7:34922/1‑147 (MQ=255)                                                                                                               ccGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACgagg                               <  3:247570/149‑1 (MQ=255)                                                                                                                 gCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACgaggag                             <  5:44347/149‑1 (MQ=255)                                                                                                                 gCCACTGGGTCAACCCACAAGGCCACAGCTATGGCCGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACgaggag                             >  1:218810/1‑149 (MQ=255)                                                                                                                     cTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTg                         <  3:364003/149‑1 (MQ=255)                                                                                                                        ggTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCg                      <  5:158503/149‑1 (MQ=255)                                                                                                                        ggTCAACCCACAAGGCCACAGCTATGGCCGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGGCg                      <  3:319860/149‑1 (MQ=255)                                                                                                                         gTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCgg                     >  8:277318/1‑149 (MQ=255)                                                                                                                            aaCCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCAc                  >  2:339990/1‑149 (MQ=255)                                                                                                                               ccACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCACCGa               >  3:163476/1‑149 (MQ=255)                                                                                                                                      gCCACAGCTATGGCCGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCACCGACACGCTc        <  2:218810/149‑1 (MQ=255)                                                                                                                                            gCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCACCGACACGCTCAGCAAg  <  7:277318/149‑1 (MQ=255)                                                                                                                                              tATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGc                            <  3:102755/121‑1 (MQ=255)                                                                                                                                              tATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGc                            >  4:102755/1‑121 (MQ=255)                                                                                                                                                    |                                                                                                                                             GGCGAAGAAGTCCCGGTGGATGTTCATCTGGAAACCCGTAAACCACGCACCATGGGCCTTGGCCTGGGCTTCTCGACCGACGTCGGGCCGCGCGGCAAGGCCAACTGGACCCGCCACTGGGTCAACCCACAAGGCCACAGCTATGGCTGGGAAACCGAACTGTCGGCCCCGCGCCAGAACGTCGGCCTGTGGTATGACATTCCCCTCGACCCGCCGCTGACCGACAAGTTGCGCTTTGCCGGCGGCTACCAGAACGAGGAGCTTGCCGGCACCGACACGCTCAGCAAG  >  NC_002947/3289737‑3290024

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