Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,087,870	A→G	F663S (TTT→TCT)	pgaB ←	poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,087,870	0	A	G	100.0%	231.8 / NA	65	F663S (TTT→TCT)	pgaB	poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/24);  total (41/24)

CGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGC  >  NC_000913/1087733‑1088004                                                                                                                                          |                                                                                                                                       cGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGAc                                                                                                                                      >  1:708553/1‑140 (MQ=255) cGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGAc                                                                                                                                      >  2:869824/1‑140 (MQ=255)       gaagaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTATTCGGATACCAGGCTGTTGAAGACTCAGGa                                                                                                                                >  2:743182/1‑140 (MQ=255)             aaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATa                                                                                                                          <  1:57050/140‑1 (MQ=255)                 aaCCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGt                                                                                                                      <  1:795797/140‑1 (MQ=255)                 aaCCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGt                                                                                                                      <  1:442065/140‑1 (MQ=255)                  aCCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCGTAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTc                                                                                                                     >  2:734646/1‑140 (MQ=255)                  aCCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTc                                                                                                                     >  1:709430/1‑140 (MQ=255)                  aCCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTc                                                                                                                     >  2:669329/1‑140 (MQ=255)                       atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTAttt                                                                                                                >  1:35757/1‑140 (MQ=255)                       atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTAttt                                                                                                                <  1:669329/140‑1 (MQ=255)                       atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTAttt                                                                                                                <  1:734646/140‑1 (MQ=255)                       atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTAttt                                                                                                                <  2:709430/140‑1 (MQ=255)                          atcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAg                                                                                                             <  2:35757/140‑1 (MQ=255)                                 cACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAggttggtt                                                                                                      >  2:206601/1‑140 (MQ=255)                                   ccAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTAt                                                                                                    >  1:399377/1‑140 (MQ=255)                                    cAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATg                                                                                                   >  2:819808/1‑140 (MQ=255)                                       agTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAg                                                                                                >  1:1045640/1‑140 (MQ=255)                                         tGAAAGGACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGaa                                                                                              >  2:541256/1‑140 (MQ=255)                                          gAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGaaa                                                                                             <  2:399377/140‑1 (MQ=255)                                           aaaGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGaaaa                                                                                            >  1:900830/1‑140 (MQ=255)                                              gTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTg                                                                                         <  1:910269/140‑1 (MQ=255)                                                acgctacgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTc                                                                                       <  1:1195316/140‑1 (MQ=255)                                                 cgctacgcACAGTGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCg                                                                                      >  2:473865/1‑140 (MQ=255)                                                  gctacgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCgg                                                                                     <  1:743182/140‑1 (MQ=255)                                                   ctacgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACAATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCggg                                                                                    <  1:1099871/140‑1 (MQ=255)                                                     acgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGat                                                                                  >  1:137075/1‑140 (MQ=255)                                                      cgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGata                                                                                 <  1:206601/140‑1 (MQ=255)                                                      cgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGata                                                                                 >  1:981380/1‑140 (MQ=255)                                                      cgcATAGGGGAATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGata                                                                                 <  1:244412/140‑1 (MQ=255)                                                           aGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAGAATTGTCGGGATAATAAc                                                                            <  2:1136734/140‑1 (MQ=255)                                                              ggTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCAt                                                                         >  1:1191773/1‑140 (MQ=255)                                                                  tGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGtt                                                                     >  1:786896/1‑140 (MQ=255)                                                                  tGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGtt                                                                     >  1:445817/1‑140 (MQ=255)                                                                       aaCACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCa                                                                >  2:336120/1‑140 (MQ=255)                                                                       aaCACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCa                                                                >  1:74412/1‑140 (MQ=255)                                                                       aaCACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCa                                                                <  1:473865/140‑1 (MQ=255)                                                                         cacCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACt                                                              <  1:247764/140‑1 (MQ=255)                                                                            cAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGttttt                                                                  <  1:258096/133‑1 (MQ=255)                                                                            cAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGttttt                                                                  >  2:258096/1‑133 (MQ=255)                                                                            cAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCa                                                           >  1:708560/1‑140 (MQ=255)                                                                             aGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCAt                                                          >  2:74210/1‑140 (MQ=255)                                                                              gACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCa                                                           <  2:717042/138‑1 (MQ=255)                                                                              gACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCa                                                           >  1:717042/1‑138 (MQ=255)                                                                              gACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCAtt                                                         >  2:769836/1‑140 (MQ=255)                                                                              gACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCAtt                                                         >  2:336688/1‑140 (MQ=255)                                                                              gACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCAtt                                                         >  1:261229/1‑140 (MQ=255)                                                                                 ataataTAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAg                                                      >  1:620052/1‑140 (MQ=255)                                                                                   aataTAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGtt                                                    >  2:1170578/1‑140 (MQ=255)                                                                                   aataTAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGtt                                                    >  2:608041/1‑140 (MQ=255)                                                                                       taAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTaa                                                <  2:137075/140‑1 (MQ=255)                                                                                       taAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTaa                                                <  2:74412/140‑1 (MQ=255)                                                                                             acgatacgatGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAg                                                      >  1:1198359/1‑128 (MQ=255)                                                                                             acgatacgatGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAg                                                      <  2:1198359/128‑1 (MQ=255)                                                                                               gatacgatGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCa                                        >  1:826691/1‑140 (MQ=255)                                                                                                           ttAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAg                               <  2:375250/137‑1 (MQ=255)                                                                                                           ttAATCATTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAg                               >  1:375250/1‑137 (MQ=255)                                                                                                                 aTTTTTCGGATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAgttgtt                      >  1:907223/1‑140 (MQ=255)                                                                                                                         gATACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCgaagaa              >  1:519364/1‑140 (MQ=255)                                                                                                                          aTACCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGaaa             <  1:869824/140‑1 (MQ=255)                                                                                                                            aCCAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATa           >  1:1110798/1‑140 (MQ=255)                                                                                                                             ccAGGCTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAg          >  1:548329/1‑140 (MQ=255)                                                                                                                                  cTGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGa     >  1:646639/1‑140 (MQ=255)                                                                                                                                   tGTTGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGAt    >  1:295361/1‑140 (MQ=255)                                                                                                                                     ttGAAGACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAATTAGCCTGATGc  >  2:1071917/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       CGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAGGTCTATTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGC  >  NC_000913/1087733‑1088004

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