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%	208.8 / NA	64	*427* (TAG→TAA)	yjiN	zinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/31);  total (33/31)

TGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATT  >  NC_000913/4565838‑4566095                                                                                                                                 |                                                                                                                                  tGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGc                                                                                                                         <  2:1215505/139‑1 (MQ=255)      aGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCTGAgg                                                                                                                    >  1:213999/1‑139 (MQ=255)           cATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGaa                                                                                                               <  1:1419502/139‑1 (MQ=255)               aTTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAg                                                                                                           >  2:163518/1‑139 (MQ=255)                  ttccttccCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGcc                                                                                                        >  2:320339/1‑139 (MQ=255)                           cGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGc                                                                                               >  2:1007955/1‑139 (MQ=255)                                cACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAg                                                                                          <  1:163518/139‑1 (MQ=255)                                 aCGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         <  1:653901/139‑1 (MQ=255)                                  cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         >  2:433259/1‑138 (MQ=255)                                  cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         <  1:433259/138‑1 (MQ=255)                                  cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                        <  1:320339/139‑1 (MQ=255)                                               aaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAACAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                           >  2:396433/1‑139 (MQ=255)                                               aaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                              >  2:1454297/1‑120 (MQ=255)                                               aaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                              <  1:1454297/120‑1 (MQ=255)                                                   agGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCa                                                                       >  2:2946/1‑139 (MQ=255)                                                    ggCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCg                                                                                                       >  1:272116/1‑106 (MQ=255)                                                    ggCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCg                                                                                                       <  2:272116/106‑1 (MQ=255)                                                         caTTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACaa                                                                 <  1:396433/139‑1 (MQ=255)                                                                   aTTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCa                                                       >  1:1388857/1‑139 (MQ=255)                                                                    ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                      >  2:129941/1‑139 (MQ=255)                                                                    ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                      >  2:1106646/1‑139 (MQ=255)                                                                    ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                      >  1:1103445/1‑139 (MQ=255)                                                                      cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTg                                                                                                                          <  2:946674/69‑1 (MQ=255)                                                                      cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTg                                                                                                                          >  1:946674/1‑69 (MQ=255)                                                                      cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                    >  2:288308/1‑139 (MQ=255)                                                                      cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                    <  2:213999/139‑1 (MQ=255)                                                                        tttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                            >  2:256299/1‑129 (MQ=255)                                                                        tttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                            <  1:256299/129‑1 (MQ=255)                                                                               tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCg                                                              <  2:510669/120‑1 (MQ=255)                                                                               tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCg                                                              >  1:510669/1‑120 (MQ=255)                                                                               tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                           <  1:1007955/139‑1 (MQ=255)                                                                                      tttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                    >  2:308449/1‑123 (MQ=255)                                                                                      tttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                    <  1:308449/123‑1 (MQ=255)                                                                                       ttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                   <  2:168720/139‑1 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                        <  2:691548/50‑1 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                        >  1:691548/1‑50 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGa                                                                                                                      >  1:459808/1‑52 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGa                                                                                                                      <  2:459808/52‑1 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCg                                                              >  1:832141/1‑108 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCg                                                              <  2:832141/108‑1 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                           <  2:640932/127‑1 (MQ=255)                                                                                           tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                           >  1:640932/1‑127 (MQ=255)                                                                                                ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAg                                                                                          <  2:664220/75‑1 (MQ=255)                                                                                                ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAg                                                                                          >  1:664220/1‑75 (MQ=255)                                                                                                ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGa                          >  1:56811/1‑139 (MQ=255)                                                                                                    cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                         >  1:697195/1‑104 (MQ=255)                                                                                                    cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                         <  2:697195/104‑1 (MQ=255)                                                                                                    cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCtt                      <  2:56811/139‑1 (MQ=255)                                                                                                         aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                           >  2:321954/1‑113 (MQ=255)                                                                                                         aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                           <  1:321954/113‑1 (MQ=255)                                                                                                         aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                 >  2:348334/1‑139 (MQ=255)                                                                                                              cTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAtt            <  2:1103445/139‑1 (MQ=255)                                                                                                              cTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAtt            <  1:289910/139‑1 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         >  2:14867/1‑59 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         <  1:14867/59‑1 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATaaa                                 <  1:494100/115‑1 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATaaa                                 >  2:494100/1‑115 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                 <  2:734077/131‑1 (MQ=255)                                                                                                                 aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                 >  1:734077/1‑131 (MQ=255)                                                                                                                       tCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTa          <  1:183723/132‑1 (MQ=255)                                                                                                                       tCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTa          >  2:183723/1‑132 (MQ=255)                                                                                                                        cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         <  1:886783/52‑1 (MQ=255)                                                                                                                        cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                         >  2:886783/1‑52 (MQ=255)                                                                                                                        cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTAtt  >  2:1370557/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                  TGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATT  >  NC_000913/4565838‑4566095

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