Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,565,966	C→T	*427* (TAG→TAA)	yjiN ←	zinc‑type alcohol dehydrogenase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,565,966	0	C	T	100.0%	248.5 / NA	75	*427* (TAG→TAA)	yjiN	zinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/37);  total (38/37)

GACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTG  >  NC_000913/4565839‑4566097                                                                                                                                |                                                                                                                                    gACAAGTGACATGATTTTCCTTCCCCGTACCACGGGCGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGAGTTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                          <  1:1089791/139‑1 (MQ=255) gACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                          <  1:755550/139‑1 (MQ=255)  aCAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                         >  2:560602/1‑139 (MQ=255)      gTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAggg                                                                                                                     >  1:1190263/1‑139 (MQ=255)                           gTACCACGGCAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                                <  2:1190263/139‑1 (MQ=255)                            tACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGa                                                                                               >  1:711820/1‑139 (MQ=255)                                 cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                          <  2:1654227/107‑1 (MQ=255)                                 cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                          >  1:1654227/1‑107 (MQ=255)                                 cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTg                                                                                                  >  2:1643610/1‑131 (MQ=255)                                 cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTg                                                                                                  <  1:1643610/131‑1 (MQ=255)                                 cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                          >  1:165096/1‑139 (MQ=255)                                  gggAGATAATGAAAAAAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCaa                                                                                         >  1:1334869/1‑139 (MQ=255)                                         aaTGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                      <  1:947150/103‑1 (MQ=255)                                         aaTGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                      >  2:947150/1‑103 (MQ=255)                                           tGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                <  2:711820/139‑1 (MQ=255)                                             aaaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTa                                                                              >  2:1260877/1‑139 (MQ=255)                                              aaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                             >  2:871134/1‑139 (MQ=255)                                                    gCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAAt                                                                       >  2:209086/1‑139 (MQ=255)                                                            attaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                               >  2:25225/1‑91 (MQ=255)                                                            attaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                               <  1:25225/91‑1 (MQ=255)                                                            attaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                           >  1:938892/1‑111 (MQ=255)                                                            attaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                           <  2:938892/111‑1 (MQ=255)                                                             ttaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                              <  2:165096/139‑1 (MQ=255)                                                                     cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                >  1:1732706/1‑81 (MQ=255)                                                                     cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                <  2:1732706/81‑1 (MQ=255)                                                                     cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                      <  1:560602/139‑1 (MQ=255)                                                                        ttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGc                                                       >  1:169112/1‑135 (MQ=255)                                                                        ttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGc                                                       <  2:169112/135‑1 (MQ=255)                                                                              tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGc                                                                                                           <  2:358119/77‑1 (MQ=255)                                                                              tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGc                                                                                                           >  1:358119/1‑77 (MQ=255)                                                                                      acagtctctctcGCTGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGATCAAGGCCGGGAGCTGCGACAGCAATTATAAAATTAGCCCAATACAACCGCCTACCAGCGTACCGTTGAGACGAATa                                     <  1:1059775/135‑1 (MQ=255)                                                                                          tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                             <  1:1272469/95‑1 (MQ=255)                                                                                          tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                             >  2:1272469/1‑95 (MQ=255)                                                                                          tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATaaa                                   >  1:1192043/1‑137 (MQ=255)                                                                                          tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATaaa                                   <  2:1192043/137‑1 (MQ=255)                                                                                          tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                                 >  2:300623/1‑139 (MQ=255)                                                                                                tcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                >  1:1283518/1‑54 (MQ=255)                                                                                                tcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                <  2:1283518/54‑1 (MQ=255)                                                                                                tcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                           <  1:300540/123‑1 (MQ=255)                                                                                                tcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                           >  2:300540/1‑123 (MQ=255)                                                                                                   cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                      <  2:921266/45‑1 (MQ=255)                                                                                                   cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                      >  1:921266/1‑45 (MQ=255)                                                                                                     gTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAACATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCCGATcctgt                      >  2:1059775/1‑135 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                               >  1:35423/1‑47 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                               <  2:35423/47‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                           >  2:419909/1‑67 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                           <  1:419909/67‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                             >  1:154748/1‑81 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                             <  2:154748/81‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                              <  1:829644/96‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                              >  2:829644/1‑96 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                      <  2:1048862/104‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                      >  1:1048862/1‑104 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                             <  1:652879/113‑1 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                             >  2:652879/1‑113 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCAAGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                >  1:1540894/1‑78 (MQ=255)                                                                                                        aGCCACTTAGTTTTTCAAGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                <  2:1540894/78‑1 (MQ=255)                                                                                                          ccACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGat                 >  2:1405874/1‑139 (MQ=255)                                                                                                             cTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAtt              <  1:1405874/139‑1 (MQ=255)                                                                                                             cTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAtt              <  1:300623/139‑1 (MQ=255)                                                                                                                aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                     <  1:1574502/81‑1 (MQ=255)                                                                                                                aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                     >  2:1574502/1‑81 (MQ=255)                                                                                                                     ttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                     <  2:492105/76‑1 (MQ=255)                                                                                                                     ttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                     >  1:492105/1‑76 (MQ=255)                                                                                                                       cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGaa                                       <  2:846666/104‑1 (MQ=255)                                                                                                                       cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGaa                                       >  1:846666/1‑104 (MQ=255)                                                                                                                         tGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTg  <  1:871134/139‑1 (MQ=255)                                                                                                                          ggATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGa                                        >  2:647601/1‑100 (MQ=255)                                                                                                                          ggATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGa                                        <  1:647601/100‑1 (MQ=255)                                                                                                                           gATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                   >  2:303954/1‑120 (MQ=255)                                                                                                                           gATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                   <  1:303954/120‑1 (MQ=255)                                                                                                                            aTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGtt                                               >  2:1453900/1‑91 (MQ=255)                                                                                                                            aTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGtt                                               <  1:1453900/91‑1 (MQ=255)                                                                                                                            aTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCa                              >  2:148768/1‑108 (MQ=255)                                                                                                                            aTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCa                              <  1:148768/108‑1 (MQ=255)                                                                                                                                |                                                                                                                                    GACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTG  >  NC_000913/4565839‑4566097

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