Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,199,963:1	+C	7.4%	coding (829/1278 nt)	PP_1944 →	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,199,963	1	.	C	7.4%	208.5 / 9.0	54	coding (829/1278 nt)	PP_1944	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (26/24);  new base C (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199821‑2200111                                                                                                                                                |                                                                                                                                                     tGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgc                                                                                                                                                >  3:263079/1‑149 (MQ=255)    cTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                             >  1:147536/1‑149 (MQ=255)    cTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                             <  1:183060/149‑1 (MQ=255)    cTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                             <  6:95839/149‑1 (MQ=255)        aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctg                                                                                                                                         >  8:23883/1‑149 (MQ=255)        aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctg                                                                                                                                         >  2:162469/1‑149 (MQ=255)         tGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgggctgg                                                                                                                                        <  1:309756/149‑1 (MQ=255)           cccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctgcga                                                                                                                                      <  1:162469/149‑4 (MQ=255)           cccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctgcga                                                                                                                                      <  7:23883/149‑4 (MQ=255)            ccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCtgggctgggc                                                                                                                                      >  3:105423/1‑148 (MQ=255)             cgcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGcc                                                                                                                                     >  7:185322/1‑148 (MQ=255)              gcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCa                                                                                                                                    >  2:132196/1‑148 (MQ=255)                     tCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCAt                                                                                                                             <  2:85989/148‑1 (MQ=255)                        aCCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                            <  4:180333/146‑1 (MQ=255)                        aCCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                            >  3:180333/1‑146 (MQ=255)                                attggtggTTATCTCAAGGAATTCCGCGCGGGTATGCGAAGTTTTGGCTCGACTTTCTATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGTTATCGCAGCCCAG‑TCGAGGTGGGCTGGGCCTATACCATCAAGTTCGATc                                                                                                                  <  2:304129/142‑1 (MQ=255)                                aTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCa                                                                                                                 <  5:188820/149‑1 (MQ=255)                                           tctcAAGGAATGCCGGGCGCGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGCAGCTTTGAGGCGAATAAGGTCAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATc                                                                                                       <  4:267073/148‑1 (MQ=255)                                               aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                   <  4:105423/148‑1 (MQ=255)                                                  gAATTCCGCGCGGGTATGCCAAGTTTTGCCACGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCaa                                                                                                                          <  4:271918/122‑1 (MQ=255)                                                  gAATTCCGCGCGGGTATGCCAAGTTTTGCCACGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCaa                                                                                                                          >  3:271918/1‑122 (MQ=255)                                                          cgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTg                                                                                        >  3:97786/1‑148 (MQ=255)                                                             gggTATGCCGAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                   >  8:291871/1‑134 (MQ=255)                                                             gggTATGCCGAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTACCGGTc                                                                                                   <  7:291871/134‑1 (MQ=255)                                                             gggTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaa                                                                                     >  2:151105/1‑148 (MQ=255)                                                               gTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         >  5:216255/1‑142 (MQ=255)                                                               gTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         <  6:216255/142‑1 (MQ=255)                                                                   gCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAgg                                                                              >  8:253703/1‑149 (MQ=255)                                                                     cAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGtt                                                                            >  1:3656/1‑149 (MQ=255)                                                                       aGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGc                                                                          <  5:57980/149‑1 (MQ=255)                                                                         ttttGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCg                                                                                                      <  6:322113/119‑1 (MQ=255)                                                                         ttttGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTg                                                                           <  4:242451/146‑1 (MQ=255)                                                                         ttttGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTg                                                                           >  3:242451/1‑146 (MQ=255)                                                                             gCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATAACGCCTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGcc                                                                    >  4:56793/1‑149 (MQ=255)                                                                              ccTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGcc                                                                    <  4:263079/148‑1 (MQ=255)                                                                                      ttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATAGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                           >  7:39091/1‑149 (MQ=255)                                                                                           tGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCg                                                                                                      <  1:102466/101‑1 (MQ=255)                                                                                           tGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCg                                                                                                      >  2:102466/1‑101 (MQ=255)                                                                                                          tGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  1:151105/147‑1 (MQ=255)                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGaa                                      <  7:253703/148‑1 (MQ=255)                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAAt                                     >  6:38039/1‑149 (MQ=255)                                                                                                                gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCGGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGt                                  >  3:21410/1‑148 (MQ=255)                                                                                                                gAATAAGGT‑AGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtga                                >  8:65977/1‑149 (MQ=255)                                                                                                                     aGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGAc                            <  2:147536/149‑1 (MQ=255)                                                                                                                      ggTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACg                           >  4:230712/1‑149 (MQ=255)                                                                                                                        tAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGt                          <  3:230712/148‑1 (MQ=255)                                                                                                                          aGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGGTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGc                       <  4:97786/149‑1 (MQ=255)                                                                                                                            cGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGa                                                                >  7:287470/1‑106 (MQ=255)                                                                                                                            cGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGa                                                                <  8:287470/106‑1 (MQ=255)                                                                                                                               ttGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATg                  <  2:3656/149‑1 (MQ=255)                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGtt                >  2:51894/1‑148 (MQ=255)                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGtt                >  2:312847/1‑148 (MQ=255)                                                                                                                                      cGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc           >  4:64468/1‑149 (MQ=255)                                                                                                                                              aG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGt   >  6:158417/1‑149 (MQ=255)                                                                                                                                               g‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGtt  >  4:106930/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     TGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199821‑2200111

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