Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	66,400	Δ1 bp	7.3%	coding (1458/1653 nt)	betA‑I →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	66,400	0	G	.	7.3%	170.7 / 7.0	55	coding (1458/1653 nt)	betA‑I	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base . (2/2);  total (26/29)
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.97e-01

AGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCAAAGGA  >  NC_002947/66252‑66547                                                                                                                                                     |                                                                                                                                                    aGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATg                                                                                                                                                     <  4:227365/149‑1 (MQ=255)   gCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGa                                                                                                                                                   >  4:48174/1‑149 (MQ=255)     cccTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACg                                                                                                                                                 >  1:320188/1‑149 (MQ=255)       cTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGtt                                                                                                                                               <  8:284939/149‑1 (MQ=255)          cccTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGATGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgt                                                                                                                                             <  5:54944/148‑1 (MQ=255)          cccTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgt                                                                                                                                             <  2:207464/148‑1 (MQ=255)          cccTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCCGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtg                                                                                                                                              <  2:9784/147‑1 (MQ=255)            cTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgtg                                                                                                                                            <  4:182571/147‑1 (MQ=255)            cTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgtg                                                                                                                                            >  3:182571/1‑147 (MQ=255)                 cgcTTTCGAGCCCAACGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGAGCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGAGCAACGGTGCCATGGACGTTGTGGATAAcc                                                                                                                                     <  7:125106/149‑1 (MQ=255)                         aGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTgcgc                                                                                                                             <  7:15503/149‑1 (MQ=255)                          gCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgt                                                                                                                                             <  1:20701/132‑1 (MQ=255)                          gCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgt                                                                                                                                             >  2:20701/1‑132 (MQ=255)                            cgaaGGAATATTTGCCGGGCCCGTCGCTGGAAAGCGAGGAAGACCTGTTGGAGGCGGCTGGCAAGATCGGCACCACCATCTTCCAGCGGGTCGGTACCTGTCGGATGGGCAACTGTGCCATGGACGTTGTGGATAACCTGCTGCGCGTg                                                                                                                          <  7:175264/147‑1 (MQ=255)                              aaGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgtg                                                                                                                                            <  1:167380/129‑1 (MQ=255)                              aaGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTtgtg                                                                                                                                            >  2:167380/1‑129 (MQ=255)                                          ccGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCgg                                                                                                            >  6:218345/1‑149 (MQ=255)                                          ccGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCg                                                                                                             <  7:36787/148‑1 (MQ=255)                                                     gcTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAg                                                                                                 <  5:218345/149‑1 (MQ=255)                                                      cTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGc                                                                                                <  2:320188/149‑1 (MQ=255)                                                         cAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCCTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGcc                                                                                               <  4:150563/147‑1 (MQ=255)                                                               gAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGGTAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTc                                                                                       >  2:287212/1‑149 (MQ=255)                                                                 ggAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGa                                                                                     >  3:272473/1‑149 (MQ=255)                                                                           ttCGAGGCCGCTGGCAAGATCGTCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCa                                                                           <  6:59605/149‑1 (MQ=255)                                                                                 gCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAAc                                                                     >  7:10539/1‑149 (MQ=255)                                                                                     cTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTcc                                                                 <  8:153286/149‑1 (MQ=255)                                                                                     cTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTcc                                                                 <  3:48174/149‑1 (MQ=255)                                                                                      tGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCg                                                                >  2:223777/1‑149 (MQ=255)                                                                                      tGGCAAGATCCGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCg                                                                >  2:236188/1‑149 (MQ=255)                                                                                           aGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAAt                                                           <  6:93727/149‑1 (MQ=255)                                                                                             aTCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATAc                                                         >  5:239378/1‑149 (MQ=255)                                                                                                    ccaccaTCTTCCACCCGGTCGGTACCTGCCGCATGGGCAACGGTGCCATCGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTcc                                                  >  6:233106/1‑149 (MQ=255)                                                                                                          tCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCAcc                                            >  1:223211/1‑149 (MQ=255)                                                                                                            ttCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCt                                          >  1:98744/1‑149 (MQ=255)                                                                                                            ttCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCt                                          >  1:7332/1‑149 (MQ=255)                                                                                                              ccACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCAT‑GACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCtgat                                       >  6:30502/1‑149 (MQ=255)                                                                                                              ccACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCAT‑GACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCtgat                                       >  6:35204/1‑149 (MQ=255)                                                                                                                     gtcggtACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCAt                                                                                                                       <  5:285324/63‑1 (MQ=255)                                                                                                                     gtcggtACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCAt                                                                                                                       >  6:285324/1‑63 (MQ=255)                                                                                                                      tcggtACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGcc                                <  4:272473/149‑1 (MQ=255)                                                                                                                         gtACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGa                              >  3:134654/1‑148 (MQ=255)                                                                                                                          tACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGaa                             >  7:214672/1‑148 (MQ=255)                                                                                                                            ccTGTCGCATGGGCAACGGTGCCAT‑GACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAg                          <  5:35204/148‑1 (MQ=255)                                                                                                                            ccTGTCGCATGGGCAACGGTGCCAT‑GACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAg                          <  5:30502/148‑1 (MQ=255)                                                                                                                             cTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAgg                         <  2:7332/149‑1 (MQ=255)                                                                                                                                tCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAggcg                       <  8:10539/148‑1 (MQ=255)                                                                                                                                  gCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAggcggc                     <  6:239378/148‑1 (MQ=255)                                                                                                                                  gCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGgcggcg                    >  8:194145/1‑149 (MQ=255)                                                                                                                                   cATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGgcggcg                    >  6:218406/1‑148 (MQ=255)                                                                                                                                    aTGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGgcgc                   >  6:92365/1‑148 (MQ=255)                                                                                                                                        gCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACt               >  2:37694/1‑148 (MQ=255)                                                                                                                                           aCGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATc           <  8:214672/149‑1 (MQ=255)                                                                                                                                                  cATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCAAAg    <  2:98744/149‑1 (MQ=255)                                                                                                                                                   aTGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCa       >  7:238943/1‑145 (MQ=255)                                                                                                                                                   aTGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCa       <  8:238943/145‑1 (MQ=255)                                                                                                                                                    tGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCAAAGGa  <  5:233106/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCCGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACCATCTTCCACCCGGTCGGTACCTGTCGCATGGGCAACGGTGCCATGGACGTTGTGGATAACCAGCTGCGCGTGCATGGCATCCCCGGTCTACGTGTAGCCGATGCTTCGATCATGCCGCAGATAACTTCCGGCAATACGTGTTCCCCCACCCTGATGATCGCCGAAAAGGCGGCGCAACTGATCCTCAAAGGA  >  NC_002947/66252‑66547

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