Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,218,367	G→C	13.1%	intergenic (‑193/+47)	PP_1958 ← / ← PP_5495	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,218,367	0	G	C	13.1%	103.8 / 6.8	38	intergenic (‑193/+47)	PP_1958/PP_5495	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/16);  new base C (3/2);  total (20/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 = 9.92e-01

GCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTCTCGAATCAT  >  NC_002947/2218224‑2218515                                                                                                                                                 |                                                                                                                                                     gcgcTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTtctc                                                                                                                                                 <  1:175874/149‑1 (MQ=255) gcgcTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTtctc                                                                                                                                                 <  3:167953/149‑1 (MQ=255)   gcTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCa                                                                                                                                               <  7:71839/149‑1 (MQ=255)       tGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCa                                                                                                                                            >  3:247677/1‑148 (MQ=255)        ggTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCAt                                                                                                                                           <  4:247677/148‑1 (MQ=255)              ggggATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGAcc                                                                                                                                     >  3:46359/1‑148 (MQ=255)                            ggAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTacac                                                                                                                      <  4:46359/149‑1 (MQ=255)                                gaagTTCGCCTAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCATGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAgg                                                                                                                  <  5:229829/149‑1 (MQ=255)                                 aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                          <  2:247659/124‑2 (MQ=38)                                 aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                          <  2:243424/124‑2 (MQ=38)                                 aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                          >  1:247659/1‑123 (MQ=38)                                 aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                          >  1:243424/1‑123 (MQ=38)                                    ttCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCATTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTgg                                                                                                              >  4:62508/1‑149 (MQ=255)                                         cAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAAc                                                                                                          >  5:187549/1‑148 (MQ=255)                                         cAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACg                                                                                                         <  6:187549/149‑1 (MQ=255)                                               gTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTAtt                                                                                                   >  2:289838/1‑149 (MQ=255)                                                        cGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCCTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTAttt                                                                                          <  3:62508/149‑1 (MQ=255)                                                                      cTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATa                                                                                                                                          <  8:18824/87‑1 (MQ=255)                                                                      cTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCACTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATa                                                                                                                                          >  7:18824/1‑87 (MQ=255)                                                                           aCTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATtcctc                                                                       <  5:48046/149‑1 (MQ=255)                                                                             tCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCg                                                                     >  4:171106/1‑149 (MQ=255)                                                                               gggTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCG‑GGATTGAAATtcctcc                                                                      <  8:137252/145‑1 (MQ=255)                                                                               gggTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCG‑GGATTGAAATtcctcc                                                                      >  7:137252/1‑145 (MQ=255)                                                                                       gTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCAcgcg                                                            >  2:64658/1‑148 (MQ=255)                                                                                            aCGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTc                                                      >  1:232523/1‑149 (MQ=255)                                                                                               aGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGa                                                   <  1:64658/149‑1 (MQ=255)                                                                                                            tttCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGGGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGt                                       >  7:302921/1‑148 (MQ=255)                                                                                                                  ccccTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGacac                                                >  5:74546/1‑133 (MQ=255)                                                                                                                  ccccTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGacac                                                <  6:74546/133‑1 (MQ=255)                                                                                                                  ccccTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAg                                  <  1:289838/147‑1 (MQ=255)                                                                                                                  ccccTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGc                                 >  4:178063/1‑148 (MQ=255)                                                                                                                    ccTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTAt                              >  4:120646/1‑149 (MQ=255)                                                                                                                       tcgtaTGCTGGAGAACCGCCATGG‑GCTACTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCg                            >  4:98589/5‑148 (MQ=38)                                                                                                                        gCAATGCTGGAGTAGC‑CCATGGTGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCgaa                          <  8:302921/149‑1 (MQ=255)                                                                                                                                 gAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTcc                 >  5:180806/1‑149 (MQ=255)                                                                                                                                 gAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTc                  >  7:198548/1‑148 (MQ=255)                                                                                                                                       c‑ccATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGtc           >  4:191959/1‑149 (MQ=255)                                                                                                                                            tGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTCTCGa       <  3:178063/149‑1 (MQ=255)                                                                                                                                                 gTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTCTCGAatcat  <  3:191959/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTCTCGAATCAT  >  NC_002947/2218224‑2218515

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