Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,331,149	C→A	pseudogene (3254/4605 nt)	yfaS ←	putative macroglobulin family protein, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,331,149	0	C	A	100.0%	214.3 / NA	60	pseudogene (3254/4605 nt)	yfaS	putative macroglobulin family protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/32);  total (28/32)

GCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAACGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAGCCGGATATTG  >  NC_000913/2331014‑2331280                                                                                                                                        |                                                                                                                                    gCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCa                                                                                                                                 >  1:175975/1‑140 (MQ=255)     cACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCg                                                                                                                             >  2:347210/1‑140 (MQ=255)             ccATCTTCCCCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAg                                                                                                                     <  2:135163/140‑1 (MQ=255)                   tCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGAt                                                                                                               <  2:530038/140‑1 (MQ=255)                          accaGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACgg                                                                                                        >  1:261503/1‑140 (MQ=255)                          accaGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGC‑GGATCAGACGGc                                                                                                       >  2:130962/1‑140 (MQ=255)                            caGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAg                                                                                                                     >  2:776597/1‑125 (MQ=255)                            caGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAg                                                                                                                     <  1:776597/125‑1 (MQ=255)                              ggTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTAc                                                                                                    <  1:616642/140‑1 (MQ=255)                              ggTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTAc                                                                                                    <  1:455504/140‑1 (MQ=255)                                tAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCg                                                                                                  <  1:243791/140‑1 (MQ=255)                                tAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCg                                                                                                  <  1:670754/140‑1 (MQ=255)                                   aGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGgtg                                                                                               >  1:849058/1‑140 (MQ=255)                                   aGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGgtg                                                                                               >  1:59583/1‑140 (MQ=255)                                   aGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGgtg                                                                                               >  1:469083/1‑140 (MQ=255)                                   aGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGgtg                                                                                               <  2:175975/140‑1 (MQ=255)                                     cgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGtt                                                                                             >  1:477199/1‑140 (MQ=255)                                       cgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGGTGATCGGCAAGCGAAAGCCAGGCGAGCCTGAGCGGGCTCAGAACGCCACCGGGGTTGa                                                                                           >  2:616642/1‑140 (MQ=255)                                       cgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGa                                                                                           >  2:960234/1‑140 (MQ=255)                                         cTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGAtt                                                                                         >  2:732338/1‑140 (MQ=255)                                         cTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCAACTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGAtt                                                                                         >  2:566801/1‑140 (MQ=255)                                          cccGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTa                                                                                        <  2:365572/139‑1 (MQ=255)                                             gggCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGc                                                                                     <  2:385006/140‑1 (MQ=255)                                               gCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACg                                                                                      >  1:1059423/1‑137 (MQ=255)                                               gCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACg                                                                                      <  2:1059423/137‑1 (MQ=255)                                                  ccGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCAccc                                                                                <  2:16732/140‑1 (MQ=255)                                                    gCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCa                                                                              >  1:819927/1‑140 (MQ=255)                                                    gCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCa                                                                              >  1:666918/1‑140 (MQ=255)                                                      cAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACg                                                                            <  1:276650/140‑1 (MQ=255)                                                            cATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCAt                                                                      <  1:960234/140‑1 (MQ=255)                                                               cAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGa                                                                   <  1:922545/140‑1 (MQ=255)                                                                  ccGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCa                                                                       <  1:794066/133‑1 (MQ=255)                                                                  ccGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCa                                                                       >  2:794066/1‑133 (MQ=255)                                                                                         tcTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGaaa                                         <  2:477199/140‑1 (MQ=255)                                                                                         tcTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGaaa                                         >  1:142526/1‑140 (MQ=255)                                                                                         tcTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGaaa                                         >  2:151613/1‑140 (MQ=255)                                                                                         tcTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGaaa                                         >  2:81954/1‑140 (MQ=255)                                                                                          cTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGaaaa                                        >  2:480630/1‑140 (MQ=255)                                                                                           tGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAAt                                       <  1:151613/140‑1 (MQ=255)                                                                                                 aTGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTg                                 >  1:140811/1‑140 (MQ=255)                                                                                                 aTGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTg                                 <  2:261503/140‑1 (MQ=255)                                                                                                    tCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTgcgg                              <  1:732338/140‑1 (MQ=255)                                                                                                       ttAGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGt                           >  2:872665/1‑140 (MQ=255)                                                                                                         aGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTgcgg                              >  2:479317/1‑135 (MQ=255)                                                                                                         aGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTgcgg                              >  1:116638/1‑135 (MQ=255)                                                                                                         aGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTgcgg                              <  2:116638/135‑1 (MQ=255)                                                                                                         aGCGGCGGCACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTgcgg                              <  1:479317/135‑1 (MQ=255)                                                                                                             gcggcACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGGTCAGACGGTTACCGGCGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCact                     >  1:292149/1‑140 (MQ=255)                                                                                                             gcggcACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCact                     <  2:469083/140‑1 (MQ=255)                                                                                                             gcggcACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCact                     >  2:731267/1‑140 (MQ=255)                                                                                                                 cACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCact                     <  2:457349/136‑1 (MQ=255)                                                                                                                 cACTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCact                     >  1:457349/1‑136 (MQ=255)                                                                                                                  aCTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACttt                <  1:566801/140‑1 (MQ=255)                                                                                                                  aCTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACttt                <  1:299396/140‑1 (MQ=255)                                                                                                                   cTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTg               <  2:59583/140‑1 (MQ=255)                                                                                                                   cTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTg               <  1:480630/140‑1 (MQ=255)                                                                                                                   cTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTg               <  2:819927/140‑1 (MQ=255)                                                                                                                   cTTTGATGATCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTg               <  2:849058/140‑1 (MQ=255)                                                                                                                      tgatgaTCGGCAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAg            <  2:666918/140‑1 (MQ=255)                                                                                                                                cAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAGCCGGATAtt   <  2:820868/139‑1 (MQ=255)                                                                                                                                cAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAGCCGGATAtt   >  1:820868/1‑139 (MQ=255)                                                                                                                                cAAGCGAAAGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAGCCGGATATTg  >  1:581895/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    GCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTTAGCGGCGGCACTTTGATGATCGGCAAGCGAACGCCAGGCGAGACTGAGCGGGATCAGACGGCTACCGGTGTTGATTACGCCACCCCACGGTTCATCGACTAACGCATCAAGATTGTTGCGGAAAATCTCCTGCGGCGTTTCACTACTTTGTAGCCGGATATTG  >  NC_000913/2331014‑2331280

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