Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,235,651	T→C	intergenic (+244/‑255)	pgi → / → yjbE	glucose‑6‑phosphate isomerase/uncharacterized protein YjbE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,235,651	0	T	C	100.0%	122.1 / NA	37	intergenic (+244/‑255)	pgi/yjbE	glucose‑6‑phosphate isomerase/uncharacterized protein YjbE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)

TAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAA  >  NC_000913/4235405‑4235890                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                caacatgtgtggATATGTAGGCCGGATAGGGCGTTCACGCCGCATCCCGCAATCGATGGCTGATTCGACGCGGTCGAGTATTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAAGGCCTCCGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTTTATTTTAACGGAAATCATACCGTGAGGTTAACCCt                                                                                                                                                                                                                                              <  2:408819/239‑1 (MQ=255)     catcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAAt                                                                                                                                                                                                                                         >  1:227788/1‑251 (MQ=255)         gtcgATATGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCAACAGGTCCATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGACGACGTCATGTCTGCATTTTTTTATCAGATAAATCCCCTTTTCTGTAATTTAAAGGTNATCCTACCTTGAGGTTAACCCTAAAAtcatt                                                                                                                                                                                                                                     >  2:266036/1‑247 (MQ=255)            gaaaggtatgCCGGATAAGGCGTTCTCGGCGCAGCCGTCAACCGATGCCTGATTCGACGCGGTCGCGTCTTATCAGGCATACAGGATGATGCCGATATGTACAGCGTATTCGGCAATTAATACATAGCTCGATTGATTAAATATCCTTAATAACAATGCCGACGTTATGTCGGCATTTTTGTATCATATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAAGAGAtttt                                                                                                                                                                                                                                  <  2:462052/242‑1 (MQ=255)             atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAttttt                                                                                                                                                                                                                                 <  1:464409/251‑1 (MQ=255)             atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTAATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAttttt                                                                                                                                                                                                                                 >  1:293104/1‑251 (MQ=255)                           tAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTANCCCTAAAATAGATTTTTAATCGTTGTTTAtt                                                                                                                                                                                                                   >  2:233990/1‑251 (MQ=255)                                  ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTCTCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCgttgggacgtgaaacc                                                                                                                                                                                                             >  1:140020/1‑238 (MQ=255)                                    cACGCCGCATCCGGCAACCGATGCCTGATGCNACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAt                                                                                                                                                                                                          <  2:253927/251‑1 (MQ=255)                                                      cGATGCCTGATGGTACGCAGTCGCGTCTTATCAGGCATACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTGTATTTCGGAAAATACGCCGATTAATTGCttt                                                                                                                                                                                        <  2:397374/251‑1 (MQ=255)                                                                      gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGt                                                                                                                                                                        >  1:139030/1‑251 (MQ=255)                                                                        ggTCGCGTCTTATCAGGCATACAGGTCGATGACGATATGTAAATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCAGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTATGTAATTTAACGGAAATAATCCCGTGAGGTTAACCCTAAAATAGTTTTTTAATCGTTGTTTATTTCGACAAATACGCAGATTAATTGATTTTGTTTTTATTTTAAGttt                                                                                                                                                                      >  2:494537/1‑251 (MQ=255)                                                                                  tagcaggcatacatgtcggtgtagATATGTACATCGTTTTCTGCAATTAATACATAGCACGATTGATTAAATAACCATAATAACAATGCCTACGTTATTTAGGCATTTTTTTATCAAATAAATCCCCTTGTCTGTTATTTAATGGAAATCATACTGTGAGGTTAACTCTAAAATAGATTTTTAATCGTTGTTTAATTCGGAAAAGACTCATATTAATTGTTTTTGATTTTATTTTAAGATTATGAtgttt                                                                                                                                                             <  2:326912/227‑1 (MQ=255)                                                                                       ggCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATGTTTAATCGTTGTGTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGGTTATGATTTTTATTGtt                                                                                                                                                       <  2:184389/251‑1 (MQ=255)                                                                                                  cGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                                                                                                                                             >  1:355687/1‑250 (MQ=255)                                                                                                  cGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                                                                                                                                             >  1:232092/1‑250 (MQ=255)                                                                                                       ccGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTg                                                                                                                                       <  1:401114/251‑1 (MQ=255)                                                                                                            atGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACt                                                                                                                                  >  1:447609/1‑251 (MQ=255)                                                                                                              gTACAGCGTTTTCGGCAATTAATACATAGAACGATTGATTAAATAACCTTAATAACAATGCCTACGTTATGTCTGCATTTTTTTATCAGATAAATACCCTTTTCTGTAATTTAACGGAAATCATACTGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTGTTATTGTTATTTAAATATAAGTTGAAACtt                                                                                                                                 <  2:425344/250‑1 (MQ=255)                                                                                                                    cGTATTCGGCAATTAATACATAGCACGATTNATTAAATAACCTTAATAACAATTCCGGCGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCATATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATAttt                                                                                                                           <  2:242311/250‑1 (MQ=255)                                                                                                                                      cATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTACATATAAGTTGCAACTTATATTTGATATTCATTCCAATTATc                                                                                                        >  2:380192/1‑251 (MQ=255)                                                                                                                                                   attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGCAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATc                                                                                           >  1:316833/1‑251 (MQ=255)                                                                                                                                                      aaaTAACCTTAAAAAGAATGCCGTCTATATGTCGGCATTTTTTTATCAGATACATCCCCTTGTCTGTACTTTAGCGGAAATAATACCGTGAGGTTACCCCTCAATTAGCTTTTTCATCGTTGTTTATTTCGGAAAATACGCAGATTACTTGTTTTTGTTTTTCTTTTAAGTTTATGATTTTTATTGTTATTTACCTATAAGTTGAAACTTATATTTGAAATTCATACACATTATCCTAAACCGCCATGGCt                                                                                        >  2:40670/1‑251 (MQ=255)                                                                                                                                                                     cAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGt                                                                         >  1:173546/1‑251 (MQ=255)                                                                                                                                                                                          aTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTccc                                                                                 >  2:219484/1‑222 (MQ=255)                                                                                                                                                                                          aTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTccc                                                                                 <  1:219484/222‑1 (MQ=255)                                                                                                                                                                                              ttctATCAGATAAATCCACTTGTCTGTAATTTAACGGAAATCATACCGTGAGTTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAAACGCAGATTAATTGCGTTTGTTTTTATTTTATGATTATGATTTTTTTTGTTATTTAAATATTAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCATTCGTTAATTTCCAGCTCCGTAATTCGTATGTTTTAGTTTTGtat                                                 >  2:471227/4‑250 (MQ=255)                                                                                                                                                                                                  aTCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTc                                             <  1:216833/250‑1 (MQ=255)                                                                                                                                                                                                          aaTCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATcccgcc                                    <  1:238013/251‑1 (MQ=255)                                                                                                                                                                                                          aaTCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATcccgcc                                    <  1:261332/251‑1 (MQ=255)                                                                                                                                                                                                                   gTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGtt                           >  1:200162/1‑251 (MQ=255)                                                                                                                                                                                                                                  gAAATCATACCGTGAGGTTAACCCTAAAATANATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTcc            <  2:266352/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCcat          <  1:93004/251‑1 (MQ=255)                                                                                                                                                                                                                                       cATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCcattc        <  2:16466/250‑1 (MQ=255)                                                                                                                                                                                                                                          aCCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCattcattc    >  2:139720/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCAGCGCTAATTCCCCGCGCCGTAATTCGCATGCTGTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCAttcattca   <  2:27864/251‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATTCGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCaa  <  1:380192/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                TAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAA  >  NC_000913/4235405‑4235890

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