Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,439	G→C	8.8%	intergenic (‑350/+338)	PP_5424 ← / ← PP_0050	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	58,439	0	G	C	8.8%	175.1 / 7.3	57	intergenic (‑350/+338)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base C (3/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

GTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTC  >  NC_002947/58295‑58576                                                                                                                                                 |                                                                                                                                           gTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTcc                                                                                                                                        >  6:148740/1‑149 (MQ=255)  tGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCct                                                                                                                                       <  5:148740/149‑1 (MQ=255)    aGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCctct                                                                                                                                     <  3:122299/149‑1 (MQ=255)      gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGa                                                                                                                                   <  8:28624/149‑1 (MQ=255)       ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                  <  3:102743/149‑1 (MQ=255)        cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                 >  1:96501/1‑149 (MQ=255)          gccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAagta                                                                                                                               <  6:324679/149‑1 (MQ=255)            cATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAagtagt                                                                                                                             <  1:335034/149‑1 (MQ=255)             aTCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTg                                                                                                                            <  2:96501/149‑1 (MQ=255)                    tGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcggc‑gg                                                                                                                    >  1:19418/1‑149 (MQ=255)                     gCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGAAATTGAGGGCTGAAAATGCGCATTTGTCGGTCCTGGGTACAAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑ggg                                                                                                                   >  1:60965/1‑149 (MQ=255)                      cATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑ggg                                                                                                                   <  2:19418/148‑1 (MQ=255)                        tCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTc                                                                                                                                         <  8:295011/125‑1 (MQ=255)                        tCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTc                                                                                                                                         >  7:295011/1‑125 (MQ=255)                         cAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCAtt                                                                                                               >  7:15189/1‑149 (MQ=255)                         cAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAA‑TTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTg                                                                                                              >  3:361992/1‑149 (MQ=255)                                 gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                  >  6:254700/1‑121 (MQ=255)                                 gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                  <  5:254700/123‑3 (MQ=255)                                 gattgaCACGCGTGCCGCTTGGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTc                                                                                                        <  4:26234/148‑1 (MQ=255)                                      acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTc                                                                                                  <  2:125320/149‑1 (MQ=255)                                       cacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCt                                                                                                 <  5:7474/149‑1 (MQ=255)                                          gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGt                                                                                                         >  8:10235/1‑138 (MQ=255)                                          gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGt                                                                                                         <  7:10235/138‑1 (MQ=255)                                             tGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAAcc                                                                                            >  6:100120/1‑148 (MQ=255)                                                  cTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGt                                                                                      >  5:95557/1‑149 (MQ=255)                                                            ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                              >  8:62940/1‑147 (MQ=255)                                                             cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                           <  8:15189/149‑1 (MQ=255)                                                                 ctaTTGAGGGCTGCAACGGCTCTTGTGTCGGGCATGGGTAGGAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAg                                                                        <  6:74863/146‑1 (MQ=255)                                                                       aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAg                                                                  <  6:52015/148‑1 (MQ=255)                                                                                       aTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActc                                                 >  3:74611/1‑149 (MQ=255)                                                                                            tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                            >  5:65255/1‑149 (MQ=255)                                                                                              gggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATc                                                              <  4:275480/129‑1 (MQ=255)                                                                                              gggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATc                                                              >  3:275480/1‑129 (MQ=255)                                                                                                    gggTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                >  5:323/1‑137 (MQ=255)                                                                                                    gggTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                <  6:323/137‑1 (MQ=255)                                                                                                              aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                          <  2:60965/149‑1 (MQ=255)                                                                                                               tGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                          >  7:233174/1‑148 (MQ=255)                                                                                                               tGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                          <  8:233174/148‑1 (MQ=255)                                                                                                                 aGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAcgc                        <  7:62940/148‑1 (MQ=255)                                                                                                                     cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcg                    >  1:91907/1‑149 (MQ=255)                                                                                                                      ccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggt                  >  1:166125/1‑149 (MQ=255)                                                                                                                      ccGCAA‑TTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtc                 <  4:361992/149‑1 (MQ=255)                                                                                                                       cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggtc                 >  2:331792/1‑149 (MQ=255)                                                                                                                                  cagcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCTTCATCCCTTAc                                     >  1:368954/1‑118 (MQ=255)                                                                                                                                  cagcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCTTCATCCCTTAc                                     <  2:368954/118‑1 (MQ=255)                                                                                                                                   ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt     >  3:229172/3‑149 (MQ=255)                                                                                                                                   agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt     >  8:298430/1‑149 (MQ=255)                                                                                                                                   agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt     >  4:371102/1‑149 (MQ=255)                                                                                                                                    gcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaaa                                <  3:172193/121‑1 (MQ=255)                                                                                                                                    gcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaaa                                >  4:172193/1‑121 (MQ=255)                                                                                                                                     cagAGGAACCAGCTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGa                                                       >  6:200680/1‑97 (MQ=255)                                                                                                                                     cagAGGAACCAGCTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGa                                                       <  5:200680/97‑1 (MQ=255)                                                                                                                                       gagGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc  >  8:338679/1‑148 (MQ=255)                                                                                                                                        agGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc  >  1:335969/1‑147 (MQ=255)                                                                                                                                        agGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc  >  6:154297/1‑147 (MQ=255)                                                                                                                                          gAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGcc                                                                                  <  3:303573/65‑1 (MQ=255)                                                                                                                                          gAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGcc                                                                                  >  4:303573/1‑65 (MQ=255)                                                                                                                                                 |                                                                                                                                           GTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGC‑GGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTC  >  NC_002947/58295‑58576

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