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

TTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGTTCGG  >  NC_000913/4077338‑4077599                                                                                                                                    |                                                                                                                                   tttCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                             >  2:865937/1‑139 (MQ=255)         aaCGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                               >  2:706927/1‑129 (MQ=255)         aaCGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                               <  1:706927/129‑1 (MQ=255)           cGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACCagaga                                                                                                                   <  2:52488/139‑1 (MQ=255)           cGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagaga                                                                                                                   <  1:843423/139‑1 (MQ=255)              tGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                                <  2:768448/139‑1 (MQ=255)                 cGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                             <  1:1829/139‑1 (MQ=255)                   tCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                           <  1:453083/139‑1 (MQ=255)                     aGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                         <  1:865937/139‑1 (MQ=255)                                                                      aTCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACt                                                        >  2:628208/1‑139 (MQ=255)                                                                         gCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                         >  2:212135/1‑87 (MQ=255)                                                                         gCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                         <  1:212135/87‑1 (MQ=255)                                                                          cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                               >  2:80846/1‑64 (MQ=255)                                                                          cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                               <  1:80846/64‑1 (MQ=255)                                                                             gTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                          <  2:808154/98‑1 (MQ=255)                                                                             gTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                          >  1:808154/1‑98 (MQ=255)                                                                             gTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGtt                                                 <  1:628208/139‑1 (MQ=255)                                                                              tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGTGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                          <  1:552470/113‑1 (MQ=255)                                                                              tCAATGATGGCCCCGTGACATTCTGGTTGCAGCTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaag                                                                          >  2:552470/1‑113 (MQ=255)                                                                                  tgatgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACaga                                                                            <  2:210563/107‑1 (MQ=255)                                                                                  tgatgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACaga                                                                            >  1:210563/1‑107 (MQ=255)                                                                                     tgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                                <  2:330745/68‑1 (MQ=255)                                                                                     tgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagagag                                                                                                                >  1:330745/1‑68 (MQ=255)                                                                                     tgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                    <  1:754670/80‑1 (MQ=255)                                                                                     tgGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                    >  2:754670/1‑80 (MQ=255)                                                                                       gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                           >  2:625425/1‑71 (MQ=255)                                                                                       gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                           <  1:625425/71‑1 (MQ=255)                                                                                       gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                       >  1:825553/1‑139 (MQ=255)                                                                                        ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                         >  1:340368/1‑72 (MQ=255)                                                                                        ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                         <  2:340368/72‑1 (MQ=255)                                                                                             tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTg                                 >  2:843188/1‑139 (MQ=255)                                                                                                       gTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                               >  2:564280/1‑35 (MQ=255)                                                                                                       gTTGCAGGTATGAGCCAGCTCCCAGGGTTATCACg                                                                                                                               <  1:564280/35‑1 (MQ=255)                                                                                                          gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGt                                                            >  1:674592/1‑99 (MQ=255)                                                                                                          gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGt                                                            <  2:674592/99‑1 (MQ=255)                                                                                                              gTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCa                >  2:90373/1‑139 (MQ=255)                                                                                                                 tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                             >  2:820520/1‑43 (MQ=255)                                                                                                                 tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                             <  1:820520/43‑1 (MQ=255)                                                                                                                   aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGt                                                                                                               >  1:752287/1‑39 (MQ=255)                                                                                                                   aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGt                                                                                                               <  2:752287/39‑1 (MQ=255)                                                                                                                   aGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCaa           >  1:897578/1‑139 (MQ=255)                                                                                                                        gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                          >  1:736808/1‑39 (MQ=38)                                                                                                                        gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                          <  2:736808/39‑1 (MQ=38)                                                                                                                         cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCg                                                             <  2:21005/83‑1 (MQ=255)                                                                                                                         cTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCg                                                             >  1:21005/1‑83 (MQ=255)                                                                                                                            ccAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAGCACCTTCGGTTTCCACAAACTGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGTTCgg  >  2:573335/1‑139 (MQ=255)                                                                                                                             cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt                  <  2:577255/122‑1 (MQ=255)                                                                                                                             cAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCAt                  >  1:577255/1‑122 (MQ=255)                                                                                                                              aGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCg                               >  2:552040/1‑108 (MQ=255)                                                                                                                              aGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCg                               <  1:552040/108‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   TTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGTTCGG  >  NC_000913/4077338‑4077599

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