Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,731,408	A→G	A774A (GCA→GCG)	lhr →	putative ATP‑dependent helicase Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,731,408	0	A	G	100.0%	178.8 / NA	50	A774A (GCA→GCG)	lhr	putative ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/22);  total (28/22)

TGATGACGGCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATC  >  NC_000913/1731288‑1731541                                                                                                                         |                                                                                                                                       tgatgaCGGCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTc                                                                                                                     >  1:1055782/1‑140 (MQ=255)   atgaCGGCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTcgc                                                                                                                   >  1:1342664/1‑140 (MQ=255)        ggCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATg                                                                                                              >  1:1053724/1‑140 (MQ=255)        ggCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATg                                                                                                              >  1:743696/1‑140 (MQ=255)           aTTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATgcgc                                                                                                           >  1:150508/1‑140 (MQ=255)           aTTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATgcgc                                                                                                           >  1:1600335/1‑140 (MQ=255)               ttgCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcg                                                                                                       <  2:1485731/140‑1 (MQ=255)               ttgCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcg                                                                                                       <  1:1415149/140‑1 (MQ=255)                 gCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                                                                     <  2:1055782/140‑1 (MQ=255)                    cGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                                                                     >  2:774784/1‑137 (MQ=255)                    cGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                                                                     <  1:774784/137‑1 (MQ=255)                    cGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                                                                  >  2:833996/1‑140 (MQ=255)                    cGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                                                                  >  2:740902/1‑140 (MQ=255)                    cGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                                                                  <  1:150627/140‑1 (MQ=255)                       aTTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatt                                                                                               >  2:1246258/1‑140 (MQ=255)                       aTTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatt                                                                                               >  1:1338587/1‑140 (MQ=255)                         tCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGcgccgc                                                                                                        <  1:168355/129‑1 (MQ=255)                         tCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGcgccgc                                                                                                        >  2:168355/1‑129 (MQ=255)                              acacCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATg                                                                                       >  1:934185/1‑140 (MQ=255)                               cacCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGc                                                                                      >  1:1285481/1‑140 (MQ=255)                                acCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGcc                                                                                     <  2:1342664/140‑1 (MQ=255)                                acCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGcc                                                                                     >  1:425379/1‑140 (MQ=255)                                 ccGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCg                                                                                    >  2:1601543/1‑140 (MQ=255)                                   gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCa                                                                                                 <  1:502318/126‑1 (MQ=255)                                   gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCa                                                                                                 >  2:502318/1‑126 (MQ=255)                                   gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCggg                                                                                  <  2:656175/140‑1 (MQ=255)                                    aTGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTTAATGCCggg                                                                                  >  1:78236/1‑140 (MQ=255)                                     tGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                                                                  <  2:689422/123‑1 (MQ=255)                                     tGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                                                                  >  1:689422/1‑123 (MQ=255)                                      ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGgcg                                                                               <  1:1246258/140‑1 (MQ=255)                                      ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGgcg                                                                               <  2:1053724/140‑1 (MQ=255)                                      ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGgcg                                                                               <  2:743696/140‑1 (MQ=255)                                      ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGgcg                                                                               <  2:150508/140‑1 (MQ=255)                                      ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGgcg                                                                               <  2:1540450/140‑1 (MQ=255)                                             tGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCg                                                                        <  2:1600335/140‑1 (MQ=255)                                               cccGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCGTTA‑TTAATGCCGGGGCGCACTCCggg                                                                      >  1:1121742/1‑140 (MQ=255)                                               cccGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCggg                                                                      >  1:1557988/1‑140 (MQ=255)                                               cccGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCggg                                                                      >  1:1070660/1‑140 (MQ=255)                                                     gCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCg                                                                <  2:1285481/140‑1 (MQ=255)                                                       cgcgATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCa                                                              >  1:1165377/1‑140 (MQ=255)                                                                 ttgtttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCGCTCCGGGCCATCGCACCCCGCTTTg                                                    <  2:1070660/140‑1 (MQ=255)                                                                              aaaaGTTGCTGCAAATTGTCCGCGAGGCGGTATGCAGGTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGAGCCATCGCACCCCGCTTTGGCAACAACGGCTg                                       >  2:1416592/1‑140 (MQ=255)                                                                               aaaGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTgc                                      <  1:1601543/140‑1 (MQ=255)                                                                                 aGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTgcgc                                    <  1:833996/140‑1 (MQ=255)                                                                                          aaaTTGTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCcagtcag                           <  2:425379/140‑1 (MQ=255)                                                                                               gTCCGCGAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCt                      <  2:1165377/140‑1 (MQ=255)                                                                                                  cgcgAGGCGGTAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGa                   >  1:266242/1‑140 (MQ=255)                                                                                                            tAGGCAGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAg         >  2:191909/1‑140 (MQ=255)                                                                                                                 aGCTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGata    >  2:672863/1‑140 (MQ=255)                                                                                                                   cTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATc  >  2:239800/1‑140 (MQ=255)                                                                                                                   cTCGGCGCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATc  >  1:469190/1‑140 (MQ=255)                                                                                                                         |                                                                                                                                       TGATGACGGCATTGTTGCACGTATTCCTGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTA‑TTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATC  >  NC_000913/1731288‑1731541

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