Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,964,330	C→A	intergenic (‑153/+31)	ppdA ← / ← thyA	conserved protein PpdA/thymidylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,964,330	0	C	A	100.0%	91.5 / NA	26	intergenic (‑153/+31)	ppdA/thyA	conserved protein PpdA/thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/14);  total (12/14)

AGGAGAAAAAAAGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTT  >  NC_000913/2964195‑2964459                                                                                                                                          |                                                                                                                                  aggagAAAAAAAGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGc                                                                                                                               >  2:485277/1‑140 (MQ=255)    agaAAAAAAGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTc                                                                                                                            >  1:467491/1‑140 (MQ=255)      aaaaaaaGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCgg                                                                                                                          <  1:317956/140‑1 (MQ=255)       aaaaaaGACGAGGGCCAGTTTCTATTTCTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCTGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCg                                                                                                                           >  2:398623/1‑140 (MQ=255)           aaGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTg                                                                                                                     <  2:467491/140‑1 (MQ=255)             gACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGc                                                                                                                   >  1:318820/1‑140 (MQ=255)                   ggCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGCTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCCCTGGCAGGATg                                                                                                             >  1:247253/1‑140 (MQ=255)                      cAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGttt                                                                                                          <  1:485277/140‑1 (MQ=255)                        gTTTCTATTTCTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCTGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGttt                                                                                                          <  1:398623/140‑1 (MQ=255)                                    cttcGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATa                                                                                              <  1:23023/140‑1 (MQ=255)                                           gcATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCg                                                                                       >  2:132926/1‑140 (MQ=255)                                              tCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGgcg                                                                                    <  1:53797/140‑1 (MQ=255)                                                        tATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGcc                                                                          <  1:283881/140‑1 (MQ=255)                                                        tATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGcc                                                                          <  2:290520/140‑1 (MQ=255)                                                          tttACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCAc                                                                                         <  2:352912/123‑1 (MQ=255)                                                          tttACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCAc                                                                                         >  1:352912/1‑123 (MQ=255)                                                            tACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGa                                                                      >  1:471583/1‑140 (MQ=255)                                                                    tGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGGAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATc                                                              >  1:313682/1‑140 (MQ=255)                                                                     gCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCg                                                             <  2:318820/140‑1 (MQ=255)                                                                        ggACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATACGGATCGTAg                                                          <  2:247253/140‑1 (MQ=255)                                                                           cGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGt                                                            <  2:174137/135‑1 (MQ=255)                                                                           cGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGt                                                            >  1:174137/1‑135 (MQ=255)                                                                           cGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCt                                                       >  1:267105/1‑140 (MQ=255)                                                                             ttGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTc                                                     <  2:419431/140‑1 (MQ=255)                                                                              tGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCa                                                    >  2:15039/1‑140 (MQ=255)                                                                                                                        tAAAAAAACCGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGt          <  1:15039/140‑1 (MQ=255)                                                                                                                                ccGACGCACAATGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGttt  >  1:29097/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  AGGAGAAAAAAAGACGAGGGCCAGTTTCTATT‑‑TCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTT  >  NC_000913/2964195‑2964459

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