Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,108,335	G→C	13.5%	W23C (TGG→TGC)	PP_5358 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,108,335	0	G	C	13.5%	117.9 / 7.7	44	W23C (TGG→TGC)	PP_5358	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base C (3/3);  total (22/22)
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

CGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGT  >  NC_002947/6108187‑6108481                                                                                                                                                     |                                                                                                                                                   cgcagcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTgg                                                                                                                                                    <  1:310355/149‑1 (MQ=255)  gcagcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTgg                                                                                                                                                    >  5:10168/1‑148 (MQ=255)        cagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAAc                                                                                                                                              >  8:187896/1‑148 (MQ=255)          gTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCAt                                                                                                                                           <  5:233097/149‑1 (MQ=255)           tATGAACTTGATTTCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTTGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCagg                                                                                                                                          <  6:54408/149‑3 (MQ=255)           tATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCAt                                                                                                                                           <  1:125259/148‑1 (MQ=255)              gAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGt                                                                                                                                       >  6:274761/1‑149 (MQ=255)                 cTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGc                                                                                                                                    >  1:131316/1‑149 (MQ=255)                   tGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                           >  6:256032/1‑139 (MQ=255)                   tGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                           <  5:256032/140‑2 (MQ=255)                    gATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCgg                                                                                                                                  <  6:10168/148‑2 (MQ=255)                          tACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAgg                                                                                                                           >  6:2868/1‑149 (MQ=255)                               gTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAAc                                                                                                                                              <  3:125114/125‑1 (MQ=255)                               gTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAAc                                                                                                                                              >  4:125114/1‑125 (MQ=255)                                tACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCt                                                                                                                       <  2:131316/147‑1 (MQ=255)                                tACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                      <  5:2868/148‑1 (MQ=255)                                          cAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCa                                                                                                            >  5:3835/1‑148 (MQ=255)                                              aaCGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAg                                                                                                       <  5:23094/149‑1 (MQ=255)                                                                cTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAgg                                                                                                                           <  6:133396/111‑1 (MQ=255)                                                                cTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGAc                                                                                     >  3:168693/1‑149 (MQ=255)                                                                 tACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGAcc                                                                                    >  3:160801/1‑149 (MQ=255)                                                                  aCCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                 <  6:32976/135‑1 (MQ=255)                                                                  aCCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                 >  5:32976/1‑135 (MQ=255)                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGAc                                                                                     <  5:274761/146‑1 (MQ=255)                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTg                                                                                  >  1:126176/1‑149 (MQ=255)                                                                       tttGGATCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCaaa                                                                                                          <  1:188500/121‑1 (MQ=255)                                                                       tttGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCaaa                                                                                                          >  2:188500/1‑121 (MQ=255)                                                                           gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGc                                                                           <  7:187896/148‑1 (MQ=255)                                                                           gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCt                                                                          >  8:175327/1‑149 (MQ=255)                                                                             gCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGc                                                                                                                          >  4:102033/1‑99 (MQ=255)                                                                              cATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGa                                                                       <  1:154302/149‑1 (MQ=255)                                                                                                cTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                     <  4:160801/149‑1 (MQ=255)                                                                                                     gaACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCACTCAAGCGTa                                                <  2:126176/149‑1 (MQ=255)                                                                                                                 ccTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTcc                                    >  7:17098/1‑149 (MQ=255)                                                                                                                               ccACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGc                      >  4:55621/1‑149 (MQ=255)                                                                                                                                   ggCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                   >  6:206251/1‑148 (MQ=255)                                                                                                                                          accTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt           >  6:161660/2‑149 (MQ=255)                                                                                                                                            atGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg         <  1:179864/148‑1 (MQ=255)                                                                                                                                            atGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg         >  2:179864/2‑149 (MQ=255)                                                                                                                                            atGGTTGTGCCACAACCATGGGTTGCGTGAAAAGCCCCTGGAACGGGAAAAACAGCAATGCGCCAAACTCGACCTGGACTTGCTTGACGATGCCGTAGCGCTCAAGAGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg         >  3:135279/2‑149 (MQ=255)                                                                                                                                              ggTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCgg                                  >  6:81600/1‑122 (MQ=255)                                                                                                                                              ggTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCgg                                  <  5:81600/122‑1 (MQ=255)                                                                                                                                               gTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                          <  3:252521/113‑1 (MQ=255)                                                                                                                                               gTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                          >  4:252521/1‑113 (MQ=255)                                                                                                                                                  gtgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc    <  7:175327/148‑1 (MQ=255)                                                                                                                                                   tgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGt  <  5:206251/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGT  >  NC_002947/6108187‑6108481

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