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

CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTA  >  NC_000913_3_pae_pgi/1731475‑1731765                                                                                                                                                   |                                                                                                                                                 cATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTGCCCGCACTGGAACGTTTGATTCGt                                                                                                                                               <  1:741742/150‑1 (MQ=255)          ccGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGCTTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTggag                                                                                                                                      >  1:25981/1‑147 (MQ=255)                    aacaacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                            >  1:545421/1‑150 (MQ=255)                            cTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAACCCCTACGCGAATGTCTGCCAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTGTCCTGACCGGTGGCCAAAtcgatata                                                                                                                    >  2:540811/1‑143 (MQ=255)                                           ttGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAcc                                                                                                     >  1:333760/1‑150 (MQ=255)                                             gCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                                                                                   <  2:803076/150‑1 (MQ=255)                                                             ggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTc                                                                                   <  2:545421/150‑1 (MQ=255)                                                               atatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACAACTACGCCCTCGCCTTTCGc                                                                                 >  1:407876/1‑150 (MQ=255)                                                                                tCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGGTTTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg                                                                <  1:928951/150‑1 (MQ=255)                                                                                 cATTCTCGAACCCCGCCGGGGATGACTGCAAGATGTTTATGATCTTCCCGCAATGGTACTTTTGCTTGCTGGCCTGAACCGTGGCGAACTTCAAACATCCGATGTATCGACCGGTTCGCCCTGGCCTTGCGCCACGAGTTTACTGTTCgg                                                               <  2:445913/150‑1 (MQ=255)                                                                                  aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGc                                                              >  1:275772/1‑150 (MQ=255)                                                                                  aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCATCGCCTGAACGGCGGCGAAATTCAAATATACGATGTACCGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGGTCGGc                                                              >  1:828466/1‑150 (MQ=255)                                                                                     ctcGCAACCCTACGCGAATGCACGCAATATTTTTATGATCTTCCCGCACTGGAACGGTTGATTCGGCGCCTGAACGGTGGCGCAATTCAAATATCCGATGTAACGACCACTACTCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAt                                                           <  2:682895/150‑1 (MQ=255)                                                                                      tcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATg                                                          >  2:543494/1‑150 (MQ=255)                                                                                      tcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATg                                                          >  1:200695/1‑150 (MQ=255)                                                                                         aaaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTGTGTcg                                                       >  1:166152/1‑150 (MQ=255)                                                                                          aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTTGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                                                                            <  2:979837/128‑1 (MQ=255)                                                                                          aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                                                                            >  1:979837/1‑128 (MQ=255)                                                                                          aaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTcgc                                                      <  2:333760/150‑1 (MQ=255)                                                                                              cTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGaa                                                  >  2:294367/1‑150 (MQ=255)                                                                                                   cgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTAt                                             <  1:56164/150‑1 (MQ=255)                                                                                                           tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAATCAAATATCCGATGTACCGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTAGGTACCaga                                     >  1:190866/1‑150 (MQ=255)                                                                                                               aGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAGCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGa                                 >  2:105886/1‑150 (MQ=255)                                                                                                                       ttGATTTTCCCGCACTGGAACGTTGGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGc                         <  2:144248/149‑1 (MQ=255)                                                                                                                        tGATTTGCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCt                        <  2:863790/150‑1 (MQ=255)                                                                                                                        tGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCAAAAGTTTATTGTTCGGCTATGTCGCGGACTTTCTGTACCAGAGCGACGCCCCCCt                        >  2:154820/1‑150 (MQ=255)                                                                                                                            tttCCCGCACTGGATAGTTTGATTCGCGCCCTGAACGGTGGCGAAATTCAAATATGCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAATTTTATTGTTCGGCTATGTCGCGGAATTTATGTACGAGAGCGACGCCCCGCTGGCa                    <  2:9150/149‑1 (MQ=255)                                                                                                                            cgtCCCGCACTGGAACGTTGGATTCGTCGCTTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCa                    <  1:372559/148‑1 (MQ=255)                                                                                                                             ttCCCGCACTGGAACGTTTGATTCGTCGCATGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTAGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTgcccg                   >  2:892217/1‑146 (MQ=255)                                                                                                                                 cGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAgcg               >  2:188347/1‑150 (MQ=255)                                                                                                                                        gAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCCGAGCGCCGGGCa        >  2:769669/1‑150 (MQ=255)                                                                                                                                              ttGATTCGTCGCCTGAACGGGGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTa  <  2:873708/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTA  >  NC_000913_3_pae_pgi/1731475‑1731765

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