Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,733,200	A→G	T1372A (ACG→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,733,200	0	A	G	100.0%	107.9 / NA	32	T1372A (ACG→GCG)	lhr	putative ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/20);  total (12/20)

ATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCG  >  NC_000913/1733067‑1733336                                                                                                                                      |                                                                                                                                         aTTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATAT‑CCTGGCGGGTTTCCATCGATGCAAGCGCTTTg                                                                                                                                   <  2:355703/140‑1 (MQ=255)   tGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGt                                                                                                                                  <  1:652203/140‑1 (MQ=255)     cAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCg                                                                                                                                <  1:498514/140‑1 (MQ=255)       ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCgga                                                                                                                              >  1:622017/1‑138 (MQ=255)       ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGaa                                                                                                                              <  1:615938/140‑1 (MQ=255)       ggTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGaa                                                                                                                              <  1:439256/140‑1 (MQ=255)          ggAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTa                                                                                                                           >  1:362477/1‑140 (MQ=255)            aGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTCCGGAATAATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCTTGCCCGGTGTCCATCGAGGCAAGTGCTTTGTCGAAgaagg                                                                                                                         >  2:301140/1‑136 (MQ=255)                    aaCGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTc                                                                                                                 <  2:596831/140‑1 (MQ=255)                                   aTGCTGGCGCTGGCGGAAAATATGTTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCg                                                                                                  >  1:16690/1‑140 (MQ=255)                                      ctggcgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAgg                                                                                               <  2:367577/140‑1 (MQ=255)                                          cgctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGt                                                                                           <  2:422102/140‑1 (MQ=255)                                            ctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGCTAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCGAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGTATTATGCGAGGTCGttt                                                                                         >  1:408878/1‑140 (MQ=255)                                                               aCCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGgc                                                                      >  2:411212/1‑140 (MQ=255)                                                                 cGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTAGTTTTGTAGAAGGTCTGGGTGgcgc                                                                    <  1:346227/140‑1 (MQ=255)                                                                            gacaggtcgtcAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTg                                                         <  1:488103/135‑1 (MQ=255)                                                                             atcaGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGa                                                        >  1:528832/1‑140 (MQ=255)                                                                                aGTCGTCAGGCGGTGATAGCCGAAATTATTCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACg                                                     <  1:321011/140‑1 (MQ=255)                                                                                  tcgtcAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTc                                                   >  1:413226/1‑140 (MQ=255)                                                                                  tcgtcAGGCGGTGATAGCCGAAAATATCCCTGGCGGGCTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTc                                                   <  1:267340/140‑1 (MQ=255)                                                                                   cgtcAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTCTAGAAGGTCTGGGTGGCGCGCGATTCGCTGAACGTCt                                                  >  1:249202/1‑140 (MQ=255)                                                                                      cAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGAc                                               <  1:378349/140‑1 (MQ=255)                                                                                      cAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGAc                                               <  2:362477/140‑1 (MQ=255)                                                                                            gTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGCCGTTTTGTCGAAGGGCTGGGTGGCGCGCAATTCGCTGAACGCCTGACtctgtc                                         >  2:488103/1‑135 (MQ=255)                                                                                                      aaaaTATCCCTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTgcgc                               >  2:21582/1‑140 (MQ=255)                                                                                                             cccTGGCGGGTTTCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTgg                        <  2:622017/140‑1 (MQ=255)                                                                                                                       tttCCATCGATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGc              >  1:509795/1‑140 (MQ=255)                                                                                                                               gATGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGAAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCaaa      <  1:612850/140‑1 (MQ=255)                                                                                                                                 tGCAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAAcg    >  1:54529/1‑140 (MQ=255)                                                                                                                                 tGCAAGCGCTTTGTCCAAGTATGGAGGACCCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAAcg    <  1:248857/140‑1 (MQ=255)                                                                                                                                   cAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAAcgcg  <  1:510729/140‑1 (MQ=255)                                                                                                                                   cAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAAcgcg  <  2:528832/140‑1 (MQ=255)                                                                                                                                   cAAGCGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAAcgcg  <  1:365428/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ATTGCAGGTGGAGCCACTAAACGATACCGAAAGGATGCTGGCGCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCG  >  NC_000913/1733067‑1733336

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