Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,462	C→G	8.9%	intergenic (+32/+45)	thiO → / ← PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	718,462	0	C	G	8.9%	101.6 / 2.7	44	intergenic (+32/+45)	thiO/PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base G (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACC  >  NC_002947/718315‑718609                                                                                                                                                     |                                                                                                                                                     aTCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCGCCGAagcg                                                                                                                                                     >  1:286545/1‑149 (MQ=255)  tCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                                      <  1:238548/148‑1 (MQ=255)     cAATGGGCTGGTGCTGGCGCCGGCTTCGTGGCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCg                                                                                                                                                 >  1:25131/1‑149 (MQ=255)      aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                                >  1:59971/1‑149 (MQ=255)      aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                                <  5:90960/149‑1 (MQ=255)          ggCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGtt                                                                                                                                             <  2:59971/148‑1 (MQ=255)          ggCTGGTGCTGGCGCCGGCTGCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGtt                                                                                                                                             <  5:266900/148‑1 (MQ=255)            cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAc                                                                                                                                          <  7:236197/149‑1 (MQ=39)               gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAcc                                                                                                                                         >  5:129350/1‑147 (MQ=39)               gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTAcc                                                                                                                                         <  6:129350/147‑1 (MQ=39)                     gcgcCGGCTTCGGGCCTGCTGTTTACCGTTTTCCTTCCCGGGGCCGAGCGTATAAGCGCCCCGGCGCCTTATG‑CCCCAGACGGGGGTTTGGGCTGAGGACCTCTTCGCGGTTAAACCCGCTCCAACGCGAGCGGGTTTACCCGCG‑AAAg                                                                                                                                <  2:255381/149‑1 (MQ=255)                       gcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAAACGACACGACACGAGGGGGTTTTCCCGAG‑aaaacg                                                                                                                              >  8:289717/1‑146 (MQ=255)                                gTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCAGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGaa                                                                                                                      <  4:167338/148‑1 (MQ=255)                                 tGCCTGCTTTTTACCGATTTGCTGCCCGGGGCCGCGCGTATTAGCGACCCGGCGCGGTATG‑CGCCAGATGGGCGTTTGGGCTGAGCACCTCTTCGCGGTTAACCCCGCTCCTACGCGCGCGGTTTTACCCGCG‑AAAGGGCCGCAGAAGc                                                                                                                    <  2:286545/149‑1 (MQ=255)                                  gCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTCCTAAGCGAGCGGGTTTAACCGCG‑AAAGGGCCGCAGAAGc                                                                                                                    >  7:172468/1‑148 (MQ=255)                                           ttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGGTCCTACGCGAGCGGGTTTTCCCGCG‑CAACGGGCCCCGAAGCCAATACAggt                                                                                                          >  5:299278/1‑147 (MQ=255)                                           ttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCTACGCGCGCGGGTTTTACCGCG‑AAAGGGCCGCAGAAGCCAATACAGAt                                                                                                          >  7:5518/1‑149 (MQ=255)                                            tACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACATCTTCGCGGGTAAAACCGCGCCTACGCGGGCCGGGTTTCCCGCG‑AAAGTGCCGCAGAACCCAATAAAGAtt                                                                                                         >  7:145098/1‑149 (MQ=255)                                               ggATTTGCTGACAGGGGCCGAGCCTATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGa                                                                                                                       <  2:223579/131‑1 (MQ=255)                                                 aTTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCAACAAACCGAGCGGGTTAAACCCCG‑AAACGGCCCTATAAGACAATAAAGAtttactg                                                                                                    >  6:167922/1‑145 (MQ=255)                                                      cTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCGACGACGCGAGCGGGGTTTACCGGG‑AAAGGGGCGCAGAAGCCAATAAAAATTTGAGGCAGt                                                                                                >  4:258153/1‑148 (MQ=255)                                                          cAGGGTCCGAGCCGATTATCTACCGGGGGCCGTATG‑CGCAAGATGGTGGTTTGGGGTTAGGACCTCTTCTCGGGTAAACCCGCTCCTACGCTTGGGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCg                                                                                           <  8:172468/149‑1 (MQ=37)                                                             gggCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGCCAAAGCGAGCGGGTTTTCCCGCG‑AAAGGGCCCCCGAAGCCAATAAAAATTTAAGCCAGTAATAGGTg                                                                                        >  2:208743/1‑149 (MQ=255)                                                                       gaTTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTg                                                                              >  1:29273/1‑149 (MQ=255)                                                                                cccGGCGCGGTATG‑CGCCTGAAGGGGGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGTGCGGTTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCgg                                                                     <  8:5518/149‑1 (MQ=39)                                                                                      gcCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGTAAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                               <  3:258153/149‑1 (MQ=255)                                                                                       ccGTATG‑CGCCAGACGGGCGTTTGGGCTGAGCACCTCATCGCGGGAAAACCCGCCCCAACGCGCGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCaa                                                                                                                 <  5:9950/98‑1 (MQ=37)                                                                                         gTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGGCAAAACGAGCGGGGGTTACCGCG‑AAAAGGCCCCAGAAACCAAAAAAAATTTAAGGAAGGAATAGGTGCTTTCAGGTGGCCGCGGGGATAAtta                                                              >  8:290266/1‑145 (MQ=255)                                                                                         gTAAG‑CGCCAGAACGGAGTTTCGGCTGAGGACCTCTTCGCGGGTAAAACCACACCAACACCACCCGGGTTACAAGCG‑AAACGGCCGCAGAAGCCAAAACAGAAATCAACCAGAAATCAATACTTTCATGTTCCGACTAGGATTAtgctg                                                            >  8:131978/1‑149 (MQ=255)                                                                                                      aGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCTACAAAGCGAGCGGGTTTACCCGCC‑AAAGGGGCGGAGAAGCCAAAACAGa                                                                                                           >  4:175835/1‑90 (MQ=255)                                                                                                      aGGGCGTTTGGGCTGAGGACCTCTTCGCGGGAAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGa                                                                                                           <  3:175835/90‑1 (MQ=255)                                                                                                       gggCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                               <  8:242458/149‑1 (MQ=255)                                                                                                       gggCGTTGGGGCGGCGCCCGTCTTCGCGGGTAAACCCGCTCCTACGCCAGGGGTTTTACCCGCG‑AAAGGGCCGTAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                               <  5:167922/149‑1 (MQ=255)                                                                                                         gCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCACCAACGCGAGCGGGGTTTCCCCCGCAAAGGGCC‑CAGAAACCCAAAAAGAATTAAGCCACTAATACGTTCTTTCAAGTGGCGGCGGGGATAAAACAGGTCAAACAccgatg                                              >  4:279323/1‑144 (MQ=255)                                                                                                                  cTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgcag                                    <  2:29273/149‑1 (MQ=255)                                                                                                                         ccTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTAACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                             <  1:208743/149‑1 (MQ=255)                                                                                                                         ccTCTTCGCGGGTAAACCCGCGCCTACGCGTGCGGTTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTg                              <  7:131978/148‑1 (MQ=255)                                                                                                                            tttCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtca                          >  6:10614/2‑149 (MQ=38)                                                                                                                             ttCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAAAACAGATTTCAGCCAGTAATCGGTGCTTGCACGTGCCGGCACGGATAATGCTGCTCATACACCTTGCACCCCCGCACCACTTGCtca                          >  4:260497/1‑148 (MQ=39)                                                                                                                                gcgGGTAAACCCTCTCTTACGCGTGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                       <  7:289717/148‑1 (MQ=255)                                                                                                                                      aaaCCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAgc                >  2:274518/1‑149 (MQ=255)                                                                                                                                        aCCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCTGCAAAgcgc              >  3:163783/1‑149 (MQ=255)                                                                                                                                        aCCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCACAgcgc              >  4:184557/1‑149 (MQ=255)                                                                                                                                          ccGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc            >  3:226565/1‑149 (MQ=255)                                                                                                                                           cGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCg           <  1:274518/149‑1 (MQ=255)                                                                                                                                             cTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGc         >  3:3704/1‑149 (MQ=255)                                                                                                                                                gaccgcgcgcgTTTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGcc        <  3:300775/144‑1 (MQ=255)                                                                                                                                                   cccgAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCATCAAAGCGCCCGGCcaccac   <  7:290266/147‑1 (MQ=255)                                                                                                                                                    gcgAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCaccacc  >  3:57016/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATG‑CGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTTACCCGCG‑AAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACC  >  NC_002947/718315‑718609

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