Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,126,647	+C	intergenic (‑101/+48)	tdcG‑II ← / ← gcvP‑I	L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,126,645	1	.	C	95.6%	147.1 / ‑0.2	45	intergenic (‑99/+50)	tdcG‑II/gcvP‑I	L‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (18/25);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126497‑1126788                                                                                                                                                        |                                                                                                                                                cTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTaaa                                                                                                                                                  >  1:65234/1‑149 (MQ=255)                 cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTAcaca                                                                                                                                    >  1:253789/1‑147 (MQ=255)                           cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTATAGCCA‑CCATAGCATTAGCGG‑GTAAACCCGCGACTACACAGAGCcccccccc                                                                                                                        >  8:116153/1‑149 (MQ=255)                            aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCCCATGAGGAGGTATATACACCCAACCAATAACCGG‑GTAAA‑CCCCAGCAACAAAAAGCCCCCCAccc                                                                                                                       >  6:289366/1‑149 (MQ=255)                              gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                     <  1:17141/149‑1 (MQ=255)                              gTCGAAGACGCTCAGTGACATGTGCACCTCCTGCCAGGTTGTTGGGTGGACATTGACCCCTGTCGGCGCGGGTTTCCCCGCGAAGGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCAc                                                                                                                         <  5:54340/145‑1 (MQ=38)                                               aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                     <  3:274391/148‑1 (MQ=255)                                                 aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                     <  7:163720/130‑1 (MQ=255)                                                   gtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGt                                                                                                              <  7:258404/135‑1 (MQ=255)                                                   gtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAAAGCGGAGGTATAGCCA‑CCATCGCATTAGCGG‑GTAAACCCGCACCAACCCAGAGCCCCCCACCCAAGCCAggg                                                                                                              >  8:258404/1‑134 (MQ=255)                                                         ctcctGCCCGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGAATAGCCA‑CAATCGCATTCGCGG‑GAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                          <  5:289366/149‑1 (MQ=255)                                                          tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACGCGAATGCGGAGGTATAAACA‑ACCAAGCAATAGCGGTGTTAA‑CAGCTCATACACAAAGACCCCCCACCAAACCAAGGGGGGGGGAGGGCTTTAAGGCg                                                                                         >  2:120681/1‑149 (MQ=255)                                                                cAGGTTGTTGGGTTGACATTGAGCCCTGTAGGCGTGGGTTAACCCTCGAATGCGGAGGTATAGCCA‑CCATCNCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                   <  1:120681/149‑1 (MQ=38)                                                                 aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                                  <  2:65234/149‑1 (MQ=255)                                                                  ggTTGTTGGGTGGACATTGAGCCCTGTAGGCGCAGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                   >  2:334325/1‑147 (MQ=255)                                                                   gttgttGGGTGGACATTGCGCCCTGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                                <  4:191259/149‑1 (MQ=255)                                                                      gttgGGTGGACATTGAGCCCCGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAGcc                                                                              <  5:197636/148‑2 (MQ=38)                                                                       ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                            <  7:116153/149‑1 (MQ=255)                                                                          ggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑TTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCAATCAGGCGTCCTGATAGCTCtcga                                                                         <  1:202257/149‑1 (MQ=255)                                                                           gggcaCATTCCGCCCTGTCGGCGCGGGTTAACCCGCTATTGCGTCGTTAACCCAA‑CCATCCCATTCGCGG‑TAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat                                                                        <  2:380135/145‑1 (MQ=255)                                                                               aCATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCGGATAGCTGTCGATCGAc                                                                    >  1:48645/1‑149 (MQ=38)                                                                                cATTGCGCCGTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg                                                                   <  8:323909/149‑1 (MQ=255)                                                                                               tcgcgGGTAAACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  >  7:81803/2‑135 (MQ=255)                                                                                               tcgcgGGTAAACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  <  8:81803/134‑1 (MQ=255)                                                                                                      tttACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                             <  1:334325/149‑1 (MQ=255)                                                                                                      tttACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                             <  3:77757/149‑1 (MQ=255)                                                                                                      tttACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                             >  4:77757/1‑149 (MQ=255)                                                                                                          cccGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                         >  6:313550/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                             <  4:188559/96‑1 (MQ=255)                                                                                                             gcgAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCTACgg                                                                  <  8:73571/121‑1 (MQ=255)                                                                                                             gcgAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  >  7:73571/1‑121 (MQ=255)                                                                                                                     ggAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                   <  1:146499/144‑1 (MQ=255)                                                                                                                     ggAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                   >  2:146499/1‑144 (MQ=255)                                                                                                                              gCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                     >  8:344675/1‑149 (MQ=255)                                                                                                                              gCCA‑CCATCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                     >  5:145404/1‑149 (MQ=255)                                                                                                                               cca‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                            <  5:396546/93‑1 (MQ=255)                                                                                                                               cca‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                            >  6:396546/1‑93 (MQ=255)                                                                                                                               cca‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  5:310083/1‑149 (MQ=255)                                                                                                                                ca‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                   <  7:241286/149‑1 (MQ=255)                                                                                                                                 a‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            <  4:131191/139‑1 (MQ=255)                                                                                                                                 a‑ccaTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            >  3:131191/1‑139 (MQ=255)                                                                                                                                     aTCGCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc               >  1:143465/1‑149 (MQ=255)                                                                                                                                        gCATTCGCGG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcgg            <  6:310083/149‑1 (MQ=255)                                                                                                                                             cgcgG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        <  5:313550/148‑1 (MQ=255)                                                                                                                                              gcgG‑GTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGACTATCAGGCGTCCTGATAGCTCTCGAACGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        >  7:200847/1‑147 (MQ=255)                                                                                                                                                   gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  4:315722/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCA‑CCATCGCATTCGCGG‑GTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126497‑1126788

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