Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	539,975	A→G	D258D (GAT→GAC)	allC ←	allantoate deiminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	539,975	0	A	G	100.0%	337.8 / NA	91	D258D (GAT→GAC)	allC	allantoate deiminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (59/32);  total (59/32)

AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGT  >  NZ_CP009273/539890‑540061                                                                                      |                                                                                       aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                    >  1:1653570/1‑90 (MQ=255) aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                    >  2:604433/1‑90 (MQ=255) aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                    >  1:281724/1‑90 (MQ=255)          tGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCttt                                                                           >  2:2912844/1‑90 (MQ=255)          tGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCttt                                                                           >  2:1053829/1‑90 (MQ=255)          tGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCttt                                                                           >  2:358677/1‑90 (MQ=255)            aaCGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTc                                                                         >  1:1998693/1‑90 (MQ=255)              cGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGc                                                                       >  2:1029049/1‑90 (MQ=255)                tGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCtt                                                                     >  2:895/1‑90 (MQ=255)                tGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCtt                                                                     >  1:1171426/1‑90 (MQ=255)                       aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                              >  2:39838/1‑90 (MQ=255)                       aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                              >  1:2001818/1‑90 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCa                                                                                 <  1:2530279/66‑1 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCa                                                                                 >  2:2530279/1‑66 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                         >  1:2725079/1‑90 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                         >  2:1640964/1‑90 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                         >  2:2074477/1‑90 (MQ=255)                            gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                         >  1:865988/1‑90 (MQ=255)                            gCACCAAATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                         >  2:2315299/1‑90 (MQ=255)                             caccacATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTg                                                        <  1:358677/90‑1 (MQ=255)                                 acaTTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATgg                                                    <  1:996617/90‑1 (MQ=255)                                    tttACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCaa                                                 >  1:396663/1‑90 (MQ=255)                                      tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAAt                                               >  1:2369478/1‑90 (MQ=255)                                       aCCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATg                                              >  2:383500/1‑90 (MQ=255)                                         cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                            >  2:2084950/1‑90 (MQ=255)                                           tATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGAc                                          >  2:1889075/1‑90 (MQ=255)                                            aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  2:903682/1‑90 (MQ=255)                                            aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  2:1407921/1‑90 (MQ=255)                                            aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  2:2696918/1‑90 (MQ=255)                                            aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  1:1350977/1‑90 (MQ=255)                                              tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:1640964/90‑1 (MQ=255)                                              tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:2912844/90‑1 (MQ=255)                                              tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:688393/90‑1 (MQ=255)                                              tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:1185206/90‑1 (MQ=255)                                              tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:39838/90‑1 (MQ=255)                                                gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                     >  2:2562875/1‑90 (MQ=255)                                                gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                     <  1:1446529/90‑1 (MQ=255)                                                gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                     <  1:736905/90‑1 (MQ=255)                                                  gcgcGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGc                                   >  1:861933/1‑90 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTg                                                        >  1:458815/1‑67 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTg                                                        <  2:458815/67‑1 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGGCTGAAAGCGt                                 >  1:2288833/1‑90 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGt                                 >  1:1773473/1‑90 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGt                                 >  1:1991873/1‑90 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGt                                 >  1:2676553/1‑90 (MQ=255)                                                    gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGt                                 >  2:2143344/1‑90 (MQ=255)                                                      ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                               >  2:329825/1‑90 (MQ=255)                                                      ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                               <  1:2315299/90‑1 (MQ=255)                                                        ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  1:1108760/90‑1 (MQ=255)                                                        ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  1:1407921/90‑1 (MQ=255)                                                        ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  1:2143344/90‑1 (MQ=255)                                                        ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  2:1350977/90‑1 (MQ=255)                                                               ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  1:1540905/83‑1 (MQ=255)                                                               ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             >  2:1540905/1‑83 (MQ=255)                                                               ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                      >  1:2412816/1‑90 (MQ=255)                                                               ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                      >  1:153424/1‑90 (MQ=255)                                                               ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                      <  1:745152/90‑1 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:689696/1‑90 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  1:647396/1‑90 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:2638502/1‑90 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:337515/1‑90 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    <  2:153424/90‑1 (MQ=255)                                                                 gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATAGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:1563125/1‑90 (MQ=255)                                                                      aGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:1020141/66‑1 (MQ=255)                                                                      aGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       >  2:1020141/1‑66 (MQ=255)                                                                      aGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGAt               <  2:2134449/90‑1 (MQ=255)                                                                       ggTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATa              <  1:319927/90‑1 (MQ=255)                                                                        gTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATaa             <  2:1991873/90‑1 (MQ=255)                                                                          cAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAAcc           >  2:2641451/1‑90 (MQ=255)                                                                           aGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  2:404238/61‑1 (MQ=255)                                                                           aGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       >  1:404238/1‑61 (MQ=255)                                                                           aGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:2390620/1‑80 (MQ=255)                                                                           aGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    <  1:2390620/80‑1 (MQ=255)                                                                             aaCCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCAt        <  1:1250317/90‑1 (MQ=255)                                                                              aCCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATc       <  1:1889075/90‑1 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGTGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  2:183879/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATa              >  1:2162063/1‑80 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATa              <  2:2162063/80‑1 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  1:1743532/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  2:1959270/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  2:18489/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  1:1532203/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  2:1793151/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  1:358029/1‑90 (MQ=255)                                                                                 aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc    >  1:2030696/1‑90 (MQ=255)                                                                                  gCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCg   >  2:572320/1‑90 (MQ=255)                                                                                  gCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCg   >  2:1670935/1‑90 (MQ=255)                                                                                   cGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCgt  <  1:2641451/90‑1 (MQ=255)                                                                                   cGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCgt  <  1:2612975/90‑1 (MQ=255)                                                                                   cGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCgt  <  2:339143/90‑1 (MQ=255)                                                                                   cGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCgt  <  2:1773473/90‑1 (MQ=255)                                                                                      |                                                                                       AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGT  >  NZ_CP009273/539890‑540061

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