Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,497	A→T	A137A (GCA→GCT)	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,729,497	0	A	T	100.0%	100.6 / NA	31	A137A (GCA→GCT)	lhr	putative ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/14);  total (17/14)

CCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGAT  >  NC_000913/1729424‑1729568                                                                          |                                                                        ccGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGaa                                                                   <  1:298285/80‑1 (MQ=255) ccGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAAc                                                                  <  2:298286/81‑1 (MQ=255)        aCGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCaa                                                           <  1:298287/81‑1 (MQ=255)          ggcggcggcGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAAc                                                         >  2:298326/1‑81 (MQ=255)          ggcggcggcGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAAc                                                         >  1:298340/1‑81 (MQ=255)             ggcggcgCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCa                                                      >  1:298321/1‑81 (MQ=255)                ggcgcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCAccc                                                   <  2:298288/81‑1 (MQ=255)                  cgcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCg                                                  >  2:298372/1‑80 (MQ=255)                   gcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGt                                                 >  2:298336/1‑80 (MQ=255)                      gcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATc                                             <  1:298289/81‑1 (MQ=255)                        gAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATccg                                           >  1:298328/1‑81 (MQ=255)                        gAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATccg                                           >  1:298342/1‑81 (MQ=255)                          aaCGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATccgcc                                         <  2:298290/81‑1 (MQ=255)                              gAAGTCAATCTTCGCGTAGGGTTCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGat                                     <  1:298291/81‑1 (MQ=255)                                aGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGatat                                   >  2:298341/1‑81 (MQ=255)                                  tCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTc                                 <  1:298292/81‑1 (MQ=255)                                   cAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCt                                >  2:298333/1‑81 (MQ=255)                                     aTCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGa                              >  1:298329/1‑81 (MQ=255)                                          cgcgTAGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATcac                          >  1:298345/1‑80 (MQ=255)                                               aGGGATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACAcc                    >  2:298327/1‑81 (MQ=255)                                                gggATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACAccc                   >  1:298349/1‑81 (MQ=255)                                                gggATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACAccc                   <  2:298293/81‑1 (MQ=255)                                                 ggATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCg                  <  2:298294/81‑1 (MQ=255)                                                  gATCCGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCg                  >  2:298353/1‑80 (MQ=255)                                                     ccGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAAt               >  2:298344/1‑80 (MQ=255)                                                      cGTACTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATc              <  1:298295/80‑1 (MQ=255)                                                          cTGGCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCActc          >  2:298296/1‑80 (MQ=255)                                                            ggCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTAt       >  2:298334/1‑81 (MQ=255)                                                             gCGATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTa        <  2:298297/79‑1 (MQ=255)                                                               gATACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTg    <  1:298298/81‑1 (MQ=255)                                                                 tACGCCTGCTCAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGAt  <  1:298299/81‑1 (MQ=255)                                                                          |                                                                        CCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGAT  >  NC_000913/1729424‑1729568

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