Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,692	G→A	A202A (GCG→GCA)	lhr →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,729,692	0	G	A	100.0%	104.0 / NA	33	A202A (GCG→GCA)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

CTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCAACCCGCCCGCAATGCGCCATCCGCAGA  >  NC_000913/1729558‑1729822                                                                                                                                       |                                                                                                                                   ctctATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCactgc                                                                                                                                >  2:807933/1‑139 (MQ=255)   ctATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGTCGTTAAGTCTGGAGCGGCTCGATGCACTGCTc                                                                                                                              <  2:850353/139‑1 (MQ=255)     aTCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCca                                                                                                                            <  2:307913/139‑1 (MQ=255)     aTCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCca                                                                                                                            <  2:600776/139‑1 (MQ=255)          aTGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCt                                                                                                                       <  1:508210/139‑1 (MQ=255)            gCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCa                                                                                                                     >  1:1143548/1‑139 (MQ=255)            gCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCa                                                                                                                     >  1:1243138/1‑139 (MQ=255)                aCCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGcac                                                                                                                 >  1:60081/1‑139 (MQ=255)                           gcgcgAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGc                                                                                                      <  1:1207915/139‑1 (MQ=255)                             gcgAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCt                                                                                                    <  2:60081/139‑1 (MQ=255)                                 aaCGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCt                                                                                                <  2:778789/139‑1 (MQ=255)                                   cGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGc                                                                                              <  1:807933/139‑1 (MQ=255)                                   cGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGc                                                                                              <  2:1135218/139‑1 (MQ=255)                                        cgcgGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACtg                                                                                         >  2:837619/1‑139 (MQ=255)                                                                 atgaAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTgcgc                                                                                    <  1:1426955/119‑1 (MQ=255)                                                                 atgaAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTgcgc                                                                                    >  2:1426955/1‑119 (MQ=255)                                                                      ttCCACGCGGTAGCGGGCGGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAGTTg                                                                                                        <  2:920934/93‑1 (MQ=255)                                                                      gTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTg                                                                                                        >  1:920934/1‑94 (MQ=255)                                                                      gTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTc                                                           >  1:1191140/1‑139 (MQ=255)                                                                      gTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTc                                                           >  1:1196260/1‑139 (MQ=255)                                                                         cACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGAtgt                                                                      <  2:1206102/125‑1 (MQ=255)                                                                         cACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGAtgt                                                                      >  1:1206102/1‑125 (MQ=255)                                                                                        aGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCgg                                         >  2:255143/1‑139 (MQ=255)                                                                                          tAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGc                                             <  2:1320741/133‑1 (MQ=255)                                                                                          tAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGc                                             >  1:1320741/1‑133 (MQ=255)                                                                                            aaCGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTAc                                     <  2:1143548/139‑1 (MQ=255)                                                                                                    gcgcATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCa                             >  2:675449/1‑139 (MQ=255)                                                                                                     cgcATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCaa                            <  2:1243138/139‑1 (MQ=255)                                                                                                       cATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCATTCTGCCACTGTGCGCTCAGCCAGCGATGTGgc                                                                   <  2:564780/98‑1 (MQ=255)                                                                                                       cATCTGGCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCATTCTGCCACTGTGCGCTCAGCCAGCGATGTGgc                                                                   >  1:564780/1‑98 (MQ=255)                                                                                                             gCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGc                                                                                                      <  1:782360/57‑1 (MQ=255)                                                                                                             gCGTTAAGTCTGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGc                                                                                                      >  2:782360/1‑57 (MQ=255)                                                                                                                       tGGAGCGGCTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCAACCCGCCCGCAATGCGCCa          >  2:610603/1‑139 (MQ=255)                                                                                                                               cTCGATGCACTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCAACCCGCCCGCAATGCGCCATCCGCAGa  >  1:897151/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                   CTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCCACACCTCAGCACAGCGAATTGGCCTTTCTGCCACTGTGCGCTCAGCCAGCGATGTGGCAGCATTTCTTGGTGGCGATCGCCCGGTTACGGTAGTCAACCCGCCCGCAATGCGCCATCCGCAGA  >  NC_000913/1729558‑1729822

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