Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,086,043	G→A	8.1%	A248T (GCC→ACC)	amaD →	D‑lysine oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,086,043	0	G	A	8.1%	142.1 / 2.2	50	A248T (GCC→ACC)	amaD	D‑lysine oxidase
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (2/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

CCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGATCTGCGACACCCAAGGCGGCTATG  >  NC_002947/4085895‑4086190                                                                                                                                                      |                                                                                                                                                    ccccGGCGCGCTCACCCGTGTCTATGCGGCGTTGTGCCTGCAAGGGGGCGGGCAGTTGGTGCATGGCGATGGGCGCAGCGTGCGCCTGGCCAATGGCCAATGGCGGGTAGAGAGGCGT‑CGCGGCGCGATCACTGCCGACGAGGttgtgg                                                                                                                                                     <  6:314092/149‑1 (MQ=255)  cccGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGgc                                                                                                                                                    >  6:332714/1‑149 (MQ=255)     ggCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGgcctg                                                                                                                                                 >  6:178281/1‑149 (MQ=255)        gcgcTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGgcctgcct                                                                                                                                              >  1:87292/1‑149 (MQ=255)           cTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTgg                                                                                                                                           >  3:196301/1‑149 (MQ=255)              aCCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTggg                                                                                                                                          <  3:350921/147‑1 (MQ=255)              aCCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTggg                                                                                                                                          <  3:345626/147‑1 (MQ=255)              aCCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTggg                                                                                                                                          <  8:29824/147‑1 (MQ=255)                  gTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAg                                                                                                                                    >  1:84212/1‑149 (MQ=255)                     gCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCa                                                                                                                                 <  4:196301/149‑1 (MQ=255)                       tATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCTGAGAATGCCGCCCAGTTGGTGGCGTGCCTTTGGGCCCAGACAGc                                                                                                                               >  1:309322/1‑149 (MQ=255)                       tATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTgg                                                                                                                                           >  8:127174/1‑137 (MQ=255)                       tATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTgg                                                                                                                                           <  7:127174/137‑1 (MQ=255)                       tATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGc                                                                                                                               >  8:157120/1‑149 (MQ=255)                       tATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGc                                                                                                                               >  3:207002/1‑149 (MQ=255)                           cggcgTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGTCCTGCCTTGGGCCACAGTCAGCCGAc                                                                                                                           <  5:332714/149‑1 (MQ=255)                             gcgTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCt                                                                                                                         <  7:157120/149‑1 (MQ=255)                                ttgttCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGtt                                                                                                                      <  1:178983/149‑1 (MQ=255)                                        gCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGcc                                                                                                              >  1:6399/1‑149 (MQ=255)                                                 cggGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTAc                                                                                                      >  2:196081/1‑148 (MQ=255)                                                  gggCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCa                                                                                                    <  4:207002/149‑1 (MQ=255)                                                     cAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGAt                                                                                                 <  5:309799/149‑1 (MQ=255)                                                           gTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGc                                                                                            <  2:84212/148‑1 (MQ=255)                                                           gTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGc                                                                                            <  2:87292/148‑1 (MQ=255)                                                           gTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCt                                                                                           >  7:272541/1‑149 (MQ=255)                                                              cATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTgg                                                                                         <  1:196081/148‑1 (MQ=255)                                                               aTGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGc                                                                                        >  6:321515/1‑148 (MQ=255)                                                                  gCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAAt                                                                                     >  1:126725/1‑148 (MQ=255)                                                                    gATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATc                                                                                    >  3:341295/1‑147 (MQ=255)                                                                        cgcgcAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGTGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCgg                                                                              <  5:321515/149‑1 (MQ=255)                                                                         gcgcAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTAc                                                                                                      <  4:267662/124‑1 (MQ=255)                                                                         gcgcAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTAc                                                                                                      >  3:267662/1‑124 (MQ=255)                                                                                 tGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACca                                                                      >  7:69093/1‑148 (MQ=255)                                                                                    gcCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGTCCGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACAt                                                                   >  1:55549/1‑149 (MQ=255)                                                                                     ccAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGc                                                                 <  2:126725/149‑1 (MQ=255)                                                                                       aGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCAc                                                               >  1:7252/1‑149 (MQ=255)                                                                                           caatggccaatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACa                                                           <  2:55549/149‑1 (MQ=255)                                                                                           caatggccaatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACa                                                           <  2:309322/149‑1 (MQ=255)                                                                                                 ccaatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGACCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTTGCAATCAAGCGGGGCTACCACATGCACTACAGCAccc                                                     >  1:172869/1‑149 (MQ=255)                                                                                                 ccaatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGACCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCAccc                                                     >  7:17215/1‑149 (MQ=255)                                                                                                   aatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGt                                                   <  8:69093/149‑1 (MQ=255)                                                                                                   aatggcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTACAGCGGCCTGGGCTACCAGATCCCGCTGCCAATCAAGAGGGGCTACCACATGCACAACAGCCCCCGt                                                   >  8:354031/1‑149 (MQ=255)                                                                                                       gcGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGAc                                                <  6:69297/148‑1 (MQ=255)                                                                                                        cGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGACCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTTGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACgg                                              <  2:172869/149‑1 (MQ=255)                                                                                                        cGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGACCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACgg                                              <  8:17215/149‑1 (MQ=255)                                                                                                           gTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGgcgc                                           >  5:215265/1‑149 (MQ=255)                                                                                                              gagagCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCa                                        <  4:341295/149‑1 (MQ=255)                                                                                                                gagCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCa                                                        <  2:138803/131‑1 (MQ=255)                                                                                                                gagCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCa                                                        >  1:138803/1‑131 (MQ=255)                                                                                                                 agCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCt                                     <  5:178281/149‑1 (MQ=255)                                                                                                                    cGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGa                                  <  8:272541/149‑1 (MQ=255)                                                                                                                              ccGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGAt                         >  5:363281/1‑149 (MQ=255)                                                                                                                                gATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGATCt                       <  6:215265/149‑1 (MQ=255)                                                                                                                                             gaGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGATCTGCGACACCCAAgg          <  2:6399/149‑1 (MQ=255)                                                                                                                                                      gcctgcctTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGATCTGCGACACCCAAGGCGGCTATg  >  5:56478/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCGT‑CGCGGCCCGATCACTGCCGACGAGGTGGTGGCCTGCCTTGGGCCACAGTCAGCCGACCTGTTCAGCGGCCTGGGCTACCAGATCCCGCTGGCAATCAAGCGGGGCTACCACATGCACTACAGCACCCGTGACGGCGCTCAGCTGGAGCACTCGATCTGCGACACCCAAGGCGGCTATG  >  NC_002947/4085895‑4086190

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