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

TGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCC  >  NC_000913/4565838‑4566099                                                                                                                                  |                                                                                                                                      tGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGc                                                                                                                             >  1:1112335/1‑139 (MQ=255)                           cGTACCACGGGAGATAATGAAAAGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGc                                                                                                   >  2:980951/1‑139 (MQ=255)                                  cGGGAGATAATGAAAAGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                            <  1:117550/139‑1 (MQ=255)                                            ataaaaGAGGCCAACATAATT‑CCATTTTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                  <  1:122146/137‑1 (MQ=255)                                               aaaGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                               <  2:41625/139‑1 (MQ=255)                                               aaaGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                               <  1:980951/139‑1 (MQ=255)                                                   agGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                                 >  1:696743/1‑101 (MQ=255)                                                   agGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                                 <  2:696743/101‑1 (MQ=255)                                                    ggCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        <  1:791823/93‑1 (MQ=255)                                                    ggCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        >  2:791823/1‑93 (MQ=255)                                                     gCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAAt                                                                         >  1:1121276/1‑139 (MQ=255)                                                     gCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAAt                                                                         >  2:697470/1‑139 (MQ=255)                                                                    aTTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        >  2:1199360/1‑78 (MQ=255)                                                                    aTTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        <  1:1199360/78‑1 (MQ=255)                                                                     ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTTATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                          >  2:902295/1‑139 (MQ=255)                                                                     ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        <  1:971226/77‑1 (MQ=255)                                                                     ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                        >  2:971226/1‑77 (MQ=255)                                                                       cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCaa                                                              <  2:955211/133‑1 (MQ=255)                                                                       cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCaa                                                              >  1:955211/1‑133 (MQ=255)                                                                       cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                        <  2:1121276/139‑1 (MQ=255)                                                                          ttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTAc                                                     >  1:696222/1‑139 (MQ=255)                                                                                  gTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                  >  2:1082140/1‑102 (MQ=255)                                                                                  gTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                             >  2:463550/1‑139 (MQ=255)                                                                                  gTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                             >  2:1320807/1‑139 (MQ=255)                                                                                  gTTTATTTGATCTCTCTCACGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                  <  1:1082140/102‑1 (MQ=255)                                                                                       tttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAAt                                        >  2:1000128/1‑139 (MQ=255)                                                                                        ttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                       >  2:439249/1‑139 (MQ=255)                                                                                           aTCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAAc                                    >  2:289860/1‑139 (MQ=255)                                                                                            tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                                   >  2:433491/1‑139 (MQ=255)                                                                                            tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                                   >  2:1508875/1‑139 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTaaaaa                                                                                                                                 <  1:392687/44‑1 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTaaaaa                                                                                                                                 >  2:392687/1‑44 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCg                                                                                                           >  2:845615/1‑66 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCg                                                                                                           <  1:845615/66‑1 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                                <  1:260412/109‑1 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                                >  2:260412/1‑109 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                               >  1:1288495/1‑126 (MQ=255)                                                                                             ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                               <  2:1288495/126‑1 (MQ=255)                                                                                                 ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                             <  2:880982/76‑1 (MQ=255)                                                                                                 ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                             >  1:880982/1‑76 (MQ=255)                                                                                                 ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                       <  2:1198763/130‑1 (MQ=255)                                                                                                 ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                       >  1:1198763/1‑130 (MQ=255)                                                                                                     cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                             >  2:394835/1‑104 (MQ=255)                                                                                                     cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                             <  1:394835/104‑1 (MQ=255)                                                                                                        ttAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGc                       <  2:1326914/139‑1 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAg                                                                                                               >  2:556444/1‑49 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAg                                                                                                               <  1:556444/49‑1 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGGCAGCAAATATAAAATTAGACCAATACAACGGCGAACAAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                     >  2:160065/1‑139 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                  >  2:1294195/1‑78 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                  <  1:1294195/78‑1 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                       <  2:1159269/89‑1 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                       >  1:1159269/1‑89 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                     <  1:289860/139‑1 (MQ=255)                                                                                                          aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                     >  1:545270/1‑139 (MQ=255)                                                                                                               cTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAtt                <  1:1320807/139‑1 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                                  >  2:268456/1‑56 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                                  <  1:268456/56‑1 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                       <  1:265091/67‑1 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                       >  2:265091/1‑67 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATc                            >  1:485476/1‑126 (MQ=255)                                                                                                                ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATc                            <  2:485476/126‑1 (MQ=255)                                                                                                                   gTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAAc            <  1:1508875/139‑1 (MQ=255)                                                                                                                      tttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                             >  1:356163/1‑71 (MQ=255)                                                                                                                      tttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                             <  2:356163/71‑1 (MQ=255)                                                                                                                      tttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTa              >  1:846944/1‑134 (MQ=255)                                                                                                                      tttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTa              <  2:846944/134‑1 (MQ=255)                                                                                                                          aTGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTAttt     >  1:978256/1‑139 (MQ=255)                                                                                                                             gATTTTTAAAAATTGCCGAGGGGGAAGAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGcc  <  2:978256/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      TGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCC‑ACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCC  >  NC_000913/4565838‑4566099

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