Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,565	A→G	Y117H (TAC→CAC)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,565	0	A	G	100.0%	146.1 / NA	42	Y117H (TAC→CAC)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/20);  total (22/20)

TAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGT  >  NC_000913/3515420‑3515706                                                                                                                                                  |                                                                                                                                              tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGt                                                                                                                                          >  1:609808/1‑151 (MQ=255)    gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCa                                                                                                                                       <  2:387498/151‑1 (MQ=255)         acaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAatat                                                                                                                                  >  2:154771/1‑151 (MQ=255)            aTAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaa                                                                                                                                >  1:319188/1‑150 (MQ=255)                aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                             >  1:458432/1‑149 (MQ=255)                           cTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGccc                                                                                                                  >  1:327991/1‑149 (MQ=255)                              cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacgg                                                                                                             >  1:141556/1‑151 (MQ=255)                              cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacgg                                                                                                             >  2:667215/1‑151 (MQ=255)                                   tttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacgg                                                                                                        >  2:418377/1‑151 (MQ=255)                                     tCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCGCCCGacggaacgg                                                                                                        >  2:403535/1‑149 (MQ=255)                                     tCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacgga                                                                                                       >  2:495568/1‑150 (MQ=255)                                     tCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAACGCCCGACGGAACGGAt                                                                                                      >  2:26620/1‑151 (MQ=255)                                        gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                                                                                   <  2:22405/151‑1 (MQ=255)                                         cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGCCGGAACGGAt                                                                                                      >  1:224145/1‑147 (MQ=255)                                          gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                                                                                 <  2:21041/151‑1 (MQ=255)                                                    aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAACGCCCGACGGAACGGGTTCTGGCTTttggtg                                                                                        >  2:424093/1‑150 (MQ=255)                                                     aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgttg                                                                                        >  1:349130/1‑149 (MQ=255)                                                          gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActc                                                                                   <  2:596704/149‑1 (MQ=255)                                                              ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgcg                                                                             <  1:456171/151‑1 (MQ=255)                                                              ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGCTCAAATATAAAAACAGTACCGCCCGACGGCACGGGTCCTCGCTTTTGTTGAACTCTTcgc                                                                              >  1:371891/1‑150 (MQ=255)                                                               cGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCgg                                                                            <  2:584692/151‑1 (MQ=255)                                                                       aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATa                                                                    <  2:544554/151‑1 (MQ=255)                                                                        gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg                                                                   <  2:406583/151‑1 (MQ=255)                                                                             ttttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTaaa                                                               <  2:404966/150‑1 (MQ=255)                                                                                   aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCCGCGGAACGGAACCTGGCTTTTGTTGAACTCTTCGCCGAACTGATAATAAACCTCGgc                                                        >  2:37871/1‑151 (MQ=255)                                                                                       gCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCaa                                                    <  1:355652/151‑1 (MQ=255)                                                                                        ccGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAt                                                   <  2:501045/151‑1 (MQ=255)                                                                                          gAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAttt                                                 <  1:541661/151‑1 (MQ=255)                                                                                            aCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGt                                               >  2:244446/1‑151 (MQ=255)                                                                                              ggCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGGACCGCCCGACGGAACGGAGCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTGttt                                             >  2:248890/1‑151 (MQ=255)                                                                                                   gTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTTTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCggg                                         >  1:43841/1‑150 (MQ=255)                                                                                                           cACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaac                                  <  1:263988/149‑1 (MQ=255)                                                                                                           cACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAg                                <  1:274819/151‑1 (MQ=255)                                                                                                               gcgcAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTcg                            <  2:161876/151‑1 (MQ=255)                                                                                                                          aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGcc                   >  1:285559/1‑149 (MQ=255)                                                                                                                          aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCt                  <  1:39913/150‑1 (MQ=255)                                                                                                                                acaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt            >  2:592898/1‑150 (MQ=255)                                                                                                                                  aGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgta           >  1:354845/1‑149 (MQ=255)                                                                                                                                  aGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACt         <  2:354849/151‑1 (MQ=255)                                                                                                                                  aGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACt         <  1:5738/151‑1 (MQ=255)                                                                                                                                         tGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAg   <  2:578946/150‑1 (MQ=255)                                                                                                                                           tacccGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGATCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGt  <  1:397498/146‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              TAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGT  >  NC_000913/3515420‑3515706

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