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%	92.0 / NA	29	*1539* (TAG→TAA)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/15);  total (14/15)

ACTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTT  >  NC_000913/1733571‑1733818                                                                                                                                     |                                                                                                                    aCTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGAGCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATaaagaaag                                                                                                               <  1:1114824/139‑1 (MQ=255)   tGCGACGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCAGGATTGGGGATAAAGAAAGGa                                                                                                             >  1:937362/1‑139 (MQ=255)            aaCCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGAt                                                                                                    <  1:797023/139‑1 (MQ=255)                       cgcTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                         <  2:937362/139‑1 (MQ=255)                                   aCAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                         >  2:818424/1‑127 (MQ=255)                                   aCAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                         <  1:818424/127‑1 (MQ=255)                                             aTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGa                                                                                                             >  2:1088894/1‑97 (MQ=255)                                             aTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGa                                                                                                             <  1:1088894/97‑1 (MQ=255)                                                 tCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGAc                                                                                                            <  2:1021755/94‑1 (MQ=255)                                                 tCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGAc                                                                                                            >  1:1021755/1‑94 (MQ=255)                                                   tACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                         <  2:1340504/111‑1 (MQ=255)                                                   tACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                         >  1:1340504/1‑111 (MQ=255)                                                    aCCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTcc                                                            >  2:216441/1‑139 (MQ=255)                                                           cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTA‑TGCAGGACGGATTATTCCGCGTCCg                                                     <  2:153441/138‑1 (MQ=255)                                                           cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTA‑TGCAGGACGGATTATTCCGCGTCCg                                                     >  1:153441/1‑138 (MQ=255)                                                                 aCGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCtt                                               <  1:216441/139‑1 (MQ=255)                                                                  cGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTc                                              <  1:16565/139‑1 (MQ=255)                                                                        tGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                 >  1:1574389/1‑98 (MQ=255)                                                                        tGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                 <  2:1574389/98‑1 (MQ=255)                                                                                   ctgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCgg                             >  1:177835/1‑139 (MQ=255)                                                                                   ctgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCgg                             >  1:1111832/1‑139 (MQ=255)                                                                                    tgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTcc                                                            <  1:222596/107‑1 (MQ=255)                                                                                    tgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTcc                                                            >  2:222596/1‑107 (MQ=255)                                                                                           ggCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGAc                                          >  1:1306041/1‑118 (MQ=255)                                                                                           ggCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGAc                                          <  2:1306041/118‑1 (MQ=255)                                                                                                gATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGAt                <  2:1111832/139‑1 (MQ=255)                                                                                                      cAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAgctg          <  2:177835/139‑1 (MQ=255)                                                                                                              ccACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACAGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGtt  >  2:1077828/1‑139 (MQ=255)                                                                                                                               gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGctct                                                >  2:1043918/1‑76 (MQ=255)                                                                                                                               gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGctct                                                <  1:1043918/76‑1 (MQ=255)                                                                                                                                     |                                                                                                                    ACTGCGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTT  >  NC_000913/1733571‑1733818

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