Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,667	A→G	10.0%	intergenic (+51/‑49)	PP_t32 → / → PP_t33	tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,271,667	0	A	G	10.0%	133.3 / 5.3	50	intergenic (+51/‑49)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (3/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810                                                                                                                                                    |                                                                                                                                                  gCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGGCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACAt                                                                                                                                                  <  6:287976/149‑1 (MQ=255)   cAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGc                                                                                                                                                >  5:75075/1‑149 (MQ=255)        agGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCCTTgt                                                                                                                                           <  5:230138/149‑2 (MQ=255)             cGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGa                                                                                                                                     <  1:233835/149‑1 (MQ=255)               cACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATc                                                                                                                                   >  4:325292/1‑149 (MQ=255)                aCTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCATTAACAGATGGGTGCGTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATc                                                                                                                                   >  1:2869/1‑148 (MQ=255)                  tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                         >  8:47899/1‑140 (MQ=255)                  tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                         <  7:41886/140‑1 (MQ=255)                  tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                         <  7:47899/140‑1 (MQ=255)                  tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                         >  8:41886/1‑140 (MQ=255)                    cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGGCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGt                                                                                                                              <  4:122902/149‑1 (MQ=255)                      cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTGCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGttt                                                                                                                            <  7:94030/149‑1 (MQ=255)                       aGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTc                                                                                                                           <  2:90158/149‑1 (MQ=255)                              gTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                     <  3:325292/148‑1 (MQ=255)                                tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                   <  6:75075/148‑1 (MQ=255)                                       gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCt                                                                                                           >  5:255240/1‑149 (MQ=255)                                                     cTGGCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGc                                                                                             <  6:312999/149‑1 (MQ=255)                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            <  6:255240/147‑1 (MQ=255)                                                          aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                         <  3:334411/148‑1 (MQ=255)                                                            gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGtt                                                                                      <  5:256689/149‑1 (MQ=255)                                                                    tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            <  7:308388/135‑1 (MQ=255)                                                                    tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            >  8:308388/1‑135 (MQ=255)                                                                      gcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCATTAAAAGATCCTTGCAAACCAACCGTCTAAAACCACAAAACACATCCCTTT‑ATCTGATAAAATTTTCAATCCGTCACAAGCATGAACGATACGCAGCGGTAGTTCAATCGGTTa                                                                            >  6:137753/1‑149 (MQ=255)                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   <  3:172186/118‑1 (MQ=255)                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   >  4:172186/1‑118 (MQ=255)                                                                                    cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGt                                                              >  2:332439/1‑149 (MQ=255)                                                                                            gcgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                      >  3:347392/1‑149 (MQ=255)                                                                                             cgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg‑g                                                    <  4:347392/149‑1 (MQ=255)                                                                                             cgcCATTCTTCAGTATCACCTGGGCGCCTGGCACCTGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCATTCGGGCTGAATCCTGATCTATTCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg‑g                                                    <  6:135979/149‑1 (MQ=25)                                                                                                cATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGt                                                 >  5:215784/1‑149 (MQ=255)                                                                                                   ttttCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACAGGCCTGTCACGCCGG‑GGGTcgc                                              >  1:228514/1‑149 (MQ=255)                                                                                                        cgcaaCAGATCGGTGCCTGGCACCGGGCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGTTa                                                                            <  6:348940/112‑1 (MQ=25)                                                                                                             cAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGt                                    >  1:281581/1‑149 (MQ=255)                                                                                                             cAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGt                                    >  1:286440/1‑149 (MQ=255)                                                                                                                aTGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTccc                                 >  2:174392/1‑149 (MQ=255)                                                                                                                  gggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGt                               >  5:100718/1‑149 (MQ=255)                                                                                                                   ggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTc                              <  1:332439/149‑1 (MQ=255)                                                                                                                               cGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCgg                                                                   <  6:255250/101‑1 (MQ=255)                                                                                                                                       aGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          <  2:286440/149‑1 (MQ=255)                                                                                                                                       aGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          <  2:281581/149‑1 (MQ=255)                                                                                                                                        agcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGttt           >  3:151036/5‑148 (MQ=14)                                                                                                                                         gcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGtt                                         <  5:260048/114‑1 (MQ=17)                                                                                                                                         gcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGtt                                         >  6:260048/4‑117 (MQ=17)                                                                                                                                         ccACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGCGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                       <  1:225757/135‑1 (MQ=255)                                                                                                                                         ccACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                       >  2:225757/1‑134 (MQ=255)                                                                                                                                          cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          >  6:144136/3‑147 (MQ=255)                                                                                                                                          cacaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg       >  5:53560/1‑149 (MQ=255)                                                                                                                                           acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  3:46372/1‑149 (MQ=255)                                                                                                                                           acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  4:212397/1‑149 (MQ=255)                                                                                                                                               ggCACATTCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACTATACGCTGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  <  2:2869/149‑1 (MQ=255)                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  >  2:233904/1‑146 (MQ=255)                                                                                                                                                  acaTGCCTTG‑CGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTAAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGCGTACGAGTCCCGTCCGCTGCGCCATCTTCGTATCAAGGTccc  >  6:250803/1‑146 (MQ=17)                                                                                                                                                    |                                                                                                                                                  GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810

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