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

TATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077336‑4077595                                                                                                                                      |                                                                                                                               tATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACg                                                                                                                           <  1:596885/139‑1 (MQ=255)     tcgtcgAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGaa                                                                                                                       >  2:97701/1‑139 (MQ=255)           aaCGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  1:97701/139‑1 (MQ=255)               cTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAagagaga                                                                                                             >  2:1086365/1‑139 (MQ=255)                  ccGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTa                                                                                                          <  1:110913/139‑1 (MQ=255)                       aGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCt                                                                                                     >  2:185689/1‑139 (MQ=255)                               aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       >  1:1062578/1‑129 (MQ=255)                               aTGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCg                                                                                                       <  2:1062578/129‑1 (MQ=255)                                       gCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCggg                                                                                     >  1:928049/1‑139 (MQ=255)                                             aGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               >  1:857107/1‑123 (MQ=255)                                             aGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               <  2:857107/123‑1 (MQ=255)                                               gACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  2:1011222/1‑114 (MQ=255)                                               gACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:1011222/114‑1 (MQ=255)                                                 cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                <  2:115259/118‑1 (MQ=255)                                                 cGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTa                                                                                                >  1:115259/1‑118 (MQ=255)                                                       gCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaaga                                                                     <  1:1086365/139‑1 (MQ=255)                                                         tGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 <  1:863088/93‑1 (MQ=255)                                                         tGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAaga                                                                                                                 >  2:863088/1‑93 (MQ=255)                                                                atGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            <  1:185689/139‑1 (MQ=255)                                                                     ggTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaagaa                                                                    >  2:895725/1‑126 (MQ=255)                                                                     ggTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACagaagaa                                                                    <  1:895725/126‑1 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTa                                                                                                    >  2:459800/1‑87 (MQ=255)                                                                            cTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTa                                                                                                    <  1:459800/87‑1 (MQ=255)                                                                               gTCACTGATGGACCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAATAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGCATTAGAGCGTTACTATCAGtt                                             >  2:659910/1‑139 (MQ=255)                                                                                tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         >  2:520370/1‑78 (MQ=255)                                                                                tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         <  1:520370/78‑1 (MQ=255)                                                                                tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  2:494058/81‑1 (MQ=255)                                                                                tCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  1:494058/1‑81 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACaaa                                                                            <  1:680608/98‑1 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACaaa                                                                            >  2:680608/1‑98 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAAt                                                                <  1:887879/110‑1 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAAt                                                                >  2:887879/1‑110 (MQ=255)                                                                                         gCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaa                                   >  1:1113829/1‑139 (MQ=255)                                                                                          ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      >  2:211456/1‑71 (MQ=255)                                                                                          ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGc                                                                                                      <  1:211456/71‑1 (MQ=255)                                                                                          ccccGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGaaa                                  >  1:911347/1‑139 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            >  1:925949/1‑109 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTaga                                                            <  2:925949/109‑1 (MQ=255)                                                                                              gTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGtt                              >  2:109441/1‑139 (MQ=255)                                                                                               tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACt                                                    >  1:644216/1‑116 (MQ=255)                                                                                               tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACt                                                    <  2:644216/116‑1 (MQ=255)                                                                                               tGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTg                             >  1:349401/1‑139 (MQ=255)                                                                                                gACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGc                            >  2:797896/1‑139 (MQ=255)                                                                                                   aTTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTa                         <  1:659910/139‑1 (MQ=255)                                                                                                        ggTTGCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGccc                    <  1:797896/139‑1 (MQ=255)                                                                                                           tGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         <  1:688030/35‑1 (MQ=37)                                                                                                           tGCAGGTATGAGCCAGCTTCCAGGGTTATCACggg                                                                                                                         >  2:688030/1‑35 (MQ=37)                                                                                                            gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                      >  1:1014289/1‑69 (MQ=255)                                                                                                            gCAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCgg                                                                                      <  2:1014289/69‑1 (MQ=255)                                                                                                             cAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               <  2:86858/59‑1 (MQ=255)                                                                                                             cAGGTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTAt                                                                                               >  1:86858/1‑59 (MQ=255)                                                                                                               ggTATGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATc             >  1:430045/1‑139 (MQ=255)                                                                                                                   tGAGCCAGCTTCCAGGGTTATCACGGGAAGCAAGAGAGAGTAt                                                                                                         >  1:418918/1‑43 (MQ=38)                                                                                                                   tGAGCCAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTAt                                                                                                         <  2:418918/43‑1 (MQ=255)                                                                                                                       ccAGCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCaaaa     >  2:409749/1‑139 (MQ=255)                                                                                                                          gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCg                                                         >  2:488779/1‑84 (MQ=255)                                                                                                                          gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCg                                                         <  1:488779/84‑1 (MQ=255)                                                                                                                          gCTTCCAGGGTTATCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGt  >  2:877514/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                               TATTTCGTCGAACGCTGCCGTCAGCAAGAGATGAACACGCAAACAGGACGCTTCGCTGCGGATATGCAGGTATCGCTGGTCAATGATGGCCCCGTGACATTCTGGTTGCAGGTATGAGCCAGCTTCCAGGGTTGTCACGGGAAACAAGAGAGAGTATCGCTATGTATCACCTTCGGGTTCCACAAACAGAAGAAGAATTAGAGCGTTACTATCAGTTTCGCTGGGAAATGTTGCGTAAGCCCCTGCATCAACCAAAAGGT  >  NC_000913/4077336‑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: ‑