Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,124,454:1	+C	12.5%	coding (689/1122 nt)	gcvT‑I ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,124,454	1	.	C	12.5%	104.8 / 10.1	32	coding (689/1122 nt)	gcvT‑I	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (10/18);  new base C (2/2);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAG  >  NC_002947/1124307‑1124555                                                                                                                                                      |                                                                                                      aaGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑g                                                                                                      <  4:244963/149‑1 (MQ=255)             gCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTgc                                                                                          >  3:232729/1‑149 (MQ=255)               aCCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCACGCGACTCGC‑Tgcgc                                                                                        >  3:315647/1‑149 (MQ=255)                cccGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCACGCGACTCGC‑TGCGcc                                                                                       <  4:315647/149‑1 (MQ=255)                   gggAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGccc                                                                                      >  2:287461/1‑147 (MQ=255)                   gggAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCaa                                                                                    <  1:266908/149‑1 (MQ=255)                   gggAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGACCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCaa                                                                                    <  5:261631/149‑1 (MQ=255)                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGc                                                                                  >  7:264750/1‑149 (MQ=255)                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGc                                                                                  <  3:160395/149‑1 (MQ=255)                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGc                                                                                  >  3:281219/1‑149 (MQ=255)                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGc                                                                                  <  4:281219/149‑1 (MQ=255)                         ccGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGAt                                                                              >  5:51331/1‑149 (MQ=255)                          cGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATc                                                                             >  7:169043/1‑149 (MQ=255)                          cGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATc                                                                             <  8:264750/149‑1 (MQ=255)                            gcagcaCGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCgg                                                                           <  1:260124/149‑1 (MQ=255)                               gcaCGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTg                                                                        <  3:162466/149‑1 (MQ=255)                                     gCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTc                                                                  <  7:19638/149‑1 (MQ=255)                                       aCCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCg                                                                <  2:40541/149‑1 (MQ=255)                                        ccGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCgg                                                               <  6:291577/149‑1 (MQ=255)                                          gTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGAGGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGAACTTCCgg                                                               >  1:308190/1‑148 (MQ=255)                                              gCACGGCGCACCTTGGAAATTTCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATAGGCTGTACTTCCGGCTCGGc                                                         <  2:308190/149‑1 (MQ=255)                                                          ttGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCACGCGACTCGC‑TGCGcct                                                                                      <  1:279839/108‑2 (MQ=255)                                                          ttGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCACGCGACTCGC‑TGCGcct                                                                                      >  2:279839/1‑107 (MQ=255)                                                                 ttGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGATCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAg                                      <  5:54220/149‑1 (MQ=255)                                                                     ccACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGt                                  <  8:169043/149‑1 (MQ=255)                                                                       acaGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCg                                <  1:287461/149‑1 (MQ=255)                                                                                    ttCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCg                   >  2:395139/1‑149 (MQ=255)                                                                                           gg‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCga            <  1:395139/149‑1 (MQ=255)                                                                                                 aGTCTCGTTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGACGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAggcggc                                              >  7:380756/1‑110 (MQ=255)                                                                                                 aGTCTCGTTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGACGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAggcggc                                              <  8:380756/110‑1 (MQ=255)                                                                                                     tcGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGa                      >  6:372540/1‑130 (MQ=255)                                                                                                     tcGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGa                      <  5:372540/130‑1 (MQ=255)                                                                                                       gCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAg  <  6:51331/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                      AAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGG‑CGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCA‑GCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAG  >  NC_002947/1124307‑1124555

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