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%	223.7 / NA	68	L7L (TTG→TTA)	yiiD	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/33);  total (35/33)

ACTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077334‑4077595                                                                                                                                        |                                                                                                                               aCTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCa                                                                                                                             >  2:847903/1‑139 (MQ=255)   tATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  2:1638833/139‑1 (MQ=255)     tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         >  2:293905/1‑139 (MQ=255)       tcgtcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaa                                                                                                                       <  1:969953/139‑1 (MQ=255)             aaCGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  1:80408/139‑1 (MQ=255)                    ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  2:664647/122‑1 (MQ=255)                    ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           >  1:664647/1‑122 (MQ=255)                    ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTa                                                                                                          >  1:1393818/1‑139 (MQ=255)                    ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTa                                                                                                          <  1:847903/139‑1 (MQ=255)                    ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTa                                                                                                          >  1:757693/1‑139 (MQ=255)                     cGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         <  2:822993/139‑1 (MQ=255)                        cAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  2:1225227/139‑1 (MQ=255)                              gagaTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                >  1:166036/1‑139 (MQ=255)                               agaTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               <  1:293905/139‑1 (MQ=255)                                gaTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  2:148295/130‑1 (MQ=255)                                gaTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  1:148295/1‑130 (MQ=255)                                   gAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCAcc                                                                                           >  1:1458341/1‑139 (MQ=255)                                    aaCACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCA‑CTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  2:1517067/1‑126 (MQ=255)                                    aaCACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCA‑CTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:1517067/126‑1 (MQ=255)                                          cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:1768970/100‑1 (MQ=255)                                          cAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           >  2:1768970/1‑100 (MQ=255)                                                ggACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCcaca                                                                              <  2:1393818/139‑1 (MQ=255)                                                   cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  2:446180/1‑111 (MQ=255)                                                   cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  1:446180/111‑1 (MQ=255)                                                     cTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  2:1447501/1‑110 (MQ=255)                                                     cTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:1447501/110‑1 (MQ=255)                                                      ttCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACaga                                                                        <  2:1458341/139‑1 (MQ=255)                                                      ttCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACaga                                                                        >  1:977192/1‑139 (MQ=255)                                                        cGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  2:166036/139‑1 (MQ=255)                                                          cTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGGGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaagaa                                                                    <  1:1086836/139‑1 (MQ=255)                                                               gATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  1:130229/1‑100 (MQ=255)                                                               gATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  2:130229/100‑1 (MQ=255)                                                                      aGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGt                                                        <  2:977192/139‑1 (MQ=255)                                                                         tATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  1:1507887/122‑1 (MQ=255)                                                                         tATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      >  2:1507887/1‑122 (MQ=255)                                                                          aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCa                                                                                                                             <  1:1346108/66‑1 (MQ=255)                                                                          aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCa                                                                                                                             >  2:1346108/1‑66 (MQ=255)                                                                          aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                >  2:1730324/1‑95 (MQ=255)                                                                          aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                <  1:1730324/95‑1 (MQ=255)                                                                          aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACt                                                    >  2:228465/1‑139 (MQ=255)                                                                             gCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATc                                                 >  2:139821/1‑139 (MQ=255)                                                                                  tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGttt                                            >  1:240370/1‑139 (MQ=255)                                                                                           gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 >  1:316656/1‑61 (MQ=255)                                                                                           gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  2:316656/61‑1 (MQ=255)                                                                                           gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                   >  2:481170/1‑139 (MQ=255)                                                                                            ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaaa                                  >  1:850667/1‑139 (MQ=255)                                                                                                gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  1:493965/56‑1 (MQ=255)                                                                                                gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 >  2:493965/1‑56 (MQ=255)                                                                                                gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGttt                                            <  1:449141/125‑1 (MQ=255)                                                                                                gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGttt                                            >  2:449141/1‑125 (MQ=255)                                                                                                 tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            <  1:934513/108‑1 (MQ=255)                                                                                                 tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            >  2:934513/1‑108 (MQ=255)                                                                                                   aCATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCg                           <  2:240370/139‑1 (MQ=255)                                                                                                     aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  1:481170/139‑1 (MQ=255)                                                                                                              gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         >  2:1737643/1‑130 (MQ=255)                                                                                                              gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  1:1737643/130‑1 (MQ=255)                                                                                                               cAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCa               >  1:292812/1‑139 (MQ=255)                                                                                                                aGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt              >  1:1046888/1‑139 (MQ=255)                                                                                                                 ggTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATc             >  2:643779/1‑139 (MQ=255)                                                                                                                 ggTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATc             >  2:659391/1‑139 (MQ=255)                                                                                                                  gTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCa            >  1:675685/1‑139 (MQ=255)                                                                                                                     tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                <  2:1291855/52‑1 (MQ=255)                                                                                                                     tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                >  1:1291855/1‑52 (MQ=255)                                                                                                                     tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAAcc         <  2:1046888/139‑1 (MQ=255)                                                                                                                       aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCaa       <  1:139821/139‑1 (MQ=255)                                                                                                                            gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGt  >  1:92010/1‑139 (MQ=255)                                                                                                                             cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  1:95623/70‑1 (MQ=255)                                                                                                                             cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      >  2:95623/1‑70 (MQ=255)                                                                                                                                        |                                                                                                                               ACTATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077334‑4077595

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