Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,094	A→G	10.8%	intergenic (+62/‑68)	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,094	0	A	G	10.8%	93.9 / 3.8	37	intergenic (+62/‑68)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base G (2/2);  total (17/20)
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.87e-01

CCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357947‑358242                                                                                                                                                      |                                                                                                                                                     ccATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAAt                                                                                                                                                     <  4:388281/149‑1 (MQ=255)   aTGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATcc                                                                                                                                                   <  2:266330/149‑1 (MQ=255)      cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                 <  3:330755/148‑1 (MQ=255)       gcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGt                                                                                                                                               <  5:284623/149‑1 (MQ=255)        cgATACCGGCGCCTTTATGCACGACATGTACTAAGAGAGATCGCGAGGCGGGCTGGCTGTCAGCGCTATAGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                              <  6:377172/149‑1 (MQ=255)           tACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGa                                                                                                                                           <  5:345476/149‑1 (MQ=37)               ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGTCTCCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGc                                                                                                                                         <  1:117227/149‑1 (MQ=25)                 cgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgggt                                                                                                                                     <  4:4644/149‑2 (MQ=255)                 cgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGc                                                                                                                                     <  4:376321/149‑1 (MQ=255)                            aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGccc                                                                                                                               <  7:211655/144‑1 (MQ=255)                            aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                          <  7:308389/149‑1 (MQ=255)                            aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCCCTCCGACAGGCACTGTTCAATCCCTGGAGGATCGGGCTTGCCCgagat                                                                                                                          >  8:211655/1‑148 (MQ=38)                             ctactaGTACAATTAGACCTGGGGCGGGGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTTTGTGATG‑‑ACCTTGTGTGCCTGTTCGCGGGCAATCCCGCTCCTACAGTAAGTTTGAAATCCCTGTAGTTGCGGGTTTGCCCGCGAag                                                                                                                         <  6:217872/144‑1 (MQ=255)                                               cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTTCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                                       >  8:19429/1‑149 (MQ=255)                                                  gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAg                                                                                                                                             >  2:250469/1‑108 (MQ=255)                                                  gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAg                                                                                                                                             <  1:250469/108‑1 (MQ=255)                                                     cggcggCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTg                                                                                                 <  7:19429/149‑1 (MQ=255)                                                                                                              gCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGTAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATg                                                           <  5:140597/132‑1 (MQ=255)                                                                                                              gCCTCTTCGCGGGCAAGCCCCCAACTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCGATGCCACTTTATTGTGCGGTGATg                                                           >  6:140597/1‑132 (MQ=37)                                                                                                                 tctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCCATCCCAGTTGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGAAGCCGGTGCCACACAACCGCCGATGCCTCGTTATGGTGCGGTGAATGGCCACCCGCCTAACGAtcg                                       >  6:163709/1‑149 (MQ=25)                                                                                                                   ctCGCGGGCACGCCCGCTCATACAGGAACTGTGCAATCCCTTTAGGCGCGGGGTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgt                                      <  4:279868/147‑1 (MQ=255)                                                                                                                      gcgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgg                                      >  3:279868/1‑144 (MQ=255)                                                                                                                        gggCAAGCCCGCTCCTACAGGATCAGTGCAATCCCTGTAGGAGCGGGCTTGACCGCGAACAGGCCCTTGAGGATGATGAAGGAGAAACACAACCGCCAATGCCTCTTTATGGTGCGGTGATGGGCGCCCCCCATACCGAAAGTCGGAt                                 >  7:59014/1‑148 (MQ=37)                                                                                                                                ccGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGATGATGCAGGTGCAACACAACCGCCTATGAAAATTTATGGTGCGGTGATGGGCGCCCCGCAAACCGAACGACGGATGACGCCAt                         >  8:244558/1‑148 (MQ=18)                                                                                                                                  gCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGACCTTGAGTCTGATGCAGGTGCAACACAACCGACTATGCATCTTTATGGTGCGGTGATGGGAGCCCCTCCTACCGATCGTCGGATTTCGCGATAAg                      >  4:209851/1‑149 (MQ=38)                                                                                                                                   cgccTCCAGCACCTGTGCAATCCCTGTCGGCGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgt                     <  1:150384/147‑1 (MQ=255)                                                                                                                                    tCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg                    >  4:30000/1‑149 (MQ=255)                                                                                                                                      cTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg                  >  8:319172/1‑149 (MQ=37)                                                                                                                                      cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGCTCGTCGGATGTCGCGATCAGtgtg                  >  8:359053/1‑149 (MQ=255)                                                                                                                                       tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                 >  4:173793/1‑149 (MQ=255)                                                                                                                                         cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc               <  3:173793/149‑1 (MQ=255)                                                                                                                                            gAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                                 <  2:168873/128‑1 (MQ=255)                                                                                                                                            gAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                                 >  1:168873/1‑128 (MQ=255)                                                                                                                                                tgtgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGGACAGTGTGGGCGTTGTCg        >  1:135799/1‑149 (MQ=255)                                                                                                                                                  tgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  6:315130/1‑149 (MQ=255)                                                                                                                                                    cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  6:363944/1‑147 (MQ=255)                                                                                                                                                    cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      <  5:363944/147‑1 (MQ=255)                                                                                                                                                     aaTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                         <  1:31591/127‑1 (MQ=255)                                                                                                                                                     aaTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTCATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                         >  2:31591/1‑127 (MQ=255)                                                                                                                                                      aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGaca  >  2:133709/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357947‑358242

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