Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	713,562	C→T	A2V (GCA→GTA)	pgm →	phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	713,562	0	C	T	100.0%	103.5 / NA	30	A2V (GCA→GTA)	pgm	phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/12);  total (18/12)

GATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGCAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAATATTAT  >  NC_000913/713484‑713638                                                                               |                                                                             gATGCAACTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTaa                                                                            >  1:134514/1‑81 (MQ=255)      aaTTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCcac                                                                       <  1:134493/81‑1 (MQ=255)      aaTTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCca                                                                        >  1:134535/1‑80 (MQ=255)       aTTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCcac                                                                       >  1:134524/1‑80 (MQ=255)             ggcggAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTg                                                                >  1:134517/1‑81 (MQ=255)               cggAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCa                                                              <  2:134494/81‑1 (MQ=255)                 gAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAgg                                                            >  2:134527/1‑81 (MQ=255)                   attgattgAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGc                                                           >  2:134523/1‑80 (MQ=255)                       attgaGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAAcc                                                      >  2:134529/1‑81 (MQ=255)                        ttgaGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCt                                                     <  2:134495/81‑1 (MQ=255)                          gagaAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTg                                                    >  2:134538/1‑80 (MQ=255)                                           tATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGtgatttga                                  >  2:134537/1‑81 (MQ=255)                                             tCTAAAACGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGtgatttgat                                 >  2:134526/1‑80 (MQ=255)                                                 aaaCGTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTa                               <  1:134496/78‑1 (MQ=255)                                                     gTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAAcgccgc                        >  2:134528/1‑81 (MQ=255)                                                     gTTGCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGc                        >  1:134533/1‑81 (MQ=255)                                                        gCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGccc                      >  1:134539/1‑80 (MQ=255)                                                        gCAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGccc                      <  1:134497/80‑1 (MQ=255)                                                         cAGACAAAGGACAAAGCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCa                     <  1:134498/80‑1 (MQ=255)                                                             caaaggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGa                >  2:134543/1‑81 (MQ=255)                                                              aaaggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGAc               <  2:134499/81‑1 (MQ=255)                                                               aaggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACg              >  1:134534/1‑81 (MQ=255)                                                                 ggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGgc            >  1:134530/1‑81 (MQ=255)                                                                 ggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGgc            <  1:134500/81‑1 (MQ=255)                                                                 ggacaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGgc            >  2:134544/1‑81 (MQ=255)                                                                   acaaagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGgcgc          <  2:134501/81‑1 (MQ=255)                                                                      aagCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAat       <  2:134502/81‑1 (MQ=255)                                                                       agCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAata      >  1:134536/1‑81 (MQ=255)                                                                        gCAATGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAatat     <  1:134503/81‑1 (MQ=255)                                                                          aaTGGTAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAa‑attat  <  2:134504/81‑1 (MQ=255)                                                                               |                                                                             GATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAACGTTGCAGACAAAGGACAAAGCAATGGCAATCCACAATCGTGCAGGCCAACCTGCACAACAGAGTGATTTGATTAACGTCGCCCAACTGACGGCGCAATATTAT  >  NC_000913/713484‑713638

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