Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,586,834	A→G	intergenic (‑73/‑115)	hsdR ← / → mrr	type I restriction enzyme EcoKI endonuclease component/methylated adenine and cytosine restriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,586,834	0	A	G	100.0%	225.0 / NA	62	intergenic (‑73/‑115)	hsdR/mrr	type I restriction enzyme EcoKI endonuclease component/methylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/34);  total (28/34)

ATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTCCTACCTATGACA  >  NC_000913/4586703‑4586970                                                                                                                                    |                                                                                                                                         aTGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCg                                                                                                                                  <  1:130059/140‑1 (MQ=255)   ggCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCggg                                                                                                                                >  2:942260/1‑140 (MQ=255)        aaGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGcc                                                                                                                                        >  2:180531/1‑127 (MQ=255)        aaGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATg                                                                                                                           <  2:414703/140‑1 (MQ=255)           tGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGAcc                                                                                                                        <  2:361586/140‑1 (MQ=255)              aGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACg                                                                                                                     <  2:540858/140‑1 (MQ=255)              aGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACg                                                                                                                     >  2:11947/1‑140 (MQ=255)               gTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACgg                                                                                                                    >  2:829449/1‑140 (MQ=255)                   ttGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTcac                                                                                                                >  1:900666/1‑140 (MQ=255)                   ttGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTcac                                                                                                                <  2:105383/140‑1 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAAc                                                                                                          <  1:942260/140‑1 (MQ=255)                            ttCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGa                                                                                                       <  2:833501/140‑1 (MQ=255)                              cAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACt                                                                                                     <  1:2421/140‑1 (MQ=255)                                ggAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCa                                                                                                   <  1:85175/140‑1 (MQ=255)                                      cAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAAc                                                                                                          <  2:1067054/127‑1 (MQ=255)                                      cAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAAc                                                                                                          >  1:1067054/1‑127 (MQ=255)                                      cAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTa                                                                                             >  1:41990/1‑140 (MQ=255)                                      cAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTa                                                                                             >  1:726711/1‑140 (MQ=255)                                         aaTTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACg                                                                                          <  2:530879/140‑1 (MQ=255)                                              gATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTc                                                                                     >  2:665615/1‑140 (MQ=255)                                              gATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTc                                                                                     <  1:1452/140‑1 (MQ=255)                                                    ttCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTAt                                                                               <  1:11947/140‑1 (MQ=255)                                                    ttCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTAt                                                                               <  1:946692/140‑1 (MQ=255)                                                        tcatTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTAttttt                                                                           >  1:76592/1‑140 (MQ=255)                                                         catTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTg                                                                          <  1:747741/140‑1 (MQ=255)                                                         catTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTg                                                                          <  1:298296/140‑1 (MQ=255)                                                          atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGt                                                                         <  2:726711/140‑1 (MQ=255)                                                          atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAACAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGt                                                                         <  2:980858/140‑1 (MQ=255)                                                           ttGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTa                                                                        <  2:638655/140‑1 (MQ=255)                                                                      ccATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTAtagta                                                             >  2:549809/1‑140 (MQ=255)                                                                       cATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTAtagtag                                                            <  1:927642/140‑1 (MQ=255)                                                                        aTTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGc                                                           >  2:902757/1‑140 (MQ=255)                                                                        aTTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGc                                                           >  1:488514/1‑140 (MQ=255)                                                                        aTTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATATTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGc                                                           >  1:148237/1‑140 (MQ=255)                                                                          tGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTg                                                                          <  1:928313/123‑1 (MQ=255)                                                                          tGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTg                                                                          >  2:928313/1‑123 (MQ=255)                                                                           gCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGa                                                        >  2:492183/1‑140 (MQ=255)                                                                               tgCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGaaaaa                                                    >  2:1062362/1‑140 (MQ=255)                                                                                 cGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACa                                                  >  1:395554/1‑140 (MQ=255)                                                                                 cGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACa                                                  >  1:413513/1‑140 (MQ=255)                                                                                  gggCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACAt                                                 <  1:492183/140‑1 (MQ=255)                                                                                   ggCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATc                                                >  2:752831/1‑140 (MQ=255)                                                                                         tCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTg                                          >  1:574221/1‑140 (MQ=255)                                                                                         tCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTg                                          >  1:455342/1‑140 (MQ=255)                                                                                            aaaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGAtt                                       <  1:665615/140‑1 (MQ=255)                                                                                                atTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGc                                   <  2:76592/140‑1 (MQ=255)                                                                                                 tttAAGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCa                                  <  1:902757/140‑1 (MQ=255)                                                                                                    aaGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               <  2:574221/140‑1 (MQ=255)                                                                                                    aaGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               <  2:488514/140‑1 (MQ=255)                                                                                                    aaGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               <  2:413513/140‑1 (MQ=255)                                                                                                    aaGGCCCATAACATCTCATATTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               <  2:148237/140‑1 (MQ=255)                                                                                                     aGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGa                              <  2:41990/140‑1 (MQ=255)                                                                                                     aGGCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGa                              <  2:1066520/140‑1 (MQ=255)                                                                                                      ggCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               >  2:327940/1‑138 (MQ=255)                                                                                                      ggCCCATAACATCTCATCTTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAgg                               <  1:327940/138‑1 (MQ=255)                                                                                                               catctcatctTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGa                    <  2:286539/140‑1 (MQ=255)                                                                                                                  ctcatctTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACgg                 >  2:377705/1‑140 (MQ=255)                                                                                                                    catctTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGtt               >  1:790455/1‑140 (MQ=255)                                                                                                                       ctTAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTcct            >  2:516523/1‑140 (MQ=255)                                                                                                                         tAGCTTTCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTcctac          <  2:723327/140‑1 (MQ=255)                                                                                                                               tCTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTCCTACCTATGa    >  2:272686/1‑140 (MQ=255)                                                                                                                                cTGTGCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTCCTACCTATGAc   >  2:224713/1‑140 (MQ=255)                                                                                                                                 tgtgCCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTCCTACCTATGACa  >  1:1114596/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                         ATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGTTCCTACCTATGACA  >  NC_000913/4586703‑4586970

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