Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,200,671	G→A	9.7%	intergenic (+19/‑2)	PP_2807 → / → PP_2808	hypothetical protein/oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,200,671	0	G	A	9.7%	107.1 / 2.5	41	intergenic (+19/‑2)	PP_2807/PP_2808	hypothetical protein/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

TGCAACGCCTGCTGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTGACGCGCAGGCCGATGTGATCGTGATCGGCGCCGGCTTT  >  NC_002947/3200525‑3200819                                                                                                                                                   |                                                                                                                                                     tGCAACGCCTGCTGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCgaa                                                                                                                                                    <  1:301940/149‑1 (MQ=255) tGCAACGCCTGCTGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCgaa                                                                                                                                                    <  2:61572/149‑1 (MQ=255)   cAACGCCTGCGGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCgaat                                                                                                                                                   <  8:13489/148‑1 (MQ=255)      cGCCTGCTGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCgaatgaat                                                                                                                                               <  8:272360/149‑1 (MQ=255)                           aGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGAt                                                                                                                          >  5:178962/1‑149 (MQ=255)                           aGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGAt                                                                                                                          >  3:149603/1‑149 (MQ=255)                            gTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATg                                                                                                                         <  1:281492/149‑1 (MQ=255)                            gTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATg                                                                                                                         <  5:162812/149‑1 (MQ=255)                               cGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGt                                                                                                                       <  4:340308/148‑1 (MQ=255)                                gACGCTGCTGGGGTTGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGa                                                                                                                     >  2:88819/1‑149 (MQ=255)                                         ggggTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTg                                                                                                                      >  5:295360/1‑139 (MQ=255)                                         ggggTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTg                                                                                                                      <  6:295360/139‑1 (MQ=255)                                         ggggTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCAcc                                                                                                            >  3:239142/1‑149 (MQ=255)                                               ggtggAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTGTGATAATTCTGGGAGAGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGaa                                                                                                      <  6:73365/149‑1 (MQ=255)                                                         gATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGACATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCgg                                                                                             <  8:47814/148‑1 (MQ=255)                                                           tGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTcgcg                                                                                              <  4:239142/145‑1 (MQ=255)                                                              ccAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTggg                                                                                       >  3:14012/1‑149 (MQ=255)                                                              ccAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTggg                                                                                       >  4:304537/1‑149 (MQ=255)                                                               cAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCAAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTgggg                                                                                      >  5:186894/1‑149 (MQ=255)                                                                     cgcgACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCAAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATg                                                                                 <  6:186894/148‑1 (MQ=255)                                                                      gcgACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATggg                                                                               <  3:304537/149‑1 (MQ=255)                                                                        gACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTa                                                                             <  4:149603/149‑1 (MQ=255)                                                                        gACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTa                                                                             <  6:178962/149‑1 (MQ=255)                                                                           ccGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTagag                                                                          >  8:173934/1‑149 (MQ=255)                                                                                       cTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGa                                                               >  5:5444/1‑148 (MQ=255)                                                                                       cTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGAt                                                              <  1:88819/149‑1 (MQ=255)                                                                                               gCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCAAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGAc                                                                      >  3:124739/1‑133 (MQ=255)                                                                                               gCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCAAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGAc                                                                      <  4:124739/133‑1 (MQ=255)                                                                                                ccGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCg                                                     <  4:14012/149‑1 (MQ=255)                                                                                                  ggCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCa                                                         <  7:183070/143‑1 (MQ=255)                                                                                                  ggCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCa                                                         >  8:183070/1‑143 (MQ=255)                                                                                                    cATTCTCCGAGCCGGGTAAAGGCTTCTGATAATGCTGGGAGTGCCGCGAATGAATGCCATCAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGGCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGccc                                                 <  4:341296/149‑1 (MQ=255)                                                                                                         tcCGAGCCGGGTAAAGGCTTCTGATAAATCTGTGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAACACCGTTAAAGTTTCGCGGTCTGGGGTGATTGGTAGAGCGACCGGCTGATCTGCAGACGGCACtag                                              >  7:52256/1‑145 (MQ=255)                                                                                                             aGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTg                                        <  7:173934/149‑1 (MQ=255)                                                                                                                 gggTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATc                                                             >  7:120738/1‑124 (MQ=255)                                                                                                                 gggTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATc                                                             <  8:120738/124‑1 (MQ=255)                                                                                                                       aGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTGACGCGCAGGc                              >  8:104118/1‑149 (MQ=255)                                                                                                                        ggCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCTATGGTGACTCAGCAGCGTTGAAGTTTCGCTGTCTGGGGGGATGGGTAGAGCGACCGGCTTAt                                                              <  3:196265/116‑1 (MQ=255)                                                                                                                        ggCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGAt                                                              >  4:196265/1‑116 (MQ=255)                                                                                                                           ttatGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCAGGCCCTCGAAGGTGAAGCGCAGGCCGAt                          >  6:352137/4‑149 (MQ=255)                                                                                                                            tCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTGACGCGCAGGCCGAtg                         >  5:258962/1‑149 (MQ=255)                                                                                                                               gATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAg                                           <  3:127997/128‑1 (MQ=255)                                                                                                                               gATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAg                                           >  4:127997/1‑128 (MQ=255)                                                                                                                               gATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTGACGCGCAGGCCGATgtga                      >  2:167814/1‑149 (MQ=255)                                                                                                                                             tGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGa                                                               >  8:280647/1‑94 (MQ=255)                                                                                                                                                   gaatgaatGCCTTGAACAATGCAACCGATGATGACTCAACACCATTAAAATTTCGCGATCTGGGGGGATGGGTAGAACGACCGGATGATATGCAGCCGGCACTAGAAGGTGAAGCGCAGGCCGATGTGATCGTGATCGGCGCCGGGttt  >  8:52549/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGCAACGCCTGCTGGGCTGGAGTGCCAGTTCGACGCTGCTGGGGTGGGTGGAAATCGATGGCCAGATGCGCGACCCGCAGGAACTCCTGGGCATGCCGGCATTCTCCGAGCCGGGTAAAGGCTTCTGATAATTCTGGGAGTGACGCGAATGAATGCCTTGAACAATGCGGCCGATGGTGACTCAGCACCGTTGAAGTTTCGCGGTCTGGGGGGATGGGTAGAGCGACCGGCTGATCTGCAGCCGGCCCTCGAAGGTGACGCGCAGGCCGATGTGATCGTGATCGGCGCCGGCTTT  >  NC_002947/3200525‑3200819

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