Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,085	A→G	15.2%	intergenic (+53/‑77)	tdcG‑I → / → PP_5441	L‑serine dehydratase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	358,085	0	A	G	15.2%	98.9 / 10.0	47	intergenic (+53/‑77)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base G (5/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357943‑358229                                                                                                                                                |                                                                                                                                                  cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACtgtg                                                                                                                                             <  7:37273/149‑1 (MQ=255) cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACtgtg                                                                                                                                             <  2:59903/149‑1 (MQ=255) cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACtgt                                                                                                                                              <  3:39328/148‑1 (MQ=255)          cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                     <  4:339129/148‑1 (MQ=255)           gcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGt                                                                                                                                   <  3:75422/149‑1 (MQ=255)            cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                  <  7:90926/149‑1 (MQ=255)               tACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGa                                                                                                                               <  8:200412/149‑1 (MQ=37)                  cGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTACTACAGGAACTGTGCAATCCCTGTAGGAGCg                                                                                                                            >  7:170736/1‑149 (MQ=255)                   ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgg                                                                                                                           <  5:314489/149‑1 (MQ=255)                   ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgg                                                                                                                           <  2:42162/149‑1 (MQ=255)                    gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGCGCggg                                                                                                                          >  5:308219/1‑149 (MQ=255)                    gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                          <  6:10015/149‑1 (MQ=255)                            atGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCcg                                                                                                                  >  2:271877/1‑149 (MQ=255)                                     aaGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGGAGGAGCGGGCTTGCGCGCGAAGAGGc                                                                                                         >  6:322035/1‑149 (MQ=255)                                             agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                                                                                 >  8:88734/1‑149 (MQ=255)                                                  ccTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GAt                                                                                           <  1:97576/149‑1 (MQ=255)                                                   cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑CAGCCCGCGCCTACAGGAACTGTGGGATACCTGTTGGCGCGGGCTTGGCCGCGCAGAGGGCCGTGGGGGG‑GATg                                                                                          >  1:317081/1‑149 (MQ=255)                                                    tCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACCGGGACTGTGCAATGCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGGT‑GATGc                                                                                         >  7:190501/1‑149 (MQ=255)                                                      gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑CAACCCGCGACGAAACGACAAGTGGAATACCAGTAGGCGCGGGGGTTCCCGCGAACAGGCCCCTGCGGCG‑Gagtcag                                                                                       >  8:127589/1‑144 (MQ=255)                                                      gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAACCCCCTACTACACGCACAGTGGCAAACCGGTAGGCGAGGGGGTGCACCCGCACAGGCCCGTGAGGCT‑GAGGCAg                                                                                       >  8:123004/1‑149 (MQ=255)                                                        gcggcgGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGt                                                                                     >  8:65420/1‑149 (MQ=255)                                                                   gCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAAGCCCGCTCCTACAGGAACTGTGCAATCCGGGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑TATGCAGGTGCAACACAAc                                                                           >  1:271092/1‑149 (MQ=38)                                                                     ggTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGc                                                                        <  7:65420/149‑1 (MQ=255)                                                                         aGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTAt                                                                    >  2:32621/1‑149 (MQ=255)                                                                           cgcTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGGAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGc                                                                  >  8:171272/1‑149 (MQ=255)                                                                           cgcTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGc                                                                  <  7:88734/149‑1 (MQ=255)                                                                                gcGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑TAGCCCGCTCGTACAGCAACTGTGCAAACCCTGTAGGAGCGGGGTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCctct                                                              <  5:322035/147‑1 (MQ=255)                                                                                    gTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATg                                                         <  1:32621/149‑1 (MQ=255)                                                                                                  aGCTTACCTTCTGGGCCTCTTCGCGGGA‑AAGCCCGCTCCTACTGTAACTGTGCAATCCCTTTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATgg                                             <  6:308219/147‑1 (MQ=38)                                                                                                         cTTCTGGGCCTCTTCGCGGGC‑AAGCCCCCTACGTCACGAACTGTGCCATCCCAGGAGGAGCCGGGTTGGCCGCGAAGAGGGCCTTGAGGGT‑GATGGAGGGGGCAAACAACCGCCTATGGCGAGTTATGGTGGGGGGCTGGGGGcccccc                                    >  2:277544/1‑147 (MQ=255)                                                                                                               ggccctttCGCGGCC‑AAGCCCGCCCCAACCGCACCTGTGCAACCCCTGACGCCGGGGGGTTGCCGGCTATTAGGCCCTTGAGGCTCCATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTa                                 <  1:271877/141‑1 (MQ=255)                                                                                                                 ccTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGAAATcc                                                                                                                                       >  3:309114/1‑43 (MQ=255)                                                                                                                             c‑AAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATgg                                             >  7:52360/1‑120 (MQ=21)                                                                                                                             c‑AAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATgg                                             <  8:52360/120‑1 (MQ=21)                                                                                                                             c‑AAGCCCGCTCCTACAGCAACTTTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTcg                <  8:170736/149‑1 (MQ=255)                                                                                                                                 gCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATgg                                             <  3:237748/117‑1 (MQ=255)                                                                                                                                 gCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATgg                                             >  4:237748/1‑117 (MQ=255)                                                                                                                                 gCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgtc                        >  4:3916/1‑138 (MQ=255)                                                                                                                                 gCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgtc                        <  3:3916/138‑1 (MQ=255)                                                                                                                                 gCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGa             <  7:171272/149‑1 (MQ=255)                                                                                                                                  cccGGTCCTACAGTACGTGTGCAATCGCTGTAGTACGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt            <  7:128184/149‑1 (MQ=255)                                                                                                                                    cGCACCTACAGAAACGGGGAACCCCCTGATGGTGGGGGCTTCCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa          <  5:294606/149‑1 (MQ=255)                                                                                                                                       tCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg       >  6:301538/1‑149 (MQ=37)                                                                                                                                       tCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg       >  8:257955/1‑149 (MQ=37)                                                                                                                                       tCCTACAGGGATTGCACAGTTCCTCTAGGAGCGGGCATGCCCGCGAAGAGGCCCTTGACTCT‑GATGCAGGTGCAACACAACCGCCAATGCCTCTTAATGGTGCGGTGATGGCCGCCCCGCCAACAGATCGTCGGATGTCGCGATCAgtg       >  6:340847/1‑149 (MQ=21)                                                                                                                                          aaccgcaaCGGTCCACTCCCTGTAGGCGCGGGCGTCCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg    <  3:185086/143‑1 (MQ=255)                                                                                                                                           aCAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTggg   >  8:308689/1‑149 (MQ=255)                                                                                                                                            cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTACAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc  >  2:175881/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                  CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGC‑AAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCT‑GATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357943‑358229

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑