Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,733,703	G→A	*1539* (TAG→TAA)	lhr →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,733,703	0	G	A	100.0%	53.3 / NA	18	*1539* (TAG→TAA)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/13);  total (5/13)

GCACTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACGCTGAT  >  NC_000913/1733569‑1733830                                                                                                                                       |                                                                                                                                gCACTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGaa                                                                                                                             <  2:460944/139‑1 (MQ=255)                    ggCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGt                                                                                                          <  1:690381/139‑1 (MQ=255)                     gCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGtt                                                                                                         <  2:673608/139‑1 (MQ=255)                            tttACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATcc                                                                                                  <  2:701729/139‑1 (MQ=255)                                      cAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTg                                                                                        <  1:670053/139‑1 (MQ=255)                                               aTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGa                                                                               >  2:541239/1‑139 (MQ=255)                                                   tCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGattat                                                                           >  1:41992/1‑139 (MQ=255)                                                            ccGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTcc                                                                  >  2:612654/1‑139 (MQ=255)                                                             cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCg                                                                 >  1:527686/1‑139 (MQ=255)                                                                       cGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGa                                                       <  2:527686/139‑1 (MQ=255)                                                                         aTGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACt                                                     <  1:535125/139‑1 (MQ=255)                                                                                cGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTAttt                                              <  2:177305/139‑1 (MQ=255)                                                                                  ctgctgCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAg                                            <  2:41992/139‑1 (MQ=255)                                                                                       gcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCa                                       <  1:612654/139‑1 (MQ=255)                                                                                        cGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAg                                      <  1:541239/139‑1 (MQ=255)                                                                                                                 cacaAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGc                                        <  2:416103/112‑1 (MQ=255)                                                                                                                 cacaAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGc                                        >  1:416103/1‑112 (MQ=255)                                                                                                                            gATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACGCTGAt  <  2:479489/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                GCACTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACGCTGAT  >  NC_000913/1733569‑1733830

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