Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,218,358:1	+G	13.5%	intergenic (‑184/+56)	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,358	1	.	G	13.5%	124.7 / 12.3	37	intergenic (‑184/+56)	PP_1958/PP_5495	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (16/16);  new base G (3/2);  total (19/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 = 8.96e-01

TATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTC  >  NC_002947/2218218‑2218506                                                                                                                                              |                                                                                                                                                     tATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCg                                                                                                                                              <  6:157900/149‑1 (MQ=255)       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                                                                                                                                 >  1:247659/1‑123 (MQ=38)                                       aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                 <  2:243424/124‑2 (MQ=38)                                       aagTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGAACCGCCATGG‑GCTACTCCAGCAtt                                                                                                                                 <  2:247659/124‑2 (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                                                             >  7:137252/1‑145 (MQ=255)                                                                                     gggTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCG‑GGATTGAAATtcctcc                                                             <  8:137252/145‑1 (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)                                                                                                                                              |                                                                                                                                                     TATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGC‑CCATGGCGGTTCTCCAGCATACGACCGAAATCCATCTACACCAGGCTGGTAACGCTTATTGGCGTATTTGCGAGGATTGAAATTCCTCCGCTTCACGCGTAATTCCGACACACTCGGAGTCAAAGCTATCGAAGAAGGCTCCGAAGTC  >  NC_002947/2218218‑2218506

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