Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	1,731,621	G→T	M845I (ATG→ATT)	lhr →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	1,731,621	0	G	T	100.0%	104.5 / NA	35	M845I (ATG→ATT)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/19);  total (16/19)

ACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGG  >  NC_000913_3_pae_pgi/1731481‑1731756                                                                                                                                             |                                                                                                                                        aCCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGGACAATGTCTATGCGCCTCCCCCCCTGGAACGTTTGATTCGTCGCCt                                                                                                                                 >  2:150699/1‑149 (MQ=255)   cccGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGcc                                                                                                                                  >  2:690159/1‑146 (MQ=255)   cccGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGcc                                                                                                                                  <  1:690159/146‑1 (MQ=255)            gCAACAACGGCTGCGCGCCAGTCAGATGCTGGAAATCGCTCAGGGATATCCGGATTATACGGTCATTCTCGAAAACCTACGCGAATGTCTGCAAGACGTTTATGATCTTCCCGCCCTGCAACGTTTGATTCGTCGCCTGACAGGTGGCg                                                                                                                      >  2:561189/1‑149 (MQ=255)               acaacGGCGGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                   <  1:239099/149‑1 (MQ=255)                   cGGCTGGGCGGTAGTCAGTTGCTGGATATCGCTCAGGGATATGCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCa                                                                                                               <  2:487193/149‑1 (MQ=255)                          gcgcCAGTCAGTTGCTGGAAATCGCTCAGTGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATcc                                                                                                        <  2:160518/149‑1 (MQ=255)                          gcgcCAGTCAGTTGCTGGAAAGCGCTCAGGGATATCCGGATTGTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATcc                                                                                                        <  1:332345/149‑1 (MQ=255)                           cgcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                   <  1:614874/137‑1 (MQ=255)                           cgcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                   >  2:614874/1‑137 (MQ=255)                           cgcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTCTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                   <  1:256326/137‑1 (MQ=255)                           cgcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTCTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                   >  2:256326/1‑137 (MQ=255)                              cagtcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATg                                                                                                    >  1:515151/1‑149 (MQ=255)                                    gatGCTGGAAATCGCTGAGGGATATCCGGCGTTCCCGGTCATTTCGGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCATTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                                                                              <  2:442064/147‑1 (MQ=255)                                                    cAGGGATATCCGGGTTTTCCCGTCATTCTCGAAACCCTACGCTAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGcc                                                                              <  1:124230/149‑1 (MQ=255)                                                     ggggATATCCGGATTTTCCGGTCATTATCGAAAGCCTACGCGAATGTCTGCAAGATGTTTATGAGCTTCCCGCACTGCAACGTTTGGTTCGTCGCCTGAGCGGTGGCGAAATTCAACTATCCTATGTACAGACCACTACGGCCTCGCCt                                                                             >  2:389102/2‑149 (MQ=255)                                                                     tCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAGGATGTTTTTTAGCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAt                                                             <  1:134444/149‑1 (MQ=255)                                                                         gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATCTTTATGATCTTCCCGAAAATCAACGGGTGATTCGTCGCCTGAACGGTGGCGAAAATCAAATATCCGATGTAACGACCACAACGCCCTCTCCTTTCGCCACAAGTTGATTGtt                                                         >  1:386597/1‑149 (MQ=255)                                                                         gTCAGTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt                                                         <  2:192231/149‑1 (MQ=255)                                                                          tCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc                                                        >  2:474117/1‑149 (MQ=255)                                                                           cATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg                                                       >  2:200790/1‑149 (MQ=255)                                                                            aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg                                                      >  2:71268/1‑149 (MQ=255)                                                                              tctcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCt                                                    <  1:227819/149‑1 (MQ=255)                                                                                tcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAt                                                  <  2:515151/149‑1 (MQ=255)                                                                                       ccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTGGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCgg                                           <  2:305339/149‑1 (MQ=255)                                                                                             cgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGCACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGTCCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTa                                     <  1:601182/149‑1 (MQ=255)                                                                                                     tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTa                                                              <  2:316469/116‑1 (MQ=255)                                                                                                     tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTAGCCACAAGTTTa                                                              >  1:316469/1‑116 (MQ=255)                                                                                                        aaGATGTTTATGATCTTCCCGCACTAGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATACGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGc                          >  1:547812/1‑149 (MQ=255)                                                                                                          gATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGACAGGTGGCGAAATTCCAATATCCGATGTAACGACCCCTACGCCCTCGCCTTTCGCCCCAAGTTTCTTCTTCGGCTATACCGCGACATTTACTTCCCAGCGCGa                        >  2:319543/1‑149 (MQ=255)                                                                                                                tATGATCTTCCCGCACTGGAAGGTTTGATTCGTCGCCTGAACCGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCCATGTCGCGGAATTTATGTACCAGAGCGAcccccc                  >  2:632729/1‑149 (MQ=255)                                                                                                                tATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGcccc                  >  1:162910/1‑149 (MQ=255)                                                                                                                        tCCCGCACTGGATCGTTTGATTCGTCGCCTGAACGGTGGCGAAAGTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCag          <  2:289591/149‑1 (MQ=255)                                                                                                                               cTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCg   <  1:550430/149‑1 (MQ=255)                                                                                                                                tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCTGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCgg  <  1:474117/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                        ACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGG  >  NC_000913_3_pae_pgi/1731481‑1731756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑