Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,206,545	T→G	T330P (ACC→CCC)	purH ←	bifunctional AICAR transformylase/IMP cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,206,545	0	T	G	100.0%	174.7 / NA	50	T330P (ACC→CCC)	purH	bifunctional AICAR transformylase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/25);  total (25/25)

GCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGTTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAA  >  NC_000913/4206420‑4206679                                                                                                                              |                                                                                                                                       gCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCt                                                                                                                          <  1:847401/140‑1 (MQ=255) gCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCt                                                                                                                          <  2:436175/140‑1 (MQ=255) gCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCt                                                                                                                          <  1:566654/140‑1 (MQ=255) gCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCt                                                                                                                          >  2:579811/1‑140 (MQ=255)   cccACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTcg                                                                                                                        <  1:111423/140‑1 (MQ=255)        tGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGtt                                                                                                                   <  1:108814/140‑1 (MQ=255)             gCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAgg                                                                                                              >  2:218712/1‑140 (MQ=255)                ggTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCaa                                                                                                           <  1:218712/140‑1 (MQ=255)                 gTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAat                                                                                                          >  2:550467/1‑140 (MQ=255)                       aCGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATgcc                                                                                                    <  1:109982/140‑1 (MQ=255)                         gcgTACGTTGTGTTTGGCGGCGGTGATTTTCAGGGTTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTTGCGGTTAAAGGCAATGATgccgc                                                                                                  <  2:459286/140‑1 (MQ=255)                         gcgTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATgccgc                                                                                                  <  2:580854/140‑1 (MQ=255)                           gTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGccgccg                                                                                                >  2:744975/1‑140 (MQ=255)                                ttCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAatga                                                                                                        >  2:845644/1‑127 (MQ=255)                                ttCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAatga                                                                                                        <  1:845644/127‑1 (MQ=255)                                    gTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAg                                                                                       >  2:597415/1‑140 (MQ=255)                                          cggcggTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTggg                                                                                 <  1:597415/140‑1 (MQ=255)                                             cggTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCg                                                                              >  2:611251/1‑140 (MQ=255)                                                       cAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTAcg                                                                    <  1:917666/140‑1 (MQ=255)                                                          ggCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTAcgcgc                                                                 >  1:924061/1‑140 (MQ=255)                                                             ttcttcGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGTTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGAt                                                              >  1:122183/1‑140 (MQ=255)                                                               cttcGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCg                                                            >  2:844460/1‑140 (MQ=255)                                                                ttcGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGt                                                           <  2:157485/140‑1 (MQ=255)                                                                       gCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATc                                                    >  2:70867/1‑140 (MQ=255)                                                                               cgcAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATgg                                            >  2:387295/1‑140 (MQ=255)                                                                               cgcAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATgg                                            >  2:469433/1‑140 (MQ=255)                                                                                 cAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGaa                                          <  2:924061/140‑1 (MQ=255)                                                                                           ttCGACAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAg                                >  2:63239/1‑140 (MQ=255)                                                                                                cAAACTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCAcg                           >  2:305588/1‑140 (MQ=255)                                                                                                    cTGACGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCAcgccgc                       <  2:926537/140‑1 (MQ=255)                                                                                                       aCGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCa                             <  1:820998/131‑1 (MQ=255)                                                                                                       aCGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCa                             >  2:820998/1‑131 (MQ=255)                                                                                                       aCGAGAAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAg                    >  1:490636/1‑140 (MQ=255)                                                                                                         gagaAATGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAggg                  >  2:855156/1‑140 (MQ=255)                                                                                                            aaaTGATGGCCTGTGCGGGTTCCGCATCCAGTTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTg               <  2:122183/140‑1 (MQ=255)                                                                                                            aaaTGATGGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTg               <  1:70867/140‑1 (MQ=255)                                                                                                               tgatgGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCg            <  2:428252/140‑1 (MQ=255)                                                                                                               tgatgGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCg            <  1:844460/140‑1 (MQ=255)                                                                                                                 atgGCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTg          <  1:305588/140‑1 (MQ=255)                                                                                                                    gCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCtt       <  1:469433/140‑1 (MQ=255)                                                                                                                    gCCTGTGCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCtt       <  1:387295/140‑1 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCg            >  1:243099/1‑130 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCg            <  2:243099/130‑1 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCtt       >  2:276686/1‑135 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCtt       <  1:276686/135‑1 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTca     >  1:655738/1‑137 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTca     <  2:655738/137‑1 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACaa  >  1:483944/1‑140 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACaa  >  1:40147/1‑140 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACaa  >  1:357777/1‑140 (MQ=255)                                                                                                                         tgCGGGTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGCCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACaa  >  2:657213/1‑140 (MQ=255)                                                                                                                              |                                                                                                                                       GCCCCACTGACCGCAGGTCAGAACGCGTACGTTCTGTTTGGCGGCGGTGATTTTCAGGGCTTCTTCGCTGGCGGACGGCGCAATAATCACTTCGACAAACTGACGAGAAATGATGGCCTGTGCGGTTTCCGCATCCAGCTCGCGGTTAAAGGCAATGATGCCGCCGAATGCGGAGGTTGGGTCGGTTTTGTACGCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAA  >  NC_000913/4206420‑4206679

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