Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	113,502	G→A	R20H (CGC→CAC)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	113,502	0	G	A	100.0%	126.7 / NA	37	R20H (CGC→CAC)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/16);  total (21/16)

CATTGATTAGTGCGTATGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCG  >  NC_000913/113368‑113638                                                                                                                                       |                                                                                                                                         cATTGATTAGTGCGTATGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCa                                                                                                                                     >  1:65433/1‑140 (MQ=255)        tAGTGCGTATGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTaa                                                                                                                              <  2:799258/140‑1 (MQ=255)               tATGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGt                                                                                                                       >  1:737470/1‑140 (MQ=255)                aTGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGtt                                                                                                                      >  1:154300/1‑140 (MQ=255)                aTGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGtt                                                                                                                      >  2:478245/1‑140 (MQ=255)                     aGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGAt                                                                                                                 >  1:457369/1‑140 (MQ=255)                      gCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATg                                                                                                                <  2:404810/140‑1 (MQ=255)                         tCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTg                                                                                                             >  2:40719/1‑140 (MQ=255)                          cACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGa                                                                                                            >  1:194315/1‑140 (MQ=255)                           aCTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGaa                                                                                                           >  1:656879/1‑140 (MQ=255)                                  ttGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAAc                                                                                                    >  1:580244/1‑140 (MQ=255)                                         ctTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAAtt                                                                                             <  1:525919/140‑1 (MQ=255)                                                  ataGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCaa                                                                                    >  1:627936/1‑140 (MQ=255)                                                  ataGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCaa                                                                                    >  2:411253/1‑140 (MQ=255)                                                    aGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAAc                                                                                  >  2:534581/1‑140 (MQ=255)                                                       cATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACAtt                                                                               <  1:521797/140‑1 (MQ=255)                                                       cATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACAtt                                                                               <  2:656879/140‑1 (MQ=255)                                                       cATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACAtt                                                                               <  2:457369/140‑1 (MQ=255)                                                        aTTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTc                                                                              >  1:45872/1‑140 (MQ=255)                                                        aTTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTc                                                                              >  1:699875/1‑140 (MQ=255)                                                                            caTGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCg                                                          <  1:478245/140‑1 (MQ=255)                                                                             aTGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGt                                                         >  1:757895/1‑140 (MQ=255)                                                                                  tATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGa                                                    >  1:130616/1‑140 (MQ=255)                                                                                         gaagaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATcgc                                             >  1:91533/1‑140 (MQ=255)                                                                                           agaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATcgccg                                           >  2:723060/1‑140 (MQ=255)                                                                                            gaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATcgccgc                                          >  1:709249/1‑140 (MQ=255)                                                                                              tCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCaa                                        <  2:613515/140‑1 (MQ=255)                                                                                              tCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCaa                                        <  1:452209/140‑1 (MQ=255)                                                                                              tCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCaa                                        >  1:10077/1‑140 (MQ=255)                                                                                                tGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAat                                      <  2:627936/140‑1 (MQ=255)                                                                                                          gttgaaaGACGTTCTCAGCGGCCCTAACCGCGCCGCGCTCAAAAGCCGTTCCGATGTTGACCTGGAACGTCAATTCCCCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCg                            <  2:339565/136‑1 (MQ=255)                                                                                                              aaaGACGTTCTCATCCGCCCTCAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAgg                        <  2:737470/140‑1 (MQ=255)                                                                                                              aaaGACGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAgg                        <  2:10077/140‑1 (MQ=255)                                                                                                                  aCGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGcaca                    <  2:65433/140‑1 (MQ=255)                                                                                                                   cGTTCTCATCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACAt                   <  2:709249/140‑1 (MQ=255)                                                                                                                          aTCCGCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTAt            <  1:723060/140‑1 (MQ=255)                                                                                                                              gCCCTAAACACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGcc        >  2:433283/1‑140 (MQ=255)                                                                                                                                    aaCACTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTgcg  <  2:757895/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         CATTGATTAGTGCGTATGATAGCGTCACTGGAGTTGCGCTCTTACCCTTATAGCCATTAACCCCAGGAATCCGCACATGCGTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCG  >  NC_000913/113368‑113638

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