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%	163.0 / NA	51	R354H (CGT→CAT)	proA	gamma‑glutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/24);  total (27/24)

TGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAAT  >  NC_000913/262429‑262692                                                                                                                                       |                                                                                                                                  tGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCAtttt                                                                                                                               >  1:1296053/1‑139 (MQ=255)     gATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTa                                                                                                                           <  2:334372/139‑1 (MQ=255)        gAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACg                                                                                                                        <  2:116183/139‑1 (MQ=255)         aGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGa                                                                                                                       >  2:1192307/1‑139 (MQ=255)         aGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGa                                                                                                                       >  1:553057/1‑139 (MQ=255)             tCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTg                                                                                                                   <  2:354721/139‑1 (MQ=255)                     tAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                           >  1:898284/1‑139 (MQ=255)                     tAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACAGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                           <  2:16766/139‑1 (MQ=255)                      aGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTcc                                                                                                          <  1:1192307/139‑1 (MQ=255)                                           cAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCctct                                                                                     <  2:1296053/139‑1 (MQ=255)                                            aGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTa                                                                                    >  2:1090714/1‑139 (MQ=255)                                                          tGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGaa                                                                                                                      >  1:667353/1‑91 (MQ=255)                                                          tGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGaa                                                                                                                      <  2:667353/91‑1 (MQ=255)                                                          tGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAc                                                                      >  1:976355/1‑139 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGttt                                                                                                   >  1:653204/1‑109 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGttt                                                                                                   <  2:653204/109‑1 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcg                                                                     >  1:1185024/1‑139 (MQ=255)                                                           gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcg                                                                     >  1:532478/1‑139 (MQ=255)                                                               tCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcggcgg                                                                 <  1:1090714/139‑1 (MQ=255)                                                                         tCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                       <  2:553057/139‑1 (MQ=255)                                                                             gAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGtt                                                                                                    <  1:591246/90‑1 (MQ=255)                                                                             gAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGtt                                                                                                    >  2:591246/1‑90 (MQ=255)                                                                               acacGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGt                                                 >  1:974240/1‑139 (MQ=255)                                                                                 acGGCACACAACTCTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGc                                               <  2:526845/139‑1 (MQ=255)                                                                                  cGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCAttt                                                                                                                                >  2:805510/1‑57 (MQ=255)                                                                                  cGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCAttt                                                                                                                                <  1:805510/57‑1 (MQ=255)                                                                                   ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTa                                                                                                                           >  1:150297/1‑61 (MQ=255)                                                                                   ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTa                                                                                                                           <  2:150297/61‑1 (MQ=255)                                                                                               tCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTaa                                 >  1:206905/1‑139 (MQ=255)                                                                                                  gATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGca                              <  1:1040952/139‑1 (MQ=255)                                                                                                  gATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGca                              <  2:976355/139‑1 (MQ=255)                                                                                                    tGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGcaca                            >  2:1250372/1‑139 (MQ=255)                                                                                                      cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTaa                                                                                                                          >  2:184905/1‑43 (MQ=255)                                                                                                      cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTaa                                                                                                                          <  1:184905/43‑1 (MQ=255)                                                                                                        aTCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACaaa                        <  2:206905/139‑1 (MQ=255)                                                                                                             gACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTg                                                                                                                   >  2:239147/1‑43 (MQ=255)                                                                                                             gACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTg                                                                                                                   <  1:239147/43‑1 (MQ=255)                                                                                                              aCCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                     <  2:1014751/104‑1 (MQ=255)                                                                                                              aCCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                     >  1:1014751/1‑104 (MQ=255)                                                                                                              aCCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTcc                  >  2:982880/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                     >  1:8648/1‑99 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                     <  2:8648/99‑1 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGaa                                           >  1:964335/1‑109 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGaa                                           <  2:964335/109‑1 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg             >  1:430166/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg             >  1:872430/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg             >  2:699283/1‑139 (MQ=255)                                                                                                                   cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg             >  2:774267/1‑139 (MQ=255)                                                                                                                     atatGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcgcg           >  2:462455/1‑139 (MQ=255)                                                                                                                     atatGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcgcg           <  1:982880/139‑1 (MQ=255)                                                                                                                              aCGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAAt  <  2:1185024/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  TGACGATGAGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAAT  >  NC_000913/262429‑262692

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