Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,991,839	G→A	W229* (TGG→TGA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,991,839	0	G	A	98.3%	177.7 / NA	58	W229* (TGG→TGA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (24/33);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACC  >  NC_000913_3_bme_pgi/3991693‑3991984                                                                                                                                                   |                                                                                                                                                   ctgctgCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACt                                                                                                                                                  >  2:414826/1‑149 (MQ=255)   gctgctTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTgg                                                                                                                                                >  2:660310/1‑149 (MQ=255)   gctcctTGACGAATTTTAGCGTACCGCCGTGCGTCGCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTgg                                                                                                                                                <  1:273801/149‑1 (MQ=255)           cagAATTTTATCGTACGGATGTTCTTTTCGCGGGTAAGCGTATTTTGTGGAATACGGTGCCGTGCGACGAAGAAGAGGATAACCACGACTATGTTATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAAGGACTGGATCTCGGt                                                                                                                                        <  2:11733/147‑1 (MQ=255)             gAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCGGACGCCAAATGAATGGCTGGATCTCGGTgg                                                                                                                                      >  1:402744/1‑149 (MQ=255)             gAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGACGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTgg                                                                                                                                      <  1:696398/149‑1 (MQ=255)                gtttATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCtt                                                                                                                                   <  1:522297/148‑1 (MQ=255)                  ttATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTaa                                                                                                                                 <  1:585147/149‑1 (MQ=255)                    tgagTATCGCGGGGCGCCTCGCTTGTAAGTGTCTTCTGTGGATTCTGGTGTCGTGCGACGAAGAAGAGAATTACGCCGGCCAGCTGATGACGCTTTTCGCGCAGGTCGTGCTGACGCCAAATCAAGGACTGGATCTCGGTGGCTTAAGc                                                                                                                               <  2:733663/146‑1 (MQ=255)                    aTCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGCCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGc                                                                                                                               >  1:324354/1‑149 (MQ=255)                       gTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTgg                                                                                                                                      >  2:368936/1‑139 (MQ=255)                       gTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTgg                                                                                                                                      <  1:368936/139‑1 (MQ=255)                                   gattcGCCGGGACGCGGAGACTGTGGAAGATGGTGCCGGGCGACGAAGGAGAGTATTACGACGACGATGTGAGGACGCTTTACGCGCAGGGCGCGCTGACGCCAAATGGCTGACTGTAGCTCGGTGGCTTAAGCTCGGTTTCTGCTGaa                                                                                                                <  2:806895/146‑1 (MQ=255)                                   gTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGaa                                                                                                                <  1:765357/149‑1 (MQ=255)                                       cGCCGGTAAGCGTTTGCTGTGTAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGt                                                                                                            <  2:488117/149‑1 (MQ=255)                                        gCCGGTAAGCGTATTCGGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTGTCTGCTGAAGAGTa                                                                                                           <  1:481279/149‑1 (MQ=255)                                         tcGGTAAGCGGATTCTGTCGAATATGGGGCCGTGCGACGTAGAAGAGTATTACGACGACTATGTGATGACGCTTTACGCTCAGGGCGTGCTGACGCCAAATGAATTACTGGATCTCGGTGGCTCAAGCTCGCTTTCTGCTGAAGAGTAc                                                                                                          <  1:636600/148‑1 (MQ=255)                                            gTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACttt                                                                                                       <  1:687007/149‑1 (MQ=255)                                                 cGTAGTCTGCGGAATATGGTGCCGTGGGCCGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGTTCTCGGTGGCTTTAGCTCGCTCTCTGCTGAAGAGTACTTTGGTGc                                                                                                  <  1:662982/149‑1 (MQ=255)                                                   tATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCa                                                                                                <  1:829693/149‑1 (MQ=255)                                                   tATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTGAAGTACCCTTTCTGCTGAAGATCGCTTTGGTcccg                                                                                                >  2:809741/1‑148 (MQ=255)                                                        tgtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCCGCCtt                                                                                           >  1:259186/1‑149 (MQ=255)                                                                gtGGTGCCGTGCGTCGAAGAAGAGCATTCCGCCGACTTTGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGAAAGct                                                                                   <  2:432073/148‑1 (MQ=255)                                                                 tGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTACGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTgtcggcgc                                                                                  >  2:513102/1‑142 (MQ=255)                                                                  ggTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGccct                                                                                 >  2:485590/1‑149 (MQ=255)                                                                   gTGCCGTGCGACGAAGAAGACCATTACGCCGACTATGTGATGACGCTTTACGCGCAGGGCGTGCGGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTa                                                                                <  1:74445/149‑1 (MQ=255)                                                                    tGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTAc                                                                               >  1:47504/1‑149 (MQ=255)                                                                        gTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAACCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAag‑a                                                                          <  1:553499/149‑1 (MQ=255)                                                                            gacgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTAt                                                                      >  1:220234/1‑149 (MQ=255)                                                                             acgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATc                                                                     >  2:546830/1‑149 (MQ=255)                                                                               gaagaaGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATCAATGACTGGATCTCGGTGGCTTTAGCGCGCTTTCTGCGGAAGAGTATTTTGGCGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGa                                                                   >  2:421293/1‑149 (MQ=255)                                                                                   aagaGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCt                                                                                                                          >  2:17160/1‑91 (MQ=255)                                                                                   aagaGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCt                                                                                                                          <  1:17160/91‑1 (MQ=255)                                                                                         aTTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCCGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATAc                                                         >  2:221466/1‑149 (MQ=255)                                                                                          ttACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGaa                                                                                                                <  2:183628/94‑1 (MQ=255)                                                                                          ttACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGaa                                                                                                                >  1:183628/1‑94 (MQ=255)                                                                                          ttACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACa                                                        >  1:178607/1‑149 (MQ=255)                                                                                             cgacgaCTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAg                                                     <  2:477338/149‑1 (MQ=255)                                                                                             cgacgaCTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAg                                                     <  1:546830/149‑1 (MQ=255)                                                                                              gacgacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAGG‑AGAATCGGTTCCCAATACaaacc                                                    >  1:189647/1‑147 (MQ=255)                                                                                              gacgacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGc                                                    >  1:724466/1‑149 (MQ=255)                                                                                               acgacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCg                                                   <  1:775462/149‑1 (MQ=255)                                                                                                cgacTATGGGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTGTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCgg                                                  <  1:129610/149‑1 (MQ=255)                                                                                                 gacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTAc                                                                               <  1:830684/120‑1 (MQ=255)                                                                                                 gacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTAc                                                                               >  2:830684/1‑120 (MQ=255)                                                                                                  acTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTa                                                <  1:458534/149‑1 (MQ=255)                                                                                                                      gcgcAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGa                            >  2:618598/1‑149 (MQ=255)                                                                                                                         cAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGc                         <  2:220234/149‑1 (MQ=255)                                                                                                                             gCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                     >  2:329793/1‑149 (MQ=255)                                                                                                                                  ccgACGCCAAATGAATGACGGGGTTTTGGTGGCTTAAGCTCGCTTTCTGCTGAAGGGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                <  2:324354/147‑1 (MQ=255)                                                                                                                                  cTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACCAGCATCATTG‑TGCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                >  1:404225/1‑149 (MQ=255)                                                                                                                                   cgaagCCACATGAATGACTGGTTTTCGGTGGCGTAAGCGCGCTTTCTGCGGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTAAAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCGGGAAGCCTATTACTgg               <  2:371461/145‑1 (MQ=255)                                                                                                                                     aCGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTcc                                                               >  1:569512/1‑99 (MQ=255)                                                                                                                                     aCGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTcc                                                               <  2:569512/99‑1 (MQ=255)                                                                                                                                       gCCAAATCAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAAt           <  1:364919/149‑1 (MQ=255)                                                                                                                                        ccAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa          <  1:340235/149‑1 (MQ=255)                                                                                                                                           aatgaatgACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc       <  1:471492/149‑1 (MQ=255)                                                                                                                                                atgACTGGATCTCGGTGGCCTAAGCTCCCTTTCAGATGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTATACAAG‑AGTATCGATTCCCCATACAAAGAGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAcc  >  2:524598/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAG‑AGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACC  >  NC_000913_3_bme_pgi/3991693‑3991984

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