Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,085	A→G	10.5%	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	10.5%	95.8 / 3.7	39	intergenic (+53/‑77)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/21);  new base G (2/2);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

ATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357940‑358229                                                                                                                                                    |                                                                                                                                                 aTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                <  8:248203/148‑1 (MQ=255) aTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                <  3:181977/148‑1 (MQ=255)    cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            <  1:98769/149‑1 (MQ=255)     gCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            >  4:330755/1‑148 (MQ=255)        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                                                                                                                                <  1:250469/108‑1 (MQ=255)                                                         gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAg                                                                                                                                >  2:250469/1‑108 (MQ=255)                                                           gcggcgGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                >  3:90443/1‑90 (MQ=255)                                                           gcggcgGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                <  4:90443/90‑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)                                                                                                                                                    |                                                                                                                                                 ATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357940‑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: ‑