Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,269,449	C→T	C12C (TGC→TGT)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,269,449	0	C	T	100.0%	81.5 / NA	28	C12C (TGC→TGT)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/11);  total (17/11)

AAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAG  >  NC_000913/4269202‑4269688                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 aaaTTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑a                                                                                                                                                                                                                                              <  1:222984/251‑1 (MQ=255)     tAATGTATCCTTACATCGAGTAATAAACATTTTTTATACCAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGt                                                                                                                                                                                                                                          >  1:509200/1‑251 (MQ=255)         gTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGc                                                                                                                                                                                                                                      >  1:250885/1‑251 (MQ=255)                  cATCGATTCATAAACAGTTTTTATACAAAAAAAGTCATGATCGTATTTACTGGGTTAAACATATTCATTCTGGTTTTAAATACAAAAAACGCAGATAATCCTACCATTCCTCGGCAGCAGGCAGTATGGACAGATGTTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACTTACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTTTAAACAg                                                                                                                                                                                                                             <  2:268853/251‑1 (MQ=255)                   aTCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTAGTTACTGGGTTAAATATAATCATCCTGCTTTTCAACACAAAAACCGCAGATAAGCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAAGCTCACAATAAAAAGGTTCAACATACTGACTATGTAGGGAAAAATATGCTCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGt                                                                                                                                                                                                                            <  1:539527/251‑1 (MQ=255)                                 aTTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGTTTAAATATAAACATCCTGCTTTTCATCTCAAAAACCGCAGCTAATCCTTCCAGTCCCCGACAGCTGGCGTTATGGTCAGATGGTTTGTGCAACAAATCTCACAATATAAAGGTTCAACATACTGCCTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAACCAGTGCCGTGTGTTT‑ATTGTTCGCTCTAAACCGTTCCGCTGTTGCCCt                                                                                                                                                                                                              <  2:499288/251‑1 (MQ=255)                                           cAAAAAAAGACAGGAACGTTTTTACTGGGTTAAATATAATCATCCTACTTTTCATCACAAAAACCGCAGATAATACTTCCTTTCCCCGGAAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAATTCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCAGTTTGTTT‑ATTGTTCGCTCTAAACAGATCCGCTGTTGCACTGACATCAtct                                                                                                                                                                                                    >  2:157635/1‑251 (MQ=255)                                           cAAAAAAAGACAGGAACGCATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATACTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAACTATGAGCAAGATCACACAGGCAATCAGTGCTGTTTGTTT‑CTGATTAGATCTAAACAGTTCCGCTGTTGCCCTGGCCTCAtct                                                                                                                                                                                                    >  2:339239/1‑251 (MQ=255)                                                                                aaTCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAgg                                                                                                                                                                <  2:114626/250‑1 (MQ=255)                                                                                  tcatcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCTCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑At                                                                                                                                                                                                                                             <  2:517510/171‑1 (MQ=255)                                                                                  tcatcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑At                                                                                                                                                                                                                                             >  1:517510/1‑171 (MQ=255)                                                                                                   ccaacacccGCAGATAATTCTTCCTTTCCCAGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGCACAGGCACCcat                                                                                                                                            >  2:51247/8‑251 (MQ=255)                                                                                                     aaaaaCCGCAGATAATCCTTCCTTTCCCCGTCAGCTGGCGTTATGGCCAGATGGTTTTTGCAACACATCTCACAATAAAAAGTTTCAACATACTGACTATTTCGGGAAAAATATGCGCAAGATCACACAGTAAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCCCCGCTAAACCCTGGTCCTAACGTTGCCCGGCTTGCTGAACAGGCCCCcata                                                                                                                                           >  2:12937/1‑249 (MQ=255)                                                                                                      cacaCCGCAGAGAATCCTTCCTTTCCCCGGCCGCTGGCGTTATGGACATATGGATTTTGCCACACCTCTCACATTAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGAGCCGTTTGTTT‑ATTGTTCGCTCTAAAACGTTCCGCTGTTGCCCTGGCATAATCTCCTTCACCGCTTACCCCTGGGACAAACGTTGCCAGGCCTGCTGACCAggaaagatcggc                                                                                                                                         >  2:205539/4‑241 (MQ=255)                                                                                                      aaaaCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGGTCCGCTGTTGCCCGGGCCTCATCGCCTTCTCCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCa                                                                                                                                                 <  1:205539/243‑1 (MQ=255)                                                                                                                                                         ggTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAAt                                                                                      >  1:142697/1‑251 (MQ=255)                                                                                                                                                             tttGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGc                                                                                   <  1:174177/250‑1 (MQ=255)                                                                                                                                                                   acaaATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTgggg                                                                            >  1:247644/1‑251 (MQ=255)                                                                                                                                                                                 aaaaaGTTCCAACAATCTGACTATTTAGGGAAAAATATGGTCAAGATCCCCCAGGCAATCAGTGCC‑TTTGTTTTATTGGTCGCTCTAATCAGATCCGGTCTGGCCCTGGCCTCACCTCCTCCGCCGCTTAACCCTTGGACTAACGTTGCCAGGCTTGCTGAACAGGCGCCCATTCATTGGGTTTCGGTCGCAGAAACTGAATATATCTACGCATGGCTTCCGCCAATGGCGGTGGGGTTTGATATCGATGa                                                              <  2:183490/251‑1 (MQ=255)                                                                                                                                                                                        ttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCCCAAATTGAAAATAGCCTCGCAGGGCGACCGCCAATGGCGGttgatcggagagccgtcgtgtagggg                                                       >  2:111071/1‑226 (MQ=255)                                                                                                                                                                                        ttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGt                                                                                <  1:111071/226‑1 (MQ=255)                                                                                                                                                                                                          tAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGGTATCGATGACCCGGTACTTTTTTCCAGTCCGGGc                                     >  1:94747/1‑251 (MQ=255)                                                                                                                                                                                                            gggAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCTGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCtt                                   <  1:358699/251‑1 (MQ=255)                                                                                                                                                                                                                aaaaTATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTTCAGTCCGGGCTTTTgg                               >  1:383407/1‑251 (MQ=255)                                                                                                                                                                                                                         gcAAGATCACACAGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGGTGACACGGTACTTTTTTTCAGTCCGGGCTTTTGGCGCGGCaa                       >  1:561721/1‑250 (MQ=255)                                                                                                                                                                                                                                   acaGGCAATCAGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCCAAAAAAACTtct            >  1:44670/1‑251 (MQ=255)                                                                                                                                                                                                                                             aGTGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCCGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAgca    >  1:462088/1‑247 (MQ=255)                                                                                                                                                                                                                                               tGCCGTTTGTTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGGACTTTTTTTCAGTCCGGGCTTCTGGCGCGGCCAAAAAACATTCTCGCCAGAAAg  >  1:210664/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 AAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTT‑ATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAG  >  NC_000913/4269202‑4269688

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