Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,096,038	G→A	W65* (TGG→TGA)	hisF →	imidazole glycerol phosphate synthase subunit HisF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,096,038	0	G	A	100.0%	125.5 / NA	37	W65* (TGG→TGA)	hisF	imidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

CGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACC  >  NC_000913/2095903‑2096171                                                                                                                                        |                                                                                                                                       cGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTAt                                                                                                                                    >  2:539297/1‑140 (MQ=255)  gTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTAtc                                                                                                                                   <  2:576167/140‑1 (MQ=255)   tACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTAtct                                                                                                                                  >  2:356875/1‑140 (MQ=255)          tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCAATGGCCGTGTGGTAGATAAAATTTGAGTATCTCGCGTgg                                                                                                                           >  2:373406/1‑140 (MQ=255)             cAACCATGAAATCATCGGCGATTTCGTGCGGCTGGCATAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCgg                                                                                                                        <  1:318801/140‑1 (MQ=255)                 cATGAAATCATTGGCGATATCGAGCTGCTAGCAAAACTCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGt                                                                                                                    <  1:328854/140‑1 (MQ=255)                  aTGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTcg                                                                                                                                <  1:851291/127‑1 (MQ=255)                  aTGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTcg                                                                                                                                >  2:851291/1‑127 (MQ=255)                  aTGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTg                                                                                                                   >  1:470871/1‑140 (MQ=255)                  aTGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTg                                                                                                                   >  2:482260/1‑140 (MQ=255)                                        gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTgg                                                                                             >  1:3701/1‑140 (MQ=255)                                           gCTGGCAAAACGCTACGCTGAAGAAGGCGATGACGAATTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCgg                                                                                          <  1:373406/140‑1 (MQ=255)                                            cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCggg                                                                                         >  2:616768/1‑140 (MQ=255)                                            cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCggg                                                                                         <  2:769907/140‑1 (MQ=255)                                             tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt                                                                                        <  2:669369/140‑1 (MQ=255)                                              ggCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTg                                                                                       >  1:522285/1‑140 (MQ=255)                                               gCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCgggggg                                                                                      >  2:508837/1‑140 (MQ=255)                                                cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTggg                                                                                     >  1:694240/1‑140 (MQ=255)                                                   aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGAtt                                                                                  <  1:356875/140‑1 (MQ=255)                                                         acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtct                                                                            <  2:694240/140‑1 (MQ=255)                                                                        cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCcgcg                                                            <  1:508837/140‑1 (MQ=255)                                                                         tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGGCTCTGGAAAAT‑GCCGCGa                                                           >  1:810888/1‑140 (MQ=255)                                                                         tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGTGGGTTAAGACTCTGGTAGAG‑Ggcgcca                                                           >  1:61911/1‑138 (MQ=255)                                                                         tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGAA‑GCCGCGa                                                           >  1:34628/1‑140 (MQ=255)                                                                            cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGaaaa                                                        >  1:593890/1‑140 (MQ=255)                                                                                 tGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCtt                                                   <  1:616768/140‑1 (MQ=255)                                                                                          aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGgc                                          <  2:593890/140‑1 (MQ=255)                                                                                               atCACCGCTTCCAGCCATGGCCCTGGGGTCGATAAAAACTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGa                                     <  1:244114/140‑1 (MQ=255)                                                                                                   ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAA‑GCCGCGGAAATTCTTTTCTTTGGCGGGGGTaaa                                 >  2:98192/1‑140 (MQ=255)                                                                                                        tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATATCTTG‑AGATAGCGCCCGAAATTCTTATCTTATGCGCGGGTAAAAtttt                            >  2:244114/1‑139 (MQ=255)                                                                                                            gCGATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCTGGAAGAT‑GCCGCGAAAATTTTTTTCTTTGGCGCGGGTAAAATTTTCATc                        >  2:633482/1‑140 (MQ=255)                                                                                                              gATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                      <  1:502924/140‑1 (MQ=255)                                                                                                              gATGGCCGTGTGGTAGATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                      <  1:98192/140‑1 (MQ=255)                                                                                                                             gATAAAAGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGc       <  1:456510/140‑1 (MQ=255)                                                                                                                                aaaaGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGa    <  2:3701/140‑1 (MQ=255)                                                                                                                                aaaaGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGa    <  2:34628/140‑1 (MQ=255)                                                                                                                                  aaGCTGAGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAcc  <  1:633482/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACC  >  NC_000913/2095903‑2096171

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