Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	262,563	G→A	R354H (CGT→CAT)	proA →	gamma‑glutamylphosphate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	262,563	0	G	A	100.0%	143.7 / NA	46	R354H (CGT→CAT)	proA	gamma‑glutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/19);  total (27/19)

TGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCC  >  NC_000913/262429‑262688                                                                                                                                       |                                                                                                                              tGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCAtttt                                                                                                                           >  2:666992/1‑139 (MQ=255)      aTGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTaa                                                                                                                      <  1:173849/139‑1 (MQ=255)      aTGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTaa                                                                                                                      <  2:188875/139‑1 (MQ=255)       tGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAAc                                                                                                                     >  1:361205/1‑139 (MQ=255)         aGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGa                                                                                                                   >  1:313906/1‑139 (MQ=255)                   aTTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCg                                                                                                         >  2:104969/1‑139 (MQ=255)                   aTTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCg                                                                                                         >  1:891489/1‑139 (MQ=255)                     tAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                       >  2:572323/1‑139 (MQ=255)                     tAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                       <  2:361205/139‑1 (MQ=255)                      aGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTcc                                                                                                      <  2:21508/139‑1 (MQ=255)                      aGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCG‑CCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                     <  1:993096/139‑1 (MQ=255)                          ttGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACTCGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTg                                                                                                  <  1:572323/139‑1 (MQ=255)                           tGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGt                                                                                                 >  1:538181/1‑139 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                     <  2:687584/103‑1 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                     >  1:687584/1‑103 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcg                                                                 >  2:997379/1‑139 (MQ=255)                                                               tCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGTAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGGTTTACCGCcggcgg                                                             >  2:471820/1‑139 (MQ=255)                                                                    cATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCctct                                                                                 >  1:72348/1‑114 (MQ=255)                                                                    cATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCctct                                                                                 <  2:72348/114‑1 (MQ=255)                                                                    cATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGt                                                        >  2:429068/1‑139 (MQ=255)                                                                      tatTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGttt                                                      >  1:338751/1‑139 (MQ=255)                                                                           gTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                 <  1:104969/139‑1 (MQ=255)                                                                            tGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTcc                                                                                                      <  2:750716/85‑1 (MQ=255)                                                                            tGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTcc                                                                                                      >  1:750716/1‑85 (MQ=255)                                                                             gAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTgg                                               <  1:997379/139‑1 (MQ=255)                                                                                 acGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGc                                           >  1:87626/1‑139 (MQ=255)                                                                                   ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGaa                                                                                                                  >  2:558007/1‑66 (MQ=255)                                                                                   ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGaa                                                                                                                  <  1:558007/66‑1 (MQ=255)                                                                                             acTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGt                               <  2:538181/139‑1 (MQ=255)                                                                                                  gATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTa                                                                                >  2:783786/1‑85 (MQ=255)                                                                                                  gATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTa                                                                                <  1:783786/85‑1 (MQ=255)                                                                                                      cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGcacaca                      >  2:64898/1‑139 (MQ=255)                                                                                                         tCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                       <  1:617392/55‑1 (MQ=255)                                                                                                         tCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                       >  2:617392/1‑55 (MQ=255)                                                                                                         tCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTAc                                                                               >  1:310597/1‑79 (MQ=255)                                                                                                         tCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTAc                                                                               <  2:310597/79‑1 (MQ=255)                                                                                                         tCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACaaaa                   >  1:622962/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTaa                                                                                        >  2:422324/1‑60 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTaa                                                                                        <  1:422324/60‑1 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg         >  1:62435/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg         <  2:354271/139‑1 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg         >  1:534558/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGA‑TCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcgc        >  1:540923/1‑139 (MQ=255)                                                                                                                          cgcAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTAGCGGTAAGCACACAAAAACTCCACGCGCGTGGcc  >  1:689992/1‑139 (MQ=255)                                                                                                                                  ccAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACggcggc                                                            <  2:454427/73‑1 (MQ=255)                                                                                                                                  ccAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACggcggc                                                            >  1:454427/1‑73 (MQ=255)                                                                                                                                       |                                                                                                                              TGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCC  >  NC_000913/262429‑262688

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