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%	98.9 / NA	31	M845I (ATG→ATT)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/15);  total (16/15)

ACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGC  >  NC_000913_3_pae_pgi/1731481‑1731758                                                                                                                                             |                                                                                                                                          aCCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCt                                                                                                                                   >  1:234842/1‑149 (MQ=255) aCCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCGCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCt                                                                                                                                   >  1:198650/1‑149 (MQ=255)          tGGCAACAACGGCTGGGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTgg                                                                                                                          <  1:293730/149‑1 (MQ=255)            gCAACAGCGGCGGCCCGCCAGTCAGTTGCTGGAAATCGCTCAGTGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGTTTCGTCGCCTGAACGGTGGCg                                                                                                                        <  2:610445/149‑1 (MQ=255)            gCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCg                                                                                                                        >  2:232237/1‑149 (MQ=255)              aacaacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGGGGGGaa                                                                                                                      >  1:109501/1‑149 (MQ=255)                caacGGTTGGGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAt                                                                                                                    <  2:102338/149‑1 (MQ=255)                 aacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTCATTCGTCGCCTGAACGGTGGCGAAAtt                                                                                                                   >  1:671307/1‑149 (MQ=255)                         tccgtccgtcagtTGCTGGAAATCGATCAGGGATATCCTGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATTTATATGATCTTCCCGCCCTGGATCCTCTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATc                                                                                                           <  1:20215/142‑1 (MQ=255)                         cgcgcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATc                                                                                                           <  1:232237/149‑1 (MQ=255)                           cgccgcgcagcTGCTGGAGATCGGTGAGGGATAGCCGGATGTGCCGGTCAGTCTGGAAACTCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGCAACGTTTGATTCGTCGCTTGAACGGTGGCGAAATTCAAATATCCg                                                                                                         <  2:369508/142‑1 (MQ=255)                                  cagtTGCTGGAAATCGCTCAGGGATAGCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTGTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAAc                                                                                                  <  1:211123/149‑1 (MQ=255)                                       gCTGGAAATCGCTCACGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGCTTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                                                                             <  1:218478/149‑1 (MQ=255)                                         tGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACt                                                                                           >  2:676078/1‑149 (MQ=255)                                              aTCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTgacgac                                                                                               >  2:683205/1‑140 (MQ=255)                                              aTCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAc                                                                                               <  1:683205/140‑1 (MQ=255)                                                             ccGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCcac                                                                       >  2:383537/1‑149 (MQ=255)                                                                                         tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGaa                                           >  2:480378/1‑149 (MQ=255)                                                                                                aTGTCTTCAAGATGTTTATGATCTTCCCGGAATGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGt                                    <  2:478996/149‑1 (MQ=255)                                                                                                aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCGCAAGGTTATTGTTCGGCTATGTCGCGGAATTTATGt                                    <  2:497005/149‑1 (MQ=255)                                                                                                aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAAACATCTCCCGGCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGt                                    <  2:115050/149‑1 (MQ=255)                                                                                                      gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCaga                              <  2:112113/149‑1 (MQ=255)                                                                                                                tATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATCTCCGCTGTAACTACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGcccc                    >  1:131881/1‑149 (MQ=255)                                                                                                                tATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGcccc                    >  1:537357/1‑149 (MQ=255)                                                                                                                tATGATCTTCACGCACTGGAACGTTTGAGTAGTCGCCTTCACGATGGGGAAATTCAAAGGTCCTAGGTGACCAACCCTACGACCACGCCTTTGGCAACAAGTTTATTGTTCGCCCATGTCGCGGAATTTATGTACACGAGCCACgagcg                    >  2:608744/1‑145 (MQ=255)                                                                                                                     tCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGCTGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGCGCGACGCCCCGCttc               >  2:148404/1‑147 (MQ=255)                                                                                                                     acctCCCGCACTGGGACGTTTGATTCGTCGCCTGAACTGTGGCGAAATTCAAATATCCGATGTACCGACCACTCCGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTgg               <  1:329371/146‑1 (MQ=255)                                                                                                                       ttCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCCGAGCGACGCCCCCCTgacg             >  2:5731/1‑146 (MQ=255)                                                                                                                           cGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAgc         >  1:375292/1‑149 (MQ=255)                                                                                                                                 ggAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCggg   <  1:480378/149‑1 (MQ=255)                                                                                                                                  gAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGc  >  2:339332/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                          ACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGC  >  NC_000913_3_pae_pgi/1731481‑1731758

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: ‑