Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,336,073	C→G	G141A (GGT→GCT)	murA ←	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,336,073	0	C	G	100.0%	226.4 / NA	56	G141A (GGT→GCT)	murA	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (32/24);  total (32/24)

TCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGA  >  NC_000913/3335926‑3336218                                                                                                                                                    |                                                                                                                                                  tCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGct                                                                                                                                                 <  1:144004/150‑1 (MQ=255)    ggTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACAGGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGcttctt                                                                                                                                             <  1:362234/151‑1 (MQ=255)     gTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCttcttc                                                                                                                                            >  2:328309/1‑151 (MQ=255)      ttCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTcc                                                                                                                                           >  1:808935/1‑151 (MQ=255)        cACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAg                                                                                                                                         <  1:755293/151‑1 (MQ=255)           gcgcTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGttt                                                                                                                                      >  2:508509/1‑151 (MQ=255)               tGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATg                                                                                                                                  >  1:415039/1‑151 (MQ=255)                 cGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGAAAGCTTCTTCCAGTTTGATGGt                                                                                                                                <  1:294773/151‑1 (MQ=255)                              cGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTg                                                                                                                    <  2:893920/150‑1 (MQ=255)                                    gCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAgg                                                                                                             <  2:17587/151‑1 (MQ=255)                                    gCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAgg                                                                                                             >  2:808639/1‑151 (MQ=255)                                      cTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGcc                                                                                                           >  1:269206/1‑151 (MQ=255)                                           gCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGaa                                                                                                       >  2:594992/1‑150 (MQ=255)                                             cAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtg                                                                                                    >  2:630057/1‑151 (MQ=255)                                              aGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtgt                                                                                                   >  1:585983/1‑151 (MQ=255)                                               gggTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtgtg                                                                                                  >  2:232067/1‑151 (MQ=255)                                                ggTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATgtgt                                                                                                 <  2:587689/151‑1 (MQ=255)                                                   tGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAg                                                                                               <  1:635626/150‑1 (MQ=255)                                                     cagcaCACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATgtgt                                                                                                 >  2:186126/1‑146 (MQ=255)                                                         acacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATc                                                                                            >  2:510960/1‑147 (MQ=255)                                                          cacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCg                                                                                       <  1:453260/151‑1 (MQ=255)                                                          cacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCg                                                                                       >  1:1030658/1‑151 (MQ=255)                                                          cacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCg                                                                                       <  1:280865/151‑1 (MQ=255)                                                          cacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCg                                                                                       <  2:441716/151‑1 (MQ=255)                                                               gatgGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGAcgc                                                                                  >  2:693813/1‑151 (MQ=255)                                                                           tGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGt                                                                       >  1:838199/1‑150 (MQ=255)                                                                            gcgcCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTa                                                                      <  1:616427/150‑1 (MQ=255)                                                                              gcCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACaa                                                                   <  1:762823/151‑1 (MQ=255)                                                                               ccAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAAc                                                                  >  2:1065126/1‑151 (MQ=255)                                                                                 aaCGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCg                                                                <  1:666139/151‑1 (MQ=255)                                                                                 aaCGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCg                                                                >  1:905294/1‑151 (MQ=255)                                                                                    gCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGcc                                                              >  2:503142/1‑150 (MQ=255)                                                                                      tGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAgg                                                           >  2:457139/1‑151 (MQ=255)                                                                                       gACTTTTTCCTTCCCAATATGTGCACGTTTCAAACAACCACCGACGGAGGCTTTACCGGAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTACAAGCCCAGAAGGGTGTAGATCAGCCGGACGCGCACCGATGGTAAAACCGCCAgg                                                           <  2:463495/150‑1 (MQ=255)                                                                                       gACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGt                                                          <  1:992791/151‑1 (MQ=255)                                                                                       gACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGt                                                          <  1:554783/151‑1 (MQ=255)                                                                                       gACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGt                                                          >  1:537396/1‑151 (MQ=255)                                                                                           ttATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGt                                                       <  1:604153/150‑1 (MQ=255)                                                                                           ttATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTg                                                      >  1:905455/1‑151 (MQ=255)                                                                                                 aTCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACtt                                                <  1:43471/151‑1 (MQ=255)                                                                                                   cACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGc                                              >  2:811166/1‑151 (MQ=255)                                                                                                              cACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCaa                                    >  2:749439/1‑150 (MQ=255)                                                                                                                 cTTTCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAgcg                                >  2:860642/1‑151 (MQ=255)                                                                                                                  tttCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAgcgc                               >  1:784187/1‑151 (MQ=255)                                                                                                                    tCAAACGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAgcgcg                              <  2:615706/150‑1 (MQ=255)                                                                                                                      aaaCGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCt                            >  1:230553/1‑150 (MQ=255)                                                                                                                       aaCGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTG‑CCCTGACCAAAGCGCGCTAcc                         >  2:1053150/1‑151 (MQ=255)                                                                                                                        aCGACCATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTAcc                         <  1:112423/151‑1 (MQ=255)                                                                                                                             cATCGACGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCgg                    >  1:598863/1‑151 (MQ=255)                                                                                                                                   cGGAAGCTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGcccc                >  1:899344/1‑149 (MQ=255)                                                                                                                                         cTTTAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGccc         >  1:281645/1‑150 (MQ=255)                                                                                                                                           ttAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAg       <  2:479344/150‑1 (MQ=255)                                                                                                                                            tAACGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGACCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGa      <  2:1010251/150‑1 (MQ=255)                                                                                                                                             aaCGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGCTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGAt     <  2:966231/150‑1 (MQ=255)                                                                                                                                               cGTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGAt     >  1:94621/1‑148 (MQ=255)                                                                                                                                                gTAAGCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGa  >  1:948286/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGA  >  NC_000913/3335926‑3336218

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