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

AGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGG  >  NC_000913/262437‑262695                                                                                                                               |                                                                                                                                     aGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCCACGCCCAGCATTTTGTTAACGa                                                                                                                          >  1:1032478/1‑139 (MQ=255)   tttCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAg                                                                                                                        >  1:708434/1‑139 (MQ=255)              aGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTcc                                                                                                             <  2:1015025/139‑1 (MQ=255)               gATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  2:997436/139‑1 (MQ=255)               gATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  1:1260162/139‑1 (MQ=255)                     aaCGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGa                                                                                                                    >  2:546639/1‑125 (MQ=255)                     aaCGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGa                                                                                                                    <  1:546639/125‑1 (MQ=255)                     aaCGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACCCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGCTTCGTCCGCTGttt                                                                                                      <  1:785355/139‑1 (MQ=255)                                   cAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCctct                                                                                        <  1:158156/139‑1 (MQ=255)                                                   gCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcg                                                                        >  1:942890/1‑139 (MQ=255)                                                     cATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGAcggc                                                                      >  1:883010/1‑139 (MQ=255)                                                              tatTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGttt                                                             >  2:630552/1‑139 (MQ=255)                                                                 tCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          <  1:88788/139‑1 (MQ=255)                                                                 tCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          <  1:94364/139‑1 (MQ=255)                                                                  cGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTc                                                         >  1:1156076/1‑139 (MQ=255)                                                                    tGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTg                                                       >  1:670251/1‑139 (MQ=255)                                                                     gAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTgg                                                      >  1:1037692/1‑139 (MQ=255)                                                                        cacGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTg                                                   >  2:542739/1‑139 (MQ=255)                                                                        cacGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTg                                                   >  1:11480/1‑139 (MQ=255)                                                                         acGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            >  2:132978/1‑81 (MQ=255)                                                                         acGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  1:132978/81‑1 (MQ=255)                                                                           ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGtt                                                                                                       >  1:1081135/1‑84 (MQ=255)                                                                           ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGtt                                                                                                       <  2:1081135/84‑1 (MQ=255)                                                                           ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCgg                                                >  2:814136/1‑139 (MQ=255)                                                                           ggCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCgg                                                >  2:1091549/1‑139 (MQ=255)                                                                             cacacaACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGaa                                              <  1:342513/139‑1 (MQ=255)                                                                                       tCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTaa                                    >  2:189609/1‑139 (MQ=255)                                                                                         cGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACg                                                                                                                           >  2:981614/1‑50 (MQ=255)                                                                                         cGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACg                                                                                                                           <  1:981614/50‑1 (MQ=255)                                                                                         cGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGc                                  >  1:741877/1‑139 (MQ=255)                                                                                         cGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGc                                  <  2:1156076/139‑1 (MQ=255)                                                                                              cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAg                                                                <  1:305158/104‑1 (MQ=255)                                                                                              cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAg                                                                >  2:305158/1‑104 (MQ=255)                                                                                              cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          <  1:1110805/110‑1 (MQ=255)                                                                                              cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          >  2:1110805/1‑110 (MQ=255)                                                                                              cGATCCTGACCCGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGcacaca                             >  1:310202/1‑139 (MQ=255)                                                                                                        ccGCGATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAc                   <  2:310202/139‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTc                                                                                                              <  1:886891/45‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            >  2:280505/1‑47 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            >  1:447277/1‑47 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  1:280505/47‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            >  1:48838/1‑47 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  2:447277/47‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCg                                                                                                            <  2:48838/47‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                        >  2:1144703/1‑99 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCt                                                        <  1:1144703/99‑1 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCAcgcg                >  2:237184/1‑139 (MQ=255)                                                                                                           cgATATGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTAGTc                                                                                                              >  2:886891/1‑45 (MQ=38)                                                                                                                tGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          <  1:585616/92‑1 (MQ=255)                                                                                                                tGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGt                                                          >  2:585616/1‑92 (MQ=255)                                                                                                                tGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTgg           >  2:381194/1‑139 (MQ=255)                                                                                                                tGCGCAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTgg           <  1:630552/139‑1 (MQ=255)                                                                                                                    cAACGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCa       >  1:1266118/1‑139 (MQ=255)                                                                                                                     aaCGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCaa      <  2:1266118/139‑1 (MQ=255)                                                                                                                     aaCGCCCAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACACCACGCGCGTGGCCCaa      <  1:809620/139‑1 (MQ=255)                                                                                                                         cccAGCATTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATggg  <  1:237184/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                     AGTTTCTGTCATTAGATTTGAACGTCAAAATCGTCAGCGATCTTGACGATGCCATCGCCCATATTCGTGAACACGGCACACAACACTCCGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGG  >  NC_000913/262437‑262695

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