Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	327,990	T→C	E248E (GAA→GAG)	betB ←	betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	327,990	0	T	C	100.0%	204.8 / NA	57	E248E (GAA→GAG)	betB	betaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/22);  total (35/22)

ATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCC  >  NC_000913/327842‑328136                                                                                                                                                     |                                                                                                                                                   aTTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCAcctc                                                                                                                                                  >  2:131379/1‑151 (MQ=255)    tCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCttt                                                                                                                                               >  1:408550/1‑151 (MQ=255)          gAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGa                                                                                                                                         <  1:334533/151‑1 (MQ=255)             cGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGa                                                                                                                                         >  1:21538/1‑148 (MQ=255)                 gACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGgcc                                                                                                                                  <  2:294711/151‑1 (MQ=255)                      gggTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGccgccg                                                                                                                              >  2:711256/1‑150 (MQ=255)                      gggTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGa                                                                                                                             <  1:718874/151‑1 (MQ=255)                        gTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGt                                                                                                                           <  1:183437/151‑1 (MQ=255)                          gCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTa                                                                                                                         <  2:381347/151‑1 (MQ=255)                             aTTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGcc                                                                                                                      >  2:806108/1‑151 (MQ=255)                               tGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCa                                                                                                                     <  2:656569/150‑1 (MQ=255)                                ggTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATc                                                                                                                   >  2:615723/1‑151 (MQ=255)                                   acacacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACt                                                                                                                >  2:509909/1‑151 (MQ=255)                                    cacacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACtt                                                                                                               <  2:599571/151‑1 (MQ=255)                                      cacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACtttt                                                                                                             >  1:252641/1‑151 (MQ=255)                                          tGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGc                                                                                                         >  2:140578/1‑151 (MQ=255)                                          tGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGc                                                                                                         >  2:58638/1‑151 (MQ=255)                                           gACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGcc                                                                                                        >  1:621929/1‑151 (MQ=255)                                           gACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTCTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGcc                                                                                                        >  1:737125/1‑151 (MQ=255)                                            aCCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCg                                                                                                       >  2:412622/1‑151 (MQ=255)                                                     gaagaaGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGac                                                                                               >  1:588148/1‑150 (MQ=255)                                                      aagaagTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGacac                                                                                             >  2:855125/1‑151 (MQ=255)                                                                ccATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTaa                                                                                    >  2:136705/1‑150 (MQ=255)                                                                      tGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCTCCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTTGCGACACCGCCGGTAAATGaca                                                                              >  2:692294/1‑150 (MQ=255)                                                                      tGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGacac                                                                             >  1:674376/1‑151 (MQ=255)                                                                       ggCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACt                                                                            >  1:10915/1‑151 (MQ=255)                                                                          gATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTg                                                                         >  2:591748/1‑151 (MQ=255)                                                                                ggcggcGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATg                                                                   >  2:690604/1‑151 (MQ=255)                                                                                   ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCCCCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGccc                                                                <  2:893267/151‑1 (MQ=255)                                                                                   ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGccc                                                                <  2:21516/151‑1 (MQ=255)                                                                                    gcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCg                                                               >  1:867739/1‑151 (MQ=255)                                                                                     cgagatcgagatCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCg                                                               <  1:695521/150‑1 (MQ=255)                                                                                           cgagaccCGCATCATCGAAAACGATCAGCGGTGATTTCCCCCCCAGTTCCATGGTCCCCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCGGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCt                                                         <  2:677849/150‑1 (MQ=255)                                                                                              gatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGt                                                     >  1:309676/1‑151 (MQ=255)                                                                                               atcCGCATCATCGAACACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTc                                                    >  2:521822/1‑151 (MQ=255)                                                                                                 ccGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTc                                                    >  2:406795/1‑149 (MQ=255)                                                                                                 ccGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAg                                                  <  1:509959/151‑1 (MQ=255)                                                                                                    catcatCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGat                                                >  1:523963/1‑150 (MQ=255)                                                                                                      tcatcGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGata                                               <  2:893310/149‑1 (MQ=255)                                                                                                       catcGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTg                                            >  2:120598/1‑151 (MQ=255)                                                                                                            aaaaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCgg                                       >  1:706477/1‑151 (MQ=255)                                                                                                              aaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtc                                     >  1:448058/1‑151 (MQ=255)                                                                                                               aCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCGGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtct                                    <  1:869085/151‑1 (MQ=255)                                                                                                                cGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtctc                                   >  1:656329/1‑151 (MQ=255)                                                                                                                 gATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTcc                                  <  2:756439/151‑1 (MQ=255)                                                                                                                  aTCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcg                                 >  1:677182/1‑151 (MQ=255)                                                                                                                   tCAGCGGTGATTTCCCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtctc                                   <  1:278012/148‑1 (MQ=255)                                                                                                                   tCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcg                                 <  1:475135/150‑1 (MQ=255)                                                                                                                        ggTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCa                           <  2:408524/151‑1 (MQ=255)                                                                                                                          tGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACg                         >  1:215343/1‑151 (MQ=255)                                                                                                                             tttACCGCCCAGTTCCAAGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGGTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGccc                      >  2:382960/1‑151 (MQ=255)                                                                                                                                    cccAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAac                >  2:127684/1‑150 (MQ=255)                                                                                                                                       aGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGCCACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGt            >  2:746933/1‑151 (MQ=255)                                                                                                                                        gTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGtt           >  1:143705/1‑151 (MQ=255)                                                                                                                                            cATGGTCCCCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAAt       <  1:833893/151‑1 (MQ=255)                                                                                                                                               ggTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACg    <  2:568812/151‑1 (MQ=255)                                                                                                                                                  cACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGcc  <  1:777116/150‑1 (MQ=255)                                                                                                                                                  cACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGcc  <  1:342035/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCC  >  NC_000913/327842‑328136

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