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

TATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077336‑4077595                                                                                                                                      |                                                                                                                               tATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCTGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:596884/139‑1 (MQ=255) tATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:870853/139‑1 (MQ=255)   tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         >  2:305546/1‑139 (MQ=255)   tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:402694/137‑1 (MQ=255)   tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           >  2:402694/1‑137 (MQ=255)    ttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGa                                                                                                                        >  2:592645/1‑139 (MQ=255)     tcgtcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaa                                                                                                                       >  2:993403/1‑139 (MQ=255)      cgtcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaaa                                                                                                                      <  1:1041495/139‑1 (MQ=255)        tcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACa                                                                                                                    <  1:699447/139‑1 (MQ=255)        tcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACa                                                                                                                    >  2:1237982/1‑139 (MQ=255)             cGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagaga                                                                                                               <  1:1300499/139‑1 (MQ=255)                tGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                            <  1:993403/139‑1 (MQ=255)                      cAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:592645/139‑1 (MQ=255)                       aGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     <  1:907313/139‑1 (MQ=255)                               aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 >  1:1252458/1‑119 (MQ=255)                               aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  2:1252458/119‑1 (MQ=255)                                 gAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCAcc                                                                                           >  2:155524/1‑139 (MQ=255)                                  aaCACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     >  1:334630/1‑128 (MQ=255)                                  aaCACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     <  2:334630/128‑1 (MQ=255)                                       gCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     <  2:743511/123‑1 (MQ=255)                                       gCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     >  1:743511/1‑123 (MQ=255)                                               gACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                            >  1:86223/1‑108 (MQ=255)                                               gACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                            <  2:86223/108‑1 (MQ=255)                                                    ttCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACaga                                                                        <  1:1237982/139‑1 (MQ=255)                                                       gCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaaga                                                                     <  1:1028230/139‑1 (MQ=255)                                                                        aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  1:85281/1‑89 (MQ=255)                                                                        aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  2:85281/89‑1 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 >  2:882762/1‑74 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  1:882762/74‑1 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               <  2:681056/92‑1 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               >  1:681056/1‑92 (MQ=255)                                                                                tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGttt                                            >  1:225216/1‑139 (MQ=255)                                                                                  aaTGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGAATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCg                                          >  1:696529/1‑139 (MQ=255)                                                                                    tgatgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCt                                        <  1:815651/139‑1 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACaaa                                                                            >  2:1003453/1‑98 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACaaa                                                                            <  1:1003453/98‑1 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                   >  1:710869/1‑139 (MQ=255)                                                                                          ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaaa                                  >  2:1161249/1‑139 (MQ=255)                                                                                            ccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATg                                <  1:1161249/139‑1 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      <  2:126224/99‑1 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                      >  1:126224/1‑99 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGtt                              >  1:1006934/1‑139 (MQ=255)                                                                                               tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAACCAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTg                             >  1:875758/1‑139 (MQ=255)                                                                                                   aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  2:698105/139‑1 (MQ=255)                                                                                                   aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  2:710869/139‑1 (MQ=255)                                                                                                        ggTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGccc                    <  1:305546/139‑1 (MQ=255)                                                                                                         gTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGcccc                   >  2:374100/1‑139 (MQ=255)                                                                                                              aGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt              <  1:155524/139‑1 (MQ=255)                                                                                                               ggTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACTAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATc             >  1:778723/1‑139 (MQ=255)                                                                                                                          gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGt  >  1:623262/1‑139 (MQ=255)                                                                                                                               cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACa                                                                          >  2:467841/1‑62 (MQ=255)                                                                                                                               cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACa                                                                          <  1:467841/62‑1 (MQ=255)                                                                                                                                      |                                                                                                                               TATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077336‑4077595

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