Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,785	C→T	*239* (TAG→TAA)	nanC ←	N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,785	0	C	T	100.0%	227.7 / NA	70	*239* (TAG→TAA)	nanC	N‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/37);  total (33/37)

GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGT  >  NC_000913/4538649‑4538914                                                                                                                                         |                                                                                                                                  gtCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                 >  1:1427229/1‑139 (MQ=255)   cGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               <  1:105640/139‑1 (MQ=255)          aaTTGCTCCGGTACCACTTTTAAATGGTACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                        <  1:1647401/139‑1 (MQ=255)          aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                        <  2:110211/139‑1 (MQ=255)             tGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                                     <  1:14278/139‑1 (MQ=255)                    gTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                   <  2:768890/134‑1 (MQ=255)                    gTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                   >  1:768890/1‑134 (MQ=255)                     tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                             <  2:1409130/139‑1 (MQ=255)                          cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                          <  2:531522/121‑1 (MQ=255)                          cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                          >  1:531522/1‑121 (MQ=255)                          cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGa                                                                                                        >  1:888763/1‑139 (MQ=255)                          cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGCTTTCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACGGTATCCTGTCTATTACAGTTTAAATGGCACACCAATGCGa                                                                                                        >  1:690475/1‑139 (MQ=255)                                    gCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtttt                                                                                              <  1:1595600/139‑1 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTc                                                                                             >  2:421147/1‑139 (MQ=255)                                           aGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaa                                                                                       <  2:1427229/139‑1 (MQ=255)                                             tttCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                        >  1:861493/1‑104 (MQ=255)                                             tttCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                        <  2:861493/104‑1 (MQ=255)                                               tCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTa                                                                                   >  1:1515954/1‑139 (MQ=255)                                                       tACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTAtct                                                                                >  2:1480582/1‑134 (MQ=255)                                                       tACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCCTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATCAATTAtct                                                                                <  1:1480582/134‑1 (MQ=255)                                                             cgcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                                     <  2:1515954/139‑1 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGttt                                                                                                                           >  1:317481/1‑84 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGttt                                                                                                                           <  2:317481/84‑1 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTaa                                                                    >  1:72817/1‑139 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                         >  2:1220555/1‑131 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                         <  1:1220555/131‑1 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                 <  2:690475/139‑1 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                 >  1:1618590/1‑139 (MQ=255)                                                                   tGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             >  1:1029452/1‑125 (MQ=255)                                                                   tGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             <  2:1029452/125‑1 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTTTCCTGTCCATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  <  2:888763/139‑1 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaa                                                                                       >  2:912095/1‑98 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaa                                                                                       <  1:912095/98‑1 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                              <  2:1618590/139‑1 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                              >  2:962941/1‑139 (MQ=255)                                                                                      atGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCa                                            >  1:509258/1‑139 (MQ=255)                                                                                       tGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           <  1:421147/139‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaa                                                                                       >  1:1478432/1‑87 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaa                                                                                       <  2:1478432/87‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                             <  1:104756/129‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                             <  1:719734/129‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                             >  2:719734/1‑129 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTCACGGTCAAGGTAGTc                                             >  2:104756/1‑129 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               <  1:720546/46‑1 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               >  2:720546/1‑46 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           <  1:1344907/130‑1 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           >  2:1344907/1‑130 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACCGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCATGACGGTCAAGGTAGTCATATTCtatgta                                 <  1:962941/139‑1 (MQ=255)                                                                                                     gTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           <  2:1561327/125‑1 (MQ=255)                                                                                                     gTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           >  1:1561327/1‑125 (MQ=255)                                                                                                       tATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             >  1:1044073/1‑89 (MQ=255)                                                                                                       tATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             <  2:1044073/89‑1 (MQ=255)                                                                                                              ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtt                                                                                                <  2:457602/63‑1 (MQ=255)                                                                                                              ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtt                                                                                                >  1:457602/1‑63 (MQ=255)                                                                                                              ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                       >  2:262762/1‑104 (MQ=255)                                                                                                              ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                       <  1:262762/104‑1 (MQ=255)                                                                                                                aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                  >  1:1620206/1‑139 (MQ=255)                                                                                                                aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                  <  2:72817/139‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                                     <  2:82952/39‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                                     >  1:82952/1‑39 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                                     <  2:1525175/87‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                                     >  1:1525175/1‑87 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGc                 >  2:650738/1‑139 (MQ=255)                                                                                                                  tCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                <  2:1620206/139‑1 (MQ=255)                                                                                                                     ttGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCa             <  2:388422/139‑1 (MQ=255)                                                                                                                          aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  <  2:978165/97‑1 (MQ=255)                                                                                                                          aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  >  1:978165/1‑97 (MQ=255)                                                                                                                           cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTTTCTCTTCTGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGt       >  1:490842/1‑139 (MQ=255)                                                                                                                               ctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                             <  2:1170839/97‑1 (MQ=255)                                                                                                                               ctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                             >  1:1170839/1‑97 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGt  >  2:74183/1‑139 (MQ=255)                                                                                                                                      ctattaCAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTa      <  2:1557449/129‑1 (MQ=255)                                                                                                                                      ctattaCAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTa      >  1:1557449/1‑129 (MQ=255)                                                                                                                                         |                                                                                                                                  GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGT  >  NC_000913/4538649‑4538914

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