Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,065,765	A→G	intergenic (‑1/+540)	cobU ← / ← insH‑6	cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,065,765	0	A	G	100.0%	222.3 / NA	63	intergenic (‑1/+540)	cobU/insH‑6	cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/37);  total (26/37)

TTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATT  >  NC_000913/2065630‑2065899                                                                                                                                        |                                                                                                                                       ttCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGa                                                                                                                                    >  2:657926/1‑140 (MQ=255)      tCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCt                                                                                                                               >  2:547808/1‑140 (MQ=255)      tCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCt                                                                                                                               <  1:551272/140‑1 (MQ=255)         gTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGAt                                                                                                                            >  1:233425/1‑140 (MQ=255)          tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGAtt                                                                                                                           <  1:418482/140‑1 (MQ=255)          tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGAtt                                                                                                                           <  1:484280/140‑1 (MQ=255)          tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGAtt                                                                                                                           <  1:735206/140‑1 (MQ=255)           gCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTc                                                                                                                          >  1:770559/1‑140 (MQ=255)           gCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGCGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTc                                                                                                                          <  2:740955/140‑1 (MQ=255)             aGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCaa                                                                                                                        <  2:770559/140‑1 (MQ=255)             aGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCaa                                                                                                                        <  1:71953/140‑1 (MQ=255)             aGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCaa                                                                                                                        <  1:655501/140‑1 (MQ=255)               ccATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAAt                                                                                                                      <  2:223439/140‑1 (MQ=255)                cATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATg                                                                                                                     >  1:495082/1‑140 (MQ=255)                  tctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTa                                                                                                                   <  1:657926/140‑1 (MQ=255)                  tctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTa                                                                                                                   <  2:233425/140‑1 (MQ=255)                    tcatcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTatg                                                                                                                 <  1:574261/140‑1 (MQ=255)                    tcatcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCAGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTatg                                                                                                                 <  2:58375/140‑1 (MQ=255)                         atcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCTACTCCTGATTCAAAt                                                                                                                      <  1:178466/130‑1 (MQ=255)                         atcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCTACTCCTGATTCAAAt                                                                                                                      >  2:178466/1‑130 (MQ=255)                         atcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGt                                                                                                            <  2:208863/140‑1 (MQ=255)                         atcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGt                                                                                                            <  1:617683/140‑1 (MQ=255)                               gATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCaa                                                                                                                        <  2:806456/122‑1 (MQ=255)                               gATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCaa                                                                                                                        >  1:806456/1‑122 (MQ=255)                                aTTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCAtt                                                                                                     >  1:231640/1‑140 (MQ=255)                                 tttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTg                                                                                                    <  2:332523/140‑1 (MQ=255)                                  ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGt                                                                                                            >  1:627407/1‑131 (MQ=255)                                  ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGt                                                                                                            <  2:627407/131‑1 (MQ=255)                                   tGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                                                  <  2:282589/140‑1 (MQ=255)                                   tGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                                                  >  2:168505/1‑140 (MQ=255)                                            gCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAAt                                                                                         <  1:168505/140‑1 (MQ=255)                                             cGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATg                                                                                        <  2:172558/140‑1 (MQ=255)                                                    cAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATcccc                                                                                 >  1:324736/1‑140 (MQ=255)                                                    cAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATcccc                                                                                 >  1:562282/1‑140 (MQ=255)                                                       aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCAc                                                                              >  1:162032/1‑140 (MQ=255)                                                          cTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGa                                                                                      <  1:817800/129‑1 (MQ=255)                                                          cTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGa                                                                                      >  2:817800/1‑129 (MQ=255)                                                            gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATcc                                                                                   >  1:632163/1‑130 (MQ=255)                                                            gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATcc                                                                                   <  2:632163/130‑1 (MQ=255)                                                            gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTg                                                                         >  2:432579/1‑140 (MQ=255)                                                             tgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc                                                                        >  1:685962/1‑140 (MQ=255)                                                              gAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt                                                                       >  2:14457/1‑140 (MQ=255)                                                                 gagTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGAc                                                                    <  2:559787/140‑1 (MQ=255)                                                                    tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGc                                                                 >  1:425601/1‑140 (MQ=255)                                                                              gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTa                                                       >  2:86227/1‑140 (MQ=255)                                                                               cgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                                      >  1:446586/1‑140 (MQ=255)                                                                               cgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                                      <  2:231640/140‑1 (MQ=255)                                                                               cgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                                      <  2:495082/140‑1 (MQ=255)                                                                                         tcgCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATa                                            <  2:562282/138‑1 (MQ=255)                                                                                         tGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATa                                            <  2:324736/140‑1 (MQ=255)                                                                                         tGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATa                                            <  1:72948/140‑1 (MQ=255)                                                                                           gCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACt                                          <  1:86227/140‑1 (MQ=255)                                                                                                 ctTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGt                                    >  2:299083/1‑140 (MQ=255)                                                                                                  ttCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTa                                   <  2:685962/140‑1 (MQ=255)                                                                                                            gTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTaa                         >  1:99317/1‑140 (MQ=255)                                                                                                            gTGCGCCCCCCGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTaa                         <  1:299083/140‑1 (MQ=255)                                                                                                                     ccGTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGc                <  1:14457/140‑1 (MQ=255)                                                                                                                       gTCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCAt              >  1:270411/1‑140 (MQ=255)                                                                                                                        tCACCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCAtt             >  2:734228/1‑140 (MQ=255)                                                                                                                          aCCAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATtct           <  1:734228/140‑1 (MQ=255)                                                                                                                            cAAAATCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCa         >  2:459837/1‑140 (MQ=255)                                                                                                                               aaTCATCATGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATc      <  1:472869/140‑1 (MQ=255)                                                                                                                                   atcatGCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCAtt  <  2:99317/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       TTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATT  >  NC_000913/2065630‑2065899

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