Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,730,559	G→A	P491P (CCG→CCA)	lhr →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,730,559	0	G	A	100.0%	139.0 / NA	37	P491P (CCG→CCA)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/21);  total (16/21)

TGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGGC  >  NC_000913/1730412‑1730697                                                                                                                                                    |                                                                                                                                           tGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATgg                                                                                                                                         >  1:159138/1‑151 (MQ=255)  gTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATgg                                                                                                                                         >  2:3728/1‑150 (MQ=255)   tAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGa                                                                                                                                        >  2:139392/1‑150 (MQ=255)    agagTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGaa                                                                                                                                       <  1:420013/150‑1 (MQ=255)           aTGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTg                                                                                                                                <  1:297450/150‑1 (MQ=255)             gTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGcc                                                                                                                              >  2:561193/1‑150 (MQ=255)                      caggcTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGtt                                                                                                                    >  1:8469/1‑151 (MQ=255)                       aggcTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGttt                                                                                                                   >  1:653161/1‑151 (MQ=255)                            ggAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGAc                                                                                                              >  2:592119/1‑151 (MQ=255)                                  ccTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACg                                                                                                        >  2:69844/1‑151 (MQ=255)                                  ccTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACg                                                                                                        <  1:179862/151‑1 (MQ=255)                                      acacCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTgg                                                                                                    <  2:349191/151‑1 (MQ=255)                                      acacCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTgg                                                                                                    <  2:15561/151‑1 (MQ=255)                                       caccacCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTgg                                                                                                    >  1:282009/1‑150 (MQ=255)                                       caccacCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTgg                                                                                                    >  2:433191/1‑150 (MQ=255)                                       caccacCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGa                                                                                                   <  2:8468/151‑1 (MQ=255)                                        accaccGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGat                                                                                                  >  1:143804/1‑151 (MQ=255)                                         ccaccGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGat                                                                                                  >  1:415424/1‑150 (MQ=255)                                            ccGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGc                                                                                              <  1:155126/151‑1 (MQ=255)                                                 tAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTcc                                                                                         >  2:208603/1‑151 (MQ=255)                                                          cGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTAt                                                                                <  2:74239/151‑1 (MQ=255)                                                           gACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATc                                                                               <  2:134633/151‑1 (MQ=255)                                                             cGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATccc                                                                             <  2:159141/151‑1 (MQ=255)                                                                   tGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGc                                                                       >  1:24550/1‑151 (MQ=255)                                                                             aCCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTg                                                             <  1:427513/151‑1 (MQ=255)                                                                               cGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCt                                                           <  1:401384/151‑1 (MQ=255)                                                                                              gatggatgCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTg                                            >  1:614035/1‑151 (MQ=255)                                                                                                       aTTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAAc                                   <  2:620248/151‑1 (MQ=255)                                                                                                          aCAGGTAGACGAATGGTAGTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCACGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACa                                  <  1:288143/149‑1 (MQ=255)                                                                                                          aCAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAg                                 <  2:338084/150‑1 (MQ=255)                                                                                                                      aTGGTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGata                    <  1:380308/151‑1 (MQ=255)                                                                                                                        ggTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATAtt                  <  2:506220/151‑1 (MQ=255)                                                                                                                         gTACTCCCGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTg                 <  2:39959/151‑1 (MQ=255)                                                                                                                               ccGCGTACGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGgaga                            >  2:410001/1‑134 (MQ=255)                                                                                                                               ccGCGTACGCCGTGCCGCACCATGGAAAGACCTGCCAAGACGTGTTTTTGACGCCTCGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGgaga                            <  1:409996/134‑1 (MQ=255)                                                                                                                                     aCGCCGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCt     <  2:199621/151‑1 (MQ=255)                                                                                                                                        ccGTGCCGCACCATGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGgc  <  1:480375/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                           TGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGGC  >  NC_000913/1730412‑1730697

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