Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,457	C→G	8.0%	intergenic (+27/+50)	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,457	0	C	G	8.0%	124.6 / 4.1	50	intergenic (+27/+50)	thiO/PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (17/28);  new base G (2/2);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

CATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC  >  NC_002947/718311‑718605                                                                                                                                                   |                                                                                                                                                       cATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCGccaa                                                                                                                                                      >  8:315460/1‑147 (MQ=255)    tATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCACCTccgaacc                                                                                                                                                   >  4:228175/1‑144 (MQ=255)     aTCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                                   >  2:296376/1‑148 (MQ=255)      tCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc                                                                                                                                                   >  7:65360/1‑147 (MQ=255)       cGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAg                                                                                                                                               >  1:82666/1‑149 (MQ=255)          aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCgg                                                                                                                                            <  3:135674/149‑1 (MQ=255)              ggCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAACCCGCGCCGAAGCGGGGCGGGtt                                                                                                                                         >  7:142363/1‑149 (MQ=255)               gCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTa                                                                                                                                       <  1:302548/149‑1 (MQ=255)                cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑c                                                                                                                                     <  1:296376/149‑1 (MQ=255)                   gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑ccc                                                                                                                                   <  5:59977/148‑1 (MQ=255)                                cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt                                                                                                                                        <  4:106776/131‑1 (MQ=255)                                cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt                                                                                                                                        >  3:106776/1‑131 (MQ=255)                                              tttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCGGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa                                                                                                       <  7:315460/149‑1 (MQ=255)                                              tttACCGATTTGCTGACAGGCGCCGACCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTCCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa                                                                                                       <  1:123977/149‑1 (MQ=39)                                                   cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGTACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  7:28754/149‑1 (MQ=255)                                                   cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGTTTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  1:95891/149‑1 (MQ=255)                                                   cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  5:186243/149‑1 (MQ=255)                                                   cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  8:65360/149‑1 (MQ=255)                                                   cGATTTGCTCACAGGGGCCCAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGGGTAACCCCGCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  7:195812/149‑1 (MQ=255)                                                    gATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGGTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGAtttacg                                                                                                 >  6:183808/1‑146 (MQ=255)                                                               aGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCGCCTACGC‑GAGCGGGTTTT‑CCCGCGAAAGGGCCGCAGTAGCCAATAAAGATTTAATCCAGTAATCgg                                                                                      >  8:75121/1‑149 (MQ=255)                                                                ggggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCACCTACGC‑GAGCGGGGTTA‑CCCCCGAAAAGGCCGCACAACCCAAAAAAGATTGAAGCCAGTAAACGGt                                                                                     >  3:75685/1‑149 (MQ=255)                                                                  ggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTACAACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAAAAGGTGc                                                                                   >  6:144671/1‑149 (MQ=255)                                                                       ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c                                                                                                                                     <  2:89735/93‑1 (MQ=25)                                                                       ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c                                                                                                                                     >  1:89735/2‑94 (MQ=25)                                                                          cgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGt                                                                           <  5:183808/149‑1 (MQ=255)                                                                            aTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGc                                                                         <  2:260830/149‑1 (MQ=255)                                                                                  gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATAAAGATTTGAGCCAGTAATAGGTGGTTGCAGGTGCCGGcgc                                                                   >  6:119289/1‑147 (MQ=255)                                                                                  gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGCCGAAGC‑GAGCGGGGTTT‑ACCGCGAAAGGGGCGCCGAAGGCCATAAAGAGTTCAGCCAGTAATAGGTGGTTGCAGGTggcggc                                                                     >  4:79136/1‑147 (MQ=255)                                                                                     cgggggcCTTATGCGCCAGAAGGGCTTTTGGGCGGTGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑TAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGa                                                                <  3:228175/144‑1 (MQ=38)                                                                                      gggCGCCTTATCCCCCAGACGCGGGTTTGCGCGCAGCACGGCTTCGCGGTTAAACCCGCCCCACCGC‑GAGCGGTTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGAt                                                               <  3:29266/148‑1 (MQ=255)                                                                                          gcCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           <  7:247735/149‑1 (MQ=255)                                                                                              tATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCTACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAATAGGTGCTTGCAGGTGCCGGCTCGGCTAAtgctgc                                                       >  2:228366/1‑149 (MQ=255)                                                                                                     atactggggGTTTGGGCTGCGAACGCCTGCG‑GGGTAAAACCGCCACTACCC‑CCGGGGGTTTACCCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATaca                                                <  1:367396/144‑1 (MQ=255)                                                                                                          gggCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                           <  2:82666/149‑1 (MQ=255)                                                                                                                   ggCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgc                                  <  1:228366/149‑1 (MQ=255)                                                                                                                     cTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgcag                                <  5:144671/149‑1 (MQ=255)                                                                                                                            ccgctgcgcgGGTAAACCCGCCCCACCCC‑CAGCCGGATTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                         <  8:142363/143‑1 (MQ=255)                                                                                                                            ccTCTTCGCGGGTAACCCCGCTCCAACGC‑TCGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                         <  3:79136/149‑1 (MQ=255)                                                                                                                              tctTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa                                                                                                  <  3:125841/74‑1 (MQ=255)                                                                                                                                   gcgGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg                  >  3:211197/1‑149 (MQ=255)                                                                                                                                   gagtgaaaaCACGCTCATACGC‑CAGCGGTTTTA‑CCCGCGAAAGTGCCGCAGAATCCAATACATATTTCAGCCATTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg                  <  6:117778/143‑1 (MQ=255)                                                                                                                                      ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                   >  4:208397/1‑145 (MQ=39)                                                                                                                                      ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca                   <  3:208397/145‑1 (MQ=39)                                                                                                                                         aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           <  4:141095/102‑1 (MQ=255)                                                                                                                                         aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                           >  3:141095/1‑102 (MQ=255)                                                                                                                                           aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc                                          <  6:68743/117‑1 (MQ=255)                                                                                                                                           aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc                                          >  5:68743/1‑117 (MQ=255)                                                                                                                                             ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        >  6:227171/1‑149 (MQ=255)                                                                                                                                             ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        >  2:272339/1‑149 (MQ=255)                                                                                                                                             ccGCTCCTACCC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc        <  7:75121/149‑1 (MQ=255)                                                                                                                                                cTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGc     <  5:119289/149‑1 (MQ=255)                                                                                                                                                   cTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGccac  >  8:80549/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       CATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC  >  NC_002947/718311‑718605

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