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	84.6%	110.1 / 16.2	39	intergenic (‑99/+50)	tdcG‑II/gcvP‑I	L‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base C (14/19);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 9.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAA‑CC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126499‑1126784                                                                                                                                                         |                                                                                                                                             ggAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATAGCA‑TTAGC‑GGG‑T‑AAA‑cc                                                                                                                                            >  7:61222/1‑149 (MQ=255)       gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTAACCCCGCATGCGGAGGTATTTCCCCCAAAGCA‑TTAGC‑GGG‑G‑AAA‑cccgcgac                                                                                                                                      >  8:178241/1‑143 (MQ=255)        ccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTcc                                                                                                                                      >  2:22636/1‑148 (MQ=255)              ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTAcaca                                                                                                                                <  2:81176/148‑1 (MQ=255)                          aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCAcct                                                                                                                   >  4:56247/1‑148 (MQ=255)                            gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCaa                                                                                                                 >  5:357480/1‑149 (MQ=255)                             tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTTTAGCCACCATAGCA‑TTCGC‑GGG‑T‑AAACCC‑GCGCCTACACAGAGCCCCCCCCCCAAg                                                                                                                >  2:96873/1‑149 (MQ=38)                              cGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                               <  8:175726/149‑1 (MQ=255)                               gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTTAACCCGAAATAGGGGGTATATCCACCCACGCC‑ATAGCGGGG‑T‑AAA‑CC‑CCTCCAACACACAGcccccccccccaccc                                                                                                              >  6:12862/1‑144 (MQ=255)                               gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTTAACCCGAAATAGGAGGTATAGCCACCATCGCA‑TTAGC‑GGG‑G‑AAACCC‑CCTCCTAAACAGACCCCCCCCCCCAAAcc                                                                                                              >  3:133999/1‑149 (MQ=255)                                         aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCACTTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                      <  1:96873/148‑1 (MQ=255)                                             aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                 <  5:106933/148‑1 (MQ=255)                                             aCATGGGGACCTCCTGCCCGGTTGTTGGGTGCACA‑TTCCCCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTAAAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                 <  1:350710/148‑1 (MQ=37)                                              cATGGGCACCTCCTGCCAGGTTGTTGGGGGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTAAAGCCACCATCCCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                               <  5:12862/149‑1 (MQ=37)                                                    gACCTCCTGCCAGGTTGTTGGGTGGACATTTGAGCCATGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTATATCCACCATAGCA‑ATAGC‑CGG‑TAAAA‑CC‑CCGCCAAAACACAGCCCCCCACCCAAACCAGGTGGCGGGGAGGGCTAt                                                                                           >  1:211147/1‑148 (MQ=255)                                                        tcctGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTAACCCCGCATGCGGGGGTTTAGCCCCCATAGCA‑ATAGC‑GGG‑G‑AAACCC‑CCGCCAAAAAGGTACCCCCCCCCCCACCCCGGTGGGGGGGAGGGCGATGAGGCg                                                                                     >  1:124330/1‑149 (MQ=255)                                                          ctGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                     <  8:61222/147‑1 (MQ=255)                                                             ccAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt                                                                                 <  6:357480/148‑1 (MQ=255)                                                             ccAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATAGCA‑TTAGC‑GGGTT‑AAC‑CC‑GCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTTTCAGGCGGCCt                                                                                 >  5:301805/1‑148 (MQ=38)                                                                  ttgttgGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATATCCACCATCGCA‑TTCGC‑GGGTT‑AAA‑CC‑GCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATAAGGCGGACTGATAGc                                                                           >  6:80024/1‑149 (MQ=38)                                                                    gttgGGTGGACT‑TTGAGCCCTGTAGGCGGGGGTTTACCCGCTAATGCGGAGTTATAGCCACCAGCGCT‑TTCGG‑GGT‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGcc                                                                          <  1:174605/148‑2 (MQ=255)                                                                    gttgGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct                                                                          <  1:22636/148‑1 (MQ=255)                                                                              cA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                <  6:176801/148‑1 (MQ=255)                                                                                   gAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGc                                                            <  3:333203/148‑1 (MQ=255)                                                                                      cccTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                        <  7:240187/149‑1 (MQ=255)                                                                                               cgcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAg                                               >  2:307975/1‑149 (MQ=255)                                                                                               cgcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑CCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCCATCGACGGGCAGGCGCACACCAg                                               >  8:263377/1‑149 (MQ=255)                                                                                                    gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGc                                          <  7:178241/149‑1 (MQ=255)                                                                                                         cccGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                     >  3:52854/1‑149 (MQ=255)                                                                                                           cgcgAATGCGGAGGTATAGCCACCATCGCA‑TTCGG‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                   <  2:211147/149‑1 (MQ=255)                                                                                                           cgcgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                   <  2:124330/149‑1 (MQ=255)                                                                                                           cgcgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCAAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                   >  4:202160/1‑149 (MQ=255)                                                                                                           cgcgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAA‑CCAGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCACGTTGCGGACGcc                                   >  8:134065/1‑149 (MQ=255)                                                                                                            gcgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                  <  7:8289/149‑1 (MQ=255)                                                                                                             cgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                            >  1:54850/1‑74 (MQ=255)                                                                                                             cgAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                            <  2:54850/74‑1 (MQ=255)                                                                                                                      aGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCTCCGAACACGTTGt                        <  7:263377/149‑1 (MQ=255)                                                                                                                          ataGCCACCATCGCA‑TTCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                    >  8:291561/1‑149 (MQ=255)                                                                                                                                          ttCGC‑GGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg     >  2:75910/1‑149 (MQ=255)                                                                                                                                             gc‑gGG‑T‑AAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc    >  2:351630/1‑147 (MQ=255)                                                                                                                                             gc‑gGG‑T‑AAA‑CCTGCTCCTACACATAGCCCGCCACCCAAGCGAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAGGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa  <  7:134065/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                             GGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACA‑TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCA‑TTCGC‑GGG‑T‑AAA‑CC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126499‑1126784

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