Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,781,351	A→G	intergenic (+12/‑44)	ydiR → / → ydiS	putative electron transfer flavoprotein, FAD‑binding subunit/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,781,351	0	A	G	100.0%	211.2 / NA	57	intergenic (+12/‑44)	ydiR/ydiS	putative electron transfer flavoprotein, FAD‑binding subunit/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/31);  total (26/31)

TTGGCGGCAACGGCGCAAAAGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAACCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCTGGATGTCCTGGTGAT  >  NC_000913/1781203‑1781495                                                                                                                                                     |                                                                                                                                                 ttGACGGCAACGGCGCAAAAGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGcc                                                                                                                                                >  2:798900/1‑151 (MQ=255)        cAACGGCGCAAAAGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACAtctc                                                                                                                                         >  2:391199/1‑151 (MQ=255)         aaCGGCGCAAAAGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTcc                                                                                                                                        >  2:621435/1‑151 (MQ=255)                 aaaaGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTaa                                                                                                                                <  1:151130/151‑1 (MQ=255)                  aaaGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTaa                                                                                                                                >  2:163278/1‑150 (MQ=255)                   aaGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAAc                                                                                                                              <  2:563991/151‑1 (MQ=255)                       gATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGTCTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAg                                                                                                                          <  1:73477/151‑1 (MQ=255)                       gATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCa                                                                                                                           <  1:645862/150‑1 (MQ=255)                        aTTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGc                                                                                                                         >  1:345639/1‑151 (MQ=255)                        aTTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGc                                                                                                                         >  1:775839/1‑151 (MQ=255)                                    aaaaaaGATAAAAAGGCGCCAATCTTCAACAAGGCCGACTACGGTCTGGGGGGGGATATCACAAAAGCCGCCCCTGCCCTGA‑TAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAg                                                                                                               <  2:273195/145‑1 (MQ=255)                                    aaTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCt                                                                                                                                    <  2:150823/128‑1 (MQ=255)                                    aaTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCt                                                                                                                                    >  1:150800/1‑128 (MQ=255)                                    aaTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGc                                                                                                              >  1:565290/1‑150 (MQ=255)                                    aaTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCa                                                                                                             <  2:582366/151‑1 (MQ=255)                                     aTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCAt                                                                                                            >  2:830157/1‑151 (MQ=255)                                     aTAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCAt                                                                                                            <  1:638215/151‑1 (MQ=255)                                        aaGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAAc                                                                                                         >  1:244649/1‑151 (MQ=255)                                         aGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACg                                                                                                        >  1:213205/1‑151 (MQ=255)                                           aTAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGc                                                                                                       >  1:704027/1‑150 (MQ=255)                                               aaaTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTc                                                                                                  >  1:552984/1‑151 (MQ=255)                                               aaaTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTc                                                                                                  >  1:561169/1‑151 (MQ=255)                                               aaaTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTc                                                                                                  <  1:482927/151‑1 (MQ=255)                                                 aTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCgg                                                                                                <  2:123933/151‑1 (MQ=255)                                                 aTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCgg                                                                                                >  2:666353/1‑151 (MQ=255)                                                 aTGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCgg                                                                                                >  2:549810/1‑151 (MQ=255)                                                   gcgcCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGAt                                                                                              <  1:203912/151‑1 (MQ=255)                                                    cgcCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATg                                                                                             >  1:65825/1‑151 (MQ=255)                                                        aaTCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACaa                                                                                         <  2:268125/151‑1 (MQ=255)                                                         aTCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACaa                                                                                         >  1:520575/1‑150 (MQ=255)                                                           cTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAAtt                                                                                      >  2:684181/1‑151 (MQ=255)                                                             tCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAAttt                                                                                     <  2:520613/150‑1 (MQ=255)                                                             tCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTg                                                                                    >  1:635292/1‑151 (MQ=255)                                                               aaCTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGAt                                                                                  <  2:694410/151‑1 (MQ=255)                                                                   aTGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGc                                                                                <  2:820236/149‑1 (MQ=255)                                                                   aTGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGc                                                                                <  1:546427/149‑1 (MQ=255)                                                                      ccGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATtg                                                                           <  1:156804/151‑1 (MQ=255)                                                                       cGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATtgt                                                                          <  2:609555/151‑1 (MQ=255)                                                                            aCGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGt                                                                       <  1:54685/149‑1 (MQ=255)                                                                                         cGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTg                                                         <  1:265702/150‑1 (MQ=255)                                                                                                 aCAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGc                                                 <  2:418227/150‑1 (MQ=255)                                                                                                  cAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGt                                               <  1:136299/151‑1 (MQ=255)                                                                                                      gtcgtcCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGcc                                           >  1:859008/1‑151 (MQ=255)                                                                                                        cgtcCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGc                                         <  2:561186/151‑1 (MQ=255)                                                                                                                ccctgATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCa                                 <  2:400918/151‑1 (MQ=255)                                                                                                                  ctgATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATg                               >  2:763445/1‑151 (MQ=255)                                                                                                                        aGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACg                          <  1:162604/150‑1 (MQ=255)                                                                                                                        aGCCAGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACg                          >  1:856297/1‑150 (MQ=255)                                                                                                                           cAGTTGAGCCGCTAATTCCCCTTTCAGCCCCATCGCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAg                        <  1:462185/149‑1 (MQ=255)                                                                                                                            aGTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGc                       <  1:185516/149‑1 (MQ=255)                                                                                                                             gTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCgg                    <  2:690965/151‑1 (MQ=255)                                                                                                                             gTTGAGCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCgg                    <  2:380022/151‑1 (MQ=255)                                                                                                                                  gCCGCTAATTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCTg                <  1:466121/150‑1 (MQ=255)                                                                                                                                        aaTTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCt                 <  2:134835/143‑1 (MQ=255)                                                                                                                                        aaTTCCCCTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCt                 >  1:134812/1‑143 (MQ=255)                                                                                                                                               cTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCTGGATGTCCTGGTGAt  >  1:100084/1‑151 (MQ=255)                                                                                                                                               cTTTCAGCCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCTGGATGTCCTGGTGAt  >  1:203544/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTGGCGGCAACGGCGCAAAAGTGATTGTCGCCATCAATAAAGATAAAAATGCGCCAATCTTCAACTATGCCGACTACGGTCTGGTGGGCGATATCTACAAAGTCGTCCCTGCCCTGATTAGCCAGTTGAGCCGCTAATTCCCCTTTCAACCACATCTCCCGCTGTAAACGCAGCGGGAAGGGAGCATAACGCATGTCGGATGACAAATTTGATGCCATTGTGGTCGGTGCGGGCGTTGCTGGTAGCGTTGCCGCACTGGTCATGGCACGAGCCGGGCTGGATGTCCTGGTGAT  >  NC_000913/1781203‑1781495

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