Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,770	A→G	H32H (CAT→CAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,770	0	A	G	100.0%	99.5 / NA	28	H32H (CAT→CAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

AAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACT  >  NC_000913/1310538‑1310999                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      aaaaCCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGAATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGt                                                                                                                                                                                                                     >  1:105572/1‑251 (MQ=255)  aaaCCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGGAACTGTAAACGTg                                                                                                                                                                                                                    >  2:90563/1‑251 (MQ=255)         ccctatAAATATTGAATCTCTATTCAAAACCCTTTCTTAAATGGTTTCCCTGAAACGTGTTCATAGCCTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAtg                                                                                                                                                                                                             <  2:323875/246‑1 (MQ=255)               aaaTATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATGCCGCAGCCCCTGGATCGTTACTGTCTACCTGTGGCATTGGAGCAGCTGTCAACAAACGACCGTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCg                                                                                                                                                                                                        <  1:150335/250‑1 (MQ=255)               aaaTATTGAAACTCTATTCAGAAAACTTTCTAAAATGGTTTCACTGAAACGTGTTCATAGACTACTGACGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCACCAAACCACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAtaggccgg                                                                                                                                                                                                       >  2:262458/1‑244 (MQ=255)                                  aGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCAGCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGtttt                                                                                                                                                                                    <  2:397918/251‑1 (MQ=255)                                              taagtGGTTTCCCTGAAACTTGTTCCTAGCCGCCTGCCGCTACGTACGGGTCAGCATCGTCCCAGGCCTGAGCTTCTTCCAGCGAGTCAAATTCAGCAATAACGGTTTGGCCAGTAAATCCCGCTGCCCCTGGATCGTTACTGTCTTCCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACTGAAAGGCGTTTTTCGAGGCTATc                                                                                                                                                                         <  2:250123/246‑1 (MQ=255)                                                     tttCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGt                                                                                                                                                                 >  1:224634/1‑251 (MQ=255)                                                           tGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTg                                                                                                                                                             >  1:229950/1‑249 (MQ=255)                                                            gAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGGCTTGa                                                                                                                                                            <  1:90563/249‑1 (MQ=255)                                                                  tgtTCATAGACTCCTGCCGCTAAGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGCTTAAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGCCCTTAATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACAGAAAGGCGATTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAAt                                                                                                                                                    >  2:332467/1‑251 (MQ=255)                                                                            cTCCTGTCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATacaac                                                                                                                                          >  1:122378/1‑251 (MQ=255)                                                                             tCCTGCCGCTACGTCCGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCACCAAAAGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGATTTTTCTTGAGCGTAAATAACATacaac                                                                                                                                          >  2:420859/1‑250 (MQ=255)                                                                                   cGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGATGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCACGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGc                                                                                                                                    >  1:381527/1‑250 (MQ=255)                                                                                                                            ccAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGa                                                                                           >  2:132234/1‑250 (MQ=255)                                                                                                                              aGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAg                                                                                        <  1:327112/251‑1 (MQ=255)                                                                                                                                cGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAggg                                                                                      >  1:149838/1‑251 (MQ=255)                                                                                                                                        aTTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATc                                                                              <  2:146630/251‑1 (MQ=255)                                                                                                                                                  aaCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGa                                                                     <  1:132234/250‑1 (MQ=255)                                                                                                                                                           gCCAGTAAATCCCGCAGCCCCTAGATCGTTACTGTCTACAGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGACGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGATGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGt                                                           >  1:343056/1‑251 (MQ=255)                                                                                                                                                                   aTCCCGCAGCCCCTGGATCGTTACTGTCTACCGATGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGtt                                                    >  1:197893/1‑250 (MQ=255)                                                                                                                                                                          aGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCAt                                            <  2:149838/251‑1 (MQ=255)                                                                                                                                                                               caggATCGTTACTGTATACCGCTGGCATTGGACCTGATGTCAACAAAAGACCTTCCTCGCGCAGTAACTGCTAACGCGCTAAATGCGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGCGCGTAAATAACATACAACACGGGCACAACTCATTGTTCTGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGcct                                        >  2:247802/3‑248 (MQ=255)                                                                                                                                                                                         aCTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGt                             <  1:262458/251‑1 (MQ=255)                                                                                                                                                                                                         aTTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCATCTTTGTCTTGAGCGTAAATAACAAACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTGATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTGTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATg             <  2:69684/251‑1 (MQ=255)                                                                                                                                                                                                            ggACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGaa           >  1:98180/1‑249 (MQ=255)                                                                                                                                                                                                              aCCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACAtt        <  1:90820/251‑1 (MQ=255)                                                                                                                                                                                                                     gTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATTTGCCGGTCGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCAct  <  2:313919/250‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      AAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACT  >  NC_000913/1310538‑1310999

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