Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,710	A→G	7.7%	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,710	0	A	G	7.7%	151.7 / 2.9	52	intergenic (+51/‑49)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base G (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAAC  >  NC_002947/2271565‑2271859                                                                                                                                                  |                                                                                                                                                       cAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATAGGTGCCTGGCACCGGTCTGCAGCCACAACGCACATGcc                                                                                                                                                    >  1:655853/1‑150 (MQ=255)   gagagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCtt                                                                                                                                                  >  2:8310/1‑150 (MQ=255)      agGAACGCACTGCGCAGCTGTAGTTCAGTCGGTTAGAATACCACCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCAGTTTTCAGTAACGGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                              <  1:488156/150‑1 (MQ=25)                       cGGTAGTTCAGTCGGTTAGAATACAGCGGCGTGACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAAt                                                                                                                             >  1:1124686/1‑150 (MQ=37)                        ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                            <  1:8310/150‑1 (MQ=255)                        ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                            <  2:714116/150‑1 (MQ=255)                         gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCCTTgtggctt                                                                                                                                           >  2:874256/1‑128 (MQ=21)                         gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCCTTgtggctt                                                                                                                                           <  1:874256/134‑7 (MQ=21)                          tAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                          >  2:727669/1‑150 (MQ=255)                          tAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                          >  1:992441/1‑150 (MQ=255)                             ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgcg                                                                                                                       >  2:1110728/1‑150 (MQ=255)                                  tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGc                                                                                                                  <  2:992441/150‑1 (MQ=255)                                                ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCggg                                                                                                                         >  1:536468/1‑129 (MQ=255)                                                ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCggg                                                                                                                         <  2:536468/129‑1 (MQ=255)                                                   cTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                                    >  2:122132/1‑131 (MQ=255)                                                   cTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                                    <  1:122132/131‑1 (MQ=255)                                                   cTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                                 <  1:1110728/150‑1 (MQ=255)                                                      tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                              >  1:617565/1‑150 (MQ=255)                                                                 gTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                                   <  2:617565/150‑1 (MQ=255)                                                                   cgcgGGTGCGAGTCCCGTCCGCCGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACAAGTCTGAAGCCACAAGGCACATGACTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGGCTGAACGATACGCAGCGGTAGTTCAGTCGGTTa                                                                                 <  2:553028/150‑1 (MQ=38)                                                                         ttCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAc                                                                           <  2:1066320/150‑1 (MQ=255)                                                                         ttCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAc                                                                           >  2:792335/1‑150 (MQ=255)                                                                          tCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAcc                                                                          <  1:727669/150‑1 (MQ=255)                                                                                 ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGt                                                                   >  1:843428/1‑150 (MQ=255)                                                                                  cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAATGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCATGCCTGAACGATACGCAGCGGTAGTTCAGTCCGTTAGAATACCGGCCTGTc                                                                  >  2:603408/1‑150 (MQ=255)                                                                                  cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                             >  1:790371/1‑107 (MQ=255)                                                                                  cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                             <  2:790371/107‑1 (MQ=255)                                                                                       gCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGcc                                                             >  1:1206821/1‑150 (MQ=255)                                                                                         tGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                           <  1:719762/150‑1 (MQ=255)                                                                                         tGCGCCATTTTCCAGGAACTGCTGGGTGCCTGGCACTGGTCCCAAGCCACAAGGCACATGCCTTG‑CGCGGTTCCCAGCTTCAAACGGGCGGAAGCCTACACGCTACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                           <  2:1015144/150‑1 (MQ=21)                                                                                             ccATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAAt                                                                                                                             >  1:1172928/1‑80 (MQ=255)                                                                                             ccATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAAt                                                                                                                             <  2:1172928/80‑1 (MQ=255)                                                                                                      aGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTc                                              >  2:719672/1‑150 (MQ=255)                                                                                                        tAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                            >  2:1025330/1‑150 (MQ=255)                                                                                                        tAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                            >  1:462072/1‑150 (MQ=255)                                                                                                         aaCAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAg                                           <  2:462072/150‑1 (MQ=255)                                                                                                                gggTGCCAGGCACCGGTCTGAAGCCACAAGTCACATGCGTTG‑GGCGGATCCCAGTCTCAATCGAGCGCAAGCCAGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTc                                    <  1:740572/150‑1 (MQ=18)                                                                                                                 ggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTc                                              <  1:438701/139‑1 (MQ=255)                                                                                                                 ggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTc                                              >  2:438701/1‑139 (MQ=255)                                                                                                                     ccTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTTAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTg                               >  2:246808/1‑150 (MQ=255)                                                                                                                       tGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                             >  1:215826/1‑150 (MQ=255)                                                                                                                       tGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                             >  2:1045673/1‑150 (MQ=255)                                                                                                                           aCCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCAt                         <  1:719672/150‑1 (MQ=255)                                                                                                                                   gAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGttt                 <  2:215826/150‑1 (MQ=255)                                                                                                                                   gAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGttt                 >  1:529282/1‑150 (MQ=255)                                                                                                                                   gAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGttt                 >  1:537901/1‑150 (MQ=255)                                                                                                                                      gCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa              <  2:1206821/150‑1 (MQ=255)                                                                                                                                      agcccaAGGCATGTGCCTTG‑TGCTGATTCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa              <  2:503548/146‑1 (MQ=14)                                                                                                                                        cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg             >  1:882479/3‑150 (MQ=14)                                                                                                                                        cacaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACACCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt                   <  1:103361/143‑1 (MQ=255)                                                                                                                                        cacaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCCGCCTGTCACACCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt                   >  2:103361/1‑143 (MQ=255)                                                                                                                                         acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGCTTCAAggc           >  2:592822/1‑149 (MQ=255)                                                                                                                                             ggaacaTGCTTTG‑CGCTGATCCCAGCTTCAATCGGGCGCAAGCCTCACGTCCTCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCt       <  2:1124686/147‑1 (MQ=17)                                                                                                                                                acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTACCTAGa    >  2:365752/1‑150 (MQ=255)                                                                                                                                                  aTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAac  >  1:882831/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAAC  >  NC_002947/2271565‑2271859

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