Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,418,632	A→G	Y332H (TAC→CAC)	ulaG ←	L‑ascorbate‑6‑phosphate lactonase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,418,632	0	A	G	100.0%	124.7 / NA	37	Y332H (TAC→CAC)	ulaG	L‑ascorbate‑6‑phosphate lactonase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/18);  total (19/18)

AATACTACGCCGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTG  >  NC_000913/4418507‑4418761                                                                                                                              |                                                                                                                                  aaTACTACGCCGGATTGCGCGACCTTCTCCTCGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAg                                                                                                                            <  1:209212/133‑1 (MQ=255)     ctacGCCGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGt                                                                                                                        >  1:401214/1‑133 (MQ=255)         gCCGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCttt                                                                                                                    >  2:459233/1‑133 (MQ=255)          ccGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTa                                                                                                                   >  2:525161/1‑133 (MQ=255)           cGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTAt                                                                                                                  >  1:607510/1‑133 (MQ=255)                 gcgcGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCg                                                                                                            <  1:471455/133‑1 (MQ=255)                  cgcgAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCgg                                                                                                           >  2:347728/1‑133 (MQ=255)                            tcCACGCTAGCCGGTATAAACGACTGTTTACAGGAATGACTGGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAAt                                                                                                 <  2:288206/133‑1 (MQ=255)                            tcCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAAt                                                                                                 >  2:311843/1‑133 (MQ=255)                                cGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGc                                                                                             <  1:525161/133‑1 (MQ=255)                                       cGGTATAAACGACAGATTACAGGAATGACTTGATCGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGccacccac                                                                                      <  2:401214/133‑1 (MQ=255)                                       cGGTATAAACCACAGTTTACAGGACTGAGTTGAACGGCAGACCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGccacccac                                                                                      <  2:199197/133‑1 (MQ=255)                                           ataAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGc                                                                                  >  2:617245/1‑133 (MQ=255)                                                   cAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATaaa                                                                          <  2:607510/133‑1 (MQ=255)                                                    aGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAAc                                                                         >  2:466247/1‑133 (MQ=255)                                                    aGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAAc                                                                         >  2:138070/1‑133 (MQ=255)                                                                tGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCa                                                             >  2:614134/1‑133 (MQ=255)                                                                tGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCa                                                             >  2:334015/1‑133 (MQ=255)                                                                 gACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCAt                                                            >  2:21751/1‑133 (MQ=255)                                                                  aCTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATa                                                           >  2:55843/1‑133 (MQ=255)                                                                               aTCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCt                                              <  1:617245/133‑1 (MQ=255)                                                                                   ggTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCttttt                                          <  1:334015/133‑1 (MQ=255)                                                                                      tCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCa                                       <  2:490390/133‑1 (MQ=255)                                                                                         aTAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCAtct                                    >  2:313520/1‑133 (MQ=255)                                                                                             tGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCca                                <  2:377985/133‑1 (MQ=255)                                                                                             tGAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCca                                <  1:459233/133‑1 (MQ=255)                                                                                              gAAGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCcac                               >  2:555238/1‑133 (MQ=255)                                                                                               aaGCAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCcaca                              >  1:207160/1‑133 (MQ=255)                                                                                                  cAATCATCGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCa                           >  1:529889/1‑133 (MQ=255)                                                                                                    atcatcGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCAcg                         >  2:82127/1‑133 (MQ=255)                                                                                                         tgAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGAt                    <  1:82127/132‑1 (MQ=255)                                                                                                         cGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGAt                    <  1:555238/133‑1 (MQ=255)                                                                                                         cGAAACCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGAt                    <  1:347728/133‑1 (MQ=255)                                                                                                             aCCGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATctct                <  1:614134/133‑1 (MQ=255)                                                                                                              ccGCGCGGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCtt               <  1:311843/133‑1 (MQ=255)                                                                                                                   cgGATAGTGGTGCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGa          >  2:350268/1‑133 (MQ=255)                                                                                                                           ggTGCTCGAAGTTGTTTTTATCCATCGGGCAGCTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTg  <  2:207160/133‑1 (MQ=255)                                                                                                                              |                                                                                                                                  AATACTACGCCGGATAGCGCGAACTTCTCCACGCTACCCGGTATAAACGACAGATTACAGGAATGACTTGAACGGCAGATCCGGTTCAATAGTGAAGCAATCATCGAAACCGCGCGGATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTG  >  NC_000913/4418507‑4418761

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