Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,821,148	G→A	P23S (CCA→TCA)	chbC ←	N,N'‑diacetylchitobiose‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,821,148	0	G	A	100.0%	59.5 / NA	18	P23S (CCA→TCA)	chbC	N,N'‑diacetylchitobiose‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/8);  total (10/8)

AATAAAGAACGGTGCCATTAAAGACATTATTCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACTTC  >  NC_000913/1820909‑1821391                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    aaTAAAGAACGGTGCCATTAAAGACATTATTCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTa                                                                                                                                                                                                                                          <  2:196279/251‑1 (MQ=255)      aGAACGGTGCCATTAAAGAAATTATTCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTAtttt                                                                                                                                                                                                                                      >  2:11260/1‑250 (MQ=255)              gCCATTAAAGACATTATTCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCaa                                                                                                                                                                                                                              >  2:106687/1‑250 (MQ=255)                          attattCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTAc                                                                                                                                                                                                                 >  1:123922/1‑251 (MQ=255)                             attCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCt                                                                                                                                                                                                               <  1:11358/250‑1 (MQ=255)                                ccTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCtttt                                                                                                                                                                                                            <  1:106687/250‑1 (MQ=255)                                                                    ttCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATAcc                                                                                                                                                                        <  1:185742/250‑1 (MQ=255)                                                                                       cAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAGCCGCAATTTaa                                                                                                                                                     <  1:55107/250‑1 (MQ=255)                                                                                                          gCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTAt                                                                                                                                 >  1:14961/1‑251 (MQ=255)                                                                                                                    aaGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAAc                                                                                                                         <  1:52546/249‑1 (MQ=255)                                                                                                                                  cGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACAGTACTCATATAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAAt                                                                                                         <  1:122777/251‑1 (MQ=255)                                                                                                                                                                            aaaCATGGCCCCCGAAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCa                                                                >  1:190256/1‑250 (MQ=255)                                                                                                                                                                                      cccGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCtt                                                     >  1:137936/1‑251 (MQ=255)                                                                                                                                                                                            aGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATAcc                                                                                                                                                                        <  2:59693/130‑1 (MQ=255)                                                                                                                                                                                            aGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATAcc                                                                                                                                                                        >  1:59693/1‑130 (MQ=255)                                                                                                                                                                                                          tAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCg                                  >  1:3259/1‑250 (MQ=255)                                                                                                                                                                                                                                     aTTAACGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCACGAGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCATTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAAtt       >  1:37313/1‑250 (MQ=255)                                                                                                                                                                                                                                          cGTGTGACTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACTTc  >  2:55899/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    AATAAAGAACGGTGCCATTAAAGACATTATTCCTAATGTTCCGTTATATACGTTGCCGCCAATACCTTTCAGACCATTAAGTGTTTCAATGGTTGAGGCGTCGAGGCGAATACCTAAGGAATAAAAAAACGACCCCTCCCCAAAGCTTAGAAAAACGTTGTTAATTAATACAAACATGGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACTTC  >  NC_000913/1820909‑1821391

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