Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,409	G→A	V29V (GTG→GTA)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,409	0	G	A	100.0%	105.3 / NA	36	V29V (GTG→GTA)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/7);  total (29/7)

ATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCC  >  NC_000913/4617187‑4617638                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         aTCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGGGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGAATGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGAACACCCTGAATGACGCCGATACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAAccc                                                                                                                                                                                                             >  1:99609/1‑249 (MQ=255)                        tGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CGGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACAGCCGCTATc                                                                                                                                                                                     >  2:43988/1‑251 (MQ=255)                                         aaTTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCtt                                                                                                                                                                    >  2:215748/1‑251 (MQ=255)                                             gtgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGCCGACACCGACGAGAAAGAAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTAACTGTATCTATCCTCGCTTTATc                                                                                                                                                                >  2:204722/1‑251 (MQ=255)                                             gtgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATc                                                                                                                                                                >  2:243380/1‑251 (MQ=255)                                                 tgtatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                             >  2:481760/5‑250 (MQ=255)                                                   tatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAtt                                                                                                                                                          >  1:171144/3‑251 (MQ=255)                                                          acgtgagtTGCAGATTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATAGGCGACAAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGAAGCCTGCGAGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCATGCCAAAACCCCGGTC‑GGCAATACCGCCTCTATCTGTCTCTATCCTCGCTTTATCCCTATTGCTCGc                                                                                                                                                    >  2:349449/7‑250 (MQ=255)                                                                       aGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGaa                                                                                                                                        >  1:432392/1‑248 (MQ=255)                                                                        gTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACAACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAag                                                                                                                                      >  2:259405/1‑249 (MQ=255)                                                                               ttAGAATACTAACAAACGCGCAAGGTGAATTTTCTTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGACTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAACACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑ACCACTGCAAGAGCAggac                                                                                                                              >  2:403661/1‑248 (MQ=255)                                                                                  aaaTACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGATAATGAAATGACTGCTCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATAACGACAACACCGACGAAAAAGAAATTGCT‑CTGGGTCATCAGGCCAAACCCCCGGTC‑GGAAAAACCGCCGCTCTCTGTATCTCTCCTCGCATTATCCCGCTTGCTCGCAA‑AAcaatcaaaaaccaggccaccc                                                                                                                          >  2:267087/2‑231 (MQ=255)                                                                                             gccTCGCAAGGGGTATTTTATTGGCGACAAACCAGGAGAATGAAAAGACTGATCTGAAAGCAAGAAGCCTGCGAGCACTGAAATTGATGGACC‑TGACCACCCTGAATGGCGACGACACCGACGCAAAAGTAATTGCT‑ATGATTCAACAGGCCAAAACCCCGGCC‑GGCAATACCGAACCTATCTGTATCTATCCTCTCTTAATCCCGATTGCTAGCAACGAC‑CTGAAAGAGCAGAGCCCCCCGGAAATcagt                                                                                                               >  2:418260/3‑248 (MQ=255)                                                                                                    aaGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCt                                                                                                         >  1:113940/1‑250 (MQ=255)                                                                                                     aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGCCTGATCTGAAAGCACGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                                       >  2:398917/1‑251 (MQ=255)                                                                                                           aTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                 >  1:466497/1‑251 (MQ=255)                                                                                                             tttATTGGCGACAAGCAAGGAGAATGAAATGACTGATCTGCAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAAAGAAGGGCACCTCGGAAATCAGTATTGCTACGGTaccc                                                                                                >  2:322958/1‑247 (MQ=255)                                                                                                                aggggcgacaggcagggAGAATGAAATGATGGATCGAAAAGCAAGCAGCCTGTGTGCACTAAAATTGTTGGACC‑TGACCACCCTGAATGACGACGACACCGACGAAAAAGTAATTGCT‑CGTTTTCATCAGGCAAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAAc                                                                                             <  1:308952/235‑1 (MQ=255)                                                                                                                               ggAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACCCGGTAacg                                                                             >  2:178662/1‑251 (MQ=255)                                                                                                                                gagaAAGAAATGACTGATCGAAAAGCAAGCAGCCTGCGTGCACTGAAATGGATGGACC‑TGACCCCCCTGAATGACGACGCCACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCAAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTTCGGTAACCAACTTCCCACACGGTAacga                                                                            <  1:459037/251‑1 (MQ=255)                                                                                                                                    agaaaatgatggATCAGAAAGCAGGCAGCATGCGTACACAAAAATTAATGGCCC‑GGCCCCCCCGGAAGGACGACCACCCCGACGAAAAATTATTTGCT‑CTGTGTCAACAGGCAAAACCCCCGGCC‑GTAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                        <  2:466497/240‑1 (MQ=255)                                                                                                                                        aatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGCCGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGAAACCAACTTCCCACACGGTAACGAtgacatcg                                                                    >  1:472046/1‑251 (MQ=255)                                                                                                                                             atGATCTAAAACCAAGCAGCCTGTGTGCACTGAAATGGA‑GGACCCTCACCCCCCTGAATGACGACGACCCCGCCGAAAAAGTAATTGCT‑CTGTTTCAGCAGGCCAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                               <  1:285173/250‑1 (MQ=255)                                                                                                                                                    gAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGCCC‑TGCCCCCCCGAAATGACGCCAACCCCGACGAAAAAGTAATTGCT‑CTGTTTCACCAGACAAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATcgc                                                             <  1:259405/246‑1 (MQ=255)                                                                                                                                                                 cTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                                           >  1:139142/1‑251 (MQ=255)                                                                                                                                                                   gCGTGCACTAAAATTTATGGACC‑TGCCCCCCCTGAATGCCGACGACCCCGACGATAAATTAATTTCT‑CTGTGTCTTTGGGCCAAACCCCCGGTC‑GCCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACTATGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                                         <  1:396346/251‑1 (MQ=255)                                                                                                                                                                     gTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CCGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                        >  2:303330/1‑250 (MQ=255)                                                                                                                                                                      tGCACTGAGATAGCTGGACC‑TGACCACCCTGAATTCCGACGACACCGACGAGAACGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAA‑AACCCTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGCCAGAAACCCGTGCGGcc                                       >  2:267937/1‑249 (MQ=255)                                                                                                                                                                                 tGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACggg                           >  1:252660/1‑250 (MQ=255)                                                                                                                                                                                                 cccTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTaa                                                                                                  >  1:346372/1‑165 (MQ=255)                                                                                                                                                                                                 cccTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTACGGCAATCGCCTACGGTGCCGATGAAGTTGAc            >  1:341694/1‑251 (MQ=255)                                                                                                                                                                                                   cTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAC‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTCCGGTGCCGATGAAGTTGACGt          >  1:301789/1‑251 (MQ=255)                                                                                                                                                                                                       aggacgacgacACCGACGAAAAAGTAAT‑GCCTCTGTGTCATC‑GGGCAAAACCCCCGTCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg      <  1:463292/249‑1 (MQ=255)                                                                                                                                                                                                       aTGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAAccc                                               >  2:6805/1‑210 (MQ=255)                                                                                                                                                                                                           cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc  >  1:236254/1‑251 (MQ=255)                                                                                                                                                                                                           cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc  >  1:121559/1‑251 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         ATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCC  >  NC_000913/4617187‑4617638

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: ‑