Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,529,094	C→A	A361A (GCG→GCT)	ligA ←	DNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,529,094	0	C	A	100.0%	109.3 / NA	35	A361A (GCG→GCT)	ligA	DNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/14);  total (21/14)

TTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGCCGTCGATATCAAAGCCCAGC  >  NC_000913/2528854‑2529340                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        tgtCAGCGGCTCCTTACGCTGCGCACCGAAAAGCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACAGACCGGACAATGCGTCAGGAATACAACCGCACGGGTATCTTCCGGGCGGTCAGAAAGCACGACGTTAACCCCCTGCGGGATCACGGCGCCAGTGCAGCGAATCACCGCTTTATCGCCAAGTCGTAAACCAAGACGTTCGAGTTCATCAGCATTgtgt                                                                                                                                                                                                                                                <  2:154039/248‑1 (MQ=255)        gACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACACTGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGc                                                                                                                                                                                                                                        >  1:198521/1‑251 (MQ=255)         aCTCTTTACGCTGCGCACCGCAAATCAGGCCGCCGTTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCg                                                                                                                                                                                                                                       >  1:51989/1‑251 (MQ=255)            ctTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGtt                                                                                                                                                                                                                                     >  1:103860/1‑250 (MQ=255)              ttACGCTGCGCCCCCCAAATCAGGCCACCTGTACAGCGTGCAACCGCGTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACACCCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCCCCCCTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTAc                                                                                                                                                                                                                                   <  2:67611/250‑1 (MQ=255)              ttACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCCGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACt                                                                                                                                                                                                                                  >  1:267682/1‑251 (MQ=255)               tacgccgagcACAGCAAATCAGGCCACCGGTCCAGATCGAAACCGCTTCACCTACCACAAGCTCAACAACAGAACCACATACCGGACAATGCTACGGGAATACAACCTCACGGGTAACTTCCGGGCGCCCAGACAGCACGGCGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCAAGACACCACATTAACGCCCATGCGTAAAACATGACGTACGATTTCAACAGCATTGTGTAAGGTTGCGATACTc                                                                                                                                                                                                                                 >  2:19010/9‑251 (MQ=255)               tACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAAAGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACTGGCATCTTCAGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTCTCGCCACTGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTc                                                                                                                                                                                                                                 >  1:54500/1‑251 (MQ=255)                           cGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCt                                                                                                                                                                                                                     <  1:293042/251‑1 (MQ=255)                                aTCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAAc                                                                                                                                                                                                                >  1:129485/1‑251 (MQ=255)                                   aGGCCACCGGGACACCGGACATCCGCTTCACCTTCCACACGCGCAACTTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGGTCAGAAAGCACGACGGTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCGGCAACATg                                                                                                                                                                                                             <  2:41039/251‑1 (MQ=255)                                           ggTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCCGCGGGATCACTTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGGTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                                                                                                                                                                                     <  2:257902/251‑1 (MQ=255)                                                            cTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCCCGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACCCGCAACAggg                                                                                                                                                                                    >  2:10474/1‑250 (MQ=255)                                                                aCCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAAGCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACGCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTaa                                                                                                                                                                                <  2:32766/251‑1 (MQ=255)                                                                     ccACACGCTCAACGCCAGAACCACATACCGTACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCTCGTCGTTATCCACCTACGGGATCACGTCGCCATCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGTTTCCATCAGCATTGTGTAAGGTTTCGTTACTCACCAGCACGCCTGCATCATGGACAGGTTCCAGACGCGCAACAGGCGTAATcgc                                                                                                                                                                            <  1:270567/250‑1 (MQ=255)                                                                               aaCGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGcc                                                                                                                                                                 <  1:215250/251‑1 (MQ=255)                                                                                            caTACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTc                                                                                                                                                    <  1:150586/251‑1 (MQ=255)                                                                                                ccGGACAAGGCGTCGGGAATCCACCCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGACCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACg                                                                                                                                                <  2:67535/251‑1 (MQ=255)                                                                                                      aaTGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGc                                                                                                                                          >  1:68183/1‑251 (MQ=255)                                                                                                             cGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGCTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCCCGCACAAAgg                                                                                                                                   >  2:63255/1‑251 (MQ=255)                                                                                                                     cAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGc                                                                                                                           >  1:118663/1‑251 (MQ=255)                                                                                                                      aaCCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCt                                                                                                                          >  2:78442/1‑251 (MQ=255)                                                                                                                      aaCCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCt                                                                                                                          >  2:185864/1‑251 (MQ=255)                                                                                                                                                    aaGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCAcc                                                                                                                                                                                                                              <  2:211946/121‑1 (MQ=255)                                                                                                                                                    aaGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCAcc                                                                                                                                                                                                                              >  1:211946/1‑121 (MQ=255)                                                                                                                                                       cacgccatcaagcaccTGCGGGACCACGCCGTCAGCGCGGCGAATAACCACTTAACCTCCCATGCGTCAACAAAGACGATCGATTTCATCAGCATTGTGAAAGGCTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGATGCACACCAGGCGTACTCGCCCCAGTATGCACAACCTGAAACTCGCCTACACGCACAACGGACATCCGCTCCGGAGCCAGGAATTTAA‑CGTCTACCGCCCAgc                                                                                        >  2:264238/13‑251 (MQ=255)                                                                                                                                                                       ggCGGGATCCCGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGAc                                                                        <  2:74316/250‑1 (MQ=255)                                                                                                                                                                                           gcggcgAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCCCGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCctgct                                                    >  2:240421/1‑251 (MQ=255)                                                                                                                                                                                             ggcgAATCCCCACTCTATCGCCAATGCTTACACCAAGACGGTCGATTTCATCAGCATTGTGTAAGGTTGCGTTCCTCACCAGCACGCCTGAAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCATCCTGAAACACGACGCCACGCACAAAGGTCATCTTCTCCTGCGCAGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCtgt                                                  <  2:263368/251‑1 (MQ=255)                                                                                                                                                                                                                gCCAATGCGTAAACCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACATGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTa                                >  2:235357/1‑250 (MQ=255)                                                                                                                                                                                                                     tGCGTAAACCAAGACGTTCGTTTTCATCAGCCTTGTGTACGGTTGAGTTACTAACCCGCACGCCTGCAACATAGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGCCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCGCCTGATGGGCCCGTGGGTTGACCTTACTcacc                          >  2:303965/1‑251 (MQ=255)                                                                                                                                                                                                                          aaaCCAAGACGTTCGATTTCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGcc                     >  1:4254/1‑251 (MQ=255)                                                                                                                                                                                                                                aGACGTTCGATTTCATCAGCATTGTGTAAGGTTACGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGCCGTCGat               <  2:231255/251‑1 (MQ=255)                                                                                                                                                                                                                                           ttCATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGACTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGCCGTCGATATCACAGccc     >  1:166143/1‑250 (MQ=255)                                                                                                                                                                                                                                             cATCAGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGCCGTCGATATCAAAGCCCAGc  >  1:103518/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        TTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTACGCCCAACCTGAAACTCGACGTCACGCACAAAGGTCATCTGCTCCTGCGCCGGGAATTTAAACG‑CTACCGCCCAGCGCGGGGCACGCGCGACAAAGCCAAGCTGCTCCTGCTGTGCCAGTGAGTTGACCTTAATCACCACGCCGTCGATATCAAAGCCCAGC  >  NC_000913/2528854‑2529340

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