Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,077,469	G→A	L7L (TTG→TTA)	yiiD →	GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,077,469	0	G	A	100.0%	227.0 / NA	70	L7L (TTG→TTA)	yiiD	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/35);  total (35/35)

CTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077335‑4077595                                                                                                                                       |                                                                                                                               cTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCa                                                                                                                             >  2:561843/1‑138 (MQ=255) cTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCa                                                                                                                             <  1:561843/138‑1 (MQ=255) cTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCAc                                                                                                                            >  1:1477164/1‑139 (MQ=255)  tATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:782856/139‑1 (MQ=255)    tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         <  1:1358706/139‑1 (MQ=255)            aaCGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  2:838183/139‑1 (MQ=255)              cGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagaga                                                                                                               <  2:214244/139‑1 (MQ=255)              cGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAACCAGCTTCCAGGGTTATCACGGGAAACAagaga                                                                                                               >  1:737086/1‑139 (MQ=255)                        aGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAACCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     <  2:737086/139‑1 (MQ=255)                                aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         <  1:1227489/127‑1 (MQ=255)                                aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         >  2:1227489/1‑127 (MQ=255)                                       cGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                      <  2:386393/139‑1 (MQ=255)                                        gCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCggg                                                                                     >  1:371340/1‑139 (MQ=255)                                         cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         <  1:366840/102‑1 (MQ=255)                                         cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         >  2:366840/1‑102 (MQ=255)                                         cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaa                                                                                                                       >  2:1388623/1‑104 (MQ=255)                                         cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaa                                                                                                                       <  1:1388623/104‑1 (MQ=255)                                         cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGt                                                                                    >  1:520076/1‑139 (MQ=255)                                                  cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  2:1198720/1‑112 (MQ=255)                                                  cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:1198720/112‑1 (MQ=255)                                                       cGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  2:556341/139‑1 (MQ=255)                                                       cGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  2:1477164/139‑1 (MQ=255)                                                        gCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaaga                                                                     >  1:1206649/1‑139 (MQ=255)                                                          tGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaagaag                                                                   >  1:261666/1‑139 (MQ=255)                                                          tGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaagaag                                                                   >  1:1456975/1‑139 (MQ=255)                                                                 atGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            >  1:216129/1‑139 (MQ=255)                                                                    cAGGTATGGCTGATCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATAACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCg                                                         <  2:371340/139‑1 (MQ=255)                                                                           cGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTAt                                                  >  1:577006/1‑139 (MQ=255)                                                                             cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           >  2:1443414/1‑64 (MQ=255)                                                                             cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:1443414/64‑1 (MQ=255)                                                                                 tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTa                                                              <  2:250140/121‑1 (MQ=255)                                                                                 tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTa                                                              >  1:250140/1‑121 (MQ=255)                                                                                 tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGttt                                            >  1:1461418/1‑139 (MQ=255)                                                                                        tgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTggg                                     >  1:849608/1‑139 (MQ=255)                                                                                          gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  2:1446606/1‑71 (MQ=255)                                                                                          gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  1:1446606/71‑1 (MQ=255)                                                                                          gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                   >  2:914749/1‑139 (MQ=255)                                                                                          gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                   >  1:1168313/1‑139 (MQ=255)                                                                                           ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                      <  1:181631/87‑1 (MQ=255)                                                                                           ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                      >  2:181631/1‑87 (MQ=255)                                                                                            cccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAAt                                 >  2:1401908/1‑139 (MQ=255)                                                                                               gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  2:1288763/1‑66 (MQ=255)                                                                                               gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  1:1288763/66‑1 (MQ=255)                                                                                                    aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  2:1206649/139‑1 (MQ=255)                                                                                                    aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  1:86445/139‑1 (MQ=255)                                                                                                      tCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGc                                                          >  2:1431501/1‑104 (MQ=255)                                                                                                      tCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGc                                                          <  1:1431501/104‑1 (MQ=255)                                                                                                        tGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGcc                     <  2:849608/139‑1 (MQ=255)                                                                                                         ggTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGccc                    <  2:577006/139‑1 (MQ=255)                                                                                                         ggTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGccc                    <  2:520076/139‑1 (MQ=255)                                                                                                          gTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  2:824205/88‑1 (MQ=255)                                                                                                          gTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      >  1:824205/1‑88 (MQ=255)                                                                                                             gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGt                                                        <  2:131835/99‑1 (MQ=255)                                                                                                             gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGt                                                        >  1:131835/1‑99 (MQ=255)                                                                                                               aGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt              <  2:216129/139‑1 (MQ=255)                                                                                                               aGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt              <  2:630439/139‑1 (MQ=255)                                                                                                                      aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         <  1:754485/41‑1 (MQ=255)                                                                                                                      aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         >  2:754485/1‑41 (MQ=255)                                                                                                                      aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTc                                           <  2:274306/103‑1 (MQ=255)                                                                                                                      aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTc                                           >  1:274306/1‑103 (MQ=255)                                                                                                                      aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCaa       >  2:380679/1‑139 (MQ=255)                                                                                                                           gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  2:639484/1‑38 (MQ=38)                                                                                                                           gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  1:639484/38‑1 (MQ=255)                                                                                                                           gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            <  1:399004/81‑1 (MQ=255)                                                                                                                           gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            >  2:399004/1‑81 (MQ=255)                                                                                                                           gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGt  >  1:1240486/1‑139 (MQ=255)                                                                                                                            cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      >  2:441527/1‑70 (MQ=255)                                                                                                                            cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  1:441527/70‑1 (MQ=255)                                                                                                                                cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                >  2:840015/1‑40 (MQ=255)                                                                                                                                cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                <  1:840015/40‑1 (MQ=255)                                                                                                                                       |                                                                                                                               CTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077335‑4077595

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: ‑