Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	300,144	C→A	R264L (CGT→CTT)	paoC ←	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	300,144	0	C	A	100.0%	134.5 / NA	42	R264L (CGT→CTT)	paoC	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/25);  total (17/25)

AGACCAGCTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGACGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGCCAGAT  >  NC_000913/300011‑300279                                                                                                                                      |                                                                                                                                        aGACCAGCTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCa                                                                                                                                   <  1:1004808/140‑1 (MQ=255)  gACCAGCTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCAc                                                                                                                                  <  1:764971/140‑1 (MQ=255)       gCTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTTAcccccc                                                                                                                             >  1:1403515/1‑136 (MQ=255)        cTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTAACCGGAAGTTTTACCGCCCg                                                                                                                            >  2:434435/1‑140 (MQ=255)         tttCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCgg                                                                                                                           <  1:1164885/140‑1 (MQ=255)                 gagaTAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGa                                                                                                                    <  1:1306575/139‑1 (MQ=255)                 gagaTAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGa                                                                                                                    >  2:1306575/1‑139 (MQ=255)                  agaTAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTTACCGGGAGTTTCACCGCTCGGGCGGCGAgg                                                                                                                  >  2:545263/1‑140 (MQ=255)                       gCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCCCGGGCGGCGAGGgccgc                                                                                                             >  2:805681/1‑140 (MQ=255)                       gCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGGGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACCTTTACCGGGAGTTTCACCGCTCGGGCGGCGAGGggcgc                                                                                                             >  2:891541/1‑140 (MQ=255)                                       tctGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCtt                                                                                             >  2:771261/1‑140 (MQ=255)                                         tGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTtct                                                                                           <  2:654924/140‑1 (MQ=255)                                           gTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCa                                                                                         >  2:937391/1‑140 (MQ=255)                                                cACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAAc                                                                                    >  1:614259/1‑140 (MQ=255)                                                 aCCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACa                                                                                   >  2:1178342/1‑140 (MQ=255)                                                   cGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCAGCACTTTTACCGGGAGGTTCACCGCCCGGGGGGCGAAGGCCGGCCGCAGGGCAACCCCTCTCCCGAACAGc                                                                                 >  2:705967/1‑140 (MQ=255)                                                     aTACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGGTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCCCTTCTCAGGAACAGCtt                                                                               >  1:593108/1‑140 (MQ=255)                                                     aTACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACCTTTACCGGAAGTTTCAACGCTCGGGGGGCGAGGGCCCCCAACCGCGCATTGCTTCTTAAGAAAAgact                                                                               >  2:1153628/1‑137 (MQ=255)                                                             gTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCCCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCaa                                                                       >  1:372147/1‑140 (MQ=255)                                                                 ggAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAAccc                                                                   <  1:1153628/139‑1 (MQ=255)                                                                  gaggggTGGGGGGGCGGTGCGTGGTGTTATTGGGAATGGAGGGGCGGGGGAGCATCACTTCAACCGGAAGTTTCACCGCTCGGGCGGCGGGGGCCGCCGGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCt                                                                  <  2:705292/137‑1 (MQ=255)                                                                  gAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCt                                                                  <  2:330071/140‑1 (MQ=255)                                                                       gTGGGGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTccgcc                                                             <  1:705967/140‑1 (MQ=255)                                                                       gTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTACCCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTccgcc                                                             <  2:1193029/140‑1 (MQ=255)                                                                       gTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTccgcc                                                             <  2:372147/140‑1 (MQ=255)                                                                        tGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTccgccg                                                            >  2:404802/1‑140 (MQ=255)                                                                        tGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCCTCACCTTAACCGGAAGGTTTACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCCGGAACAGCTTGCCGCCGAACCCTccgccc                                                            >  2:752477/1‑139 (MQ=255)                                                                        tGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTTACCGCTCGGGCGGCGGGGGCCGCCCGCCACGCATCGCTTCTCAGGAACAACTTTGCGCCAAACCCTcccccc                                                            >  2:382762/1‑139 (MQ=255)                                                                        tGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTTACCCCCCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCCTGCCGCCAAACCCTccgccg                                                            >  2:1024938/1‑140 (MQ=255)                                                                        tGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTccgccg                                                            >  2:1383200/1‑140 (MQ=255)                                                                           cggggcggTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGata                                                         <  1:434435/140‑1 (MQ=255)                                                                             gggcggTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGAtata                                                       >  1:146093/1‑140 (MQ=255)                                                                             gggcggTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGAtata                                                       <  2:614259/140‑1 (MQ=255)                                                                               gcggTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACg                                                     <  2:146093/140‑1 (MQ=255)                                                                                gggTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACgg                                                    <  1:382762/139‑1 (MQ=255)                                                                                  gTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGGGCATCACTTTTAACGGGAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACgggg                                                  >  2:472180/1‑140 (MQ=255)                                                                                  gTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACgggg                                                  <  1:891541/140‑1 (MQ=255)                                                                                  gTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACgggg                                                  >  1:205056/1‑140 (MQ=255)                                                                                      gtggtgTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGat                                              <  1:692793/140‑1 (MQ=255)                                                                                                    aaTAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATtct                                >  1:29771/1‑140 (MQ=255)                                                                                                      tAGAGGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTcc                              >  2:381418/1‑140 (MQ=255)                                                                                                       agagGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCa                             <  2:793751/140‑1 (MQ=255)                                                                                                         agGGGCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACg                           >  2:354610/1‑140 (MQ=255)                                                                                                         agGGGCGGGGGAGCATCACTTTAACCCGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAAACCTCCGCCGATATACGGGGAGGTAATACGCACATTCTCCCCg                           >  2:683882/1‑140 (MQ=255)                                                                                                          ggggCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACgg                          <  1:752477/140‑1 (MQ=255)                                                                                                           gggCGGGGGAGCATCAATTTAACCGGAAGTTTTACCGCTCGGGCGGCGAGGGCCGCCAGCAGGGCATCGCTTCTCAGGAACAGCTTGCCCCCAAACCCTCCGGCGATATACGGGGAGATAAAAAGACCATTCTCCACggg                         >  2:1310506/1‑140 (MQ=255)                                                                                                           gggCGGGGAAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACggg                         <  1:683882/140‑1 (MQ=255)                                                                                                            ggCGGGGGAGCATCGCTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGa                        <  2:1403515/140‑1 (MQ=255)                                                                                                            ggCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATtctc                               >  1:486991/1‑133 (MQ=255)                                                                                                            ggCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATtctc                               <  2:486991/133‑1 (MQ=255)                                                                                                            ggCGGGGGAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGa                        >  2:88088/1‑140 (MQ=255)                                                                                                                ggggAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACtt                    >  1:317932/1‑140 (MQ=255)                                                                                                                ggggAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCACTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACtt                    >  1:674867/1‑140 (MQ=255)                                                                                                                   gAGCATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCa                 <  1:354610/140‑1 (MQ=255)                                                                                                                      cATCACTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCg              >  1:336437/1‑140 (MQ=255)                                                                                                                           cTTTAACCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGttttt         >  1:633525/1‑140 (MQ=255)                                                                                                                             ttAACCGGAAGTTTCACCGCCCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGc       <  2:317932/140‑1 (MQ=255)                                                                                                                                aCCGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGCCAg    >  2:345514/1‑140 (MQ=255)                                                                                                                                  cGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGCCAGAt  <  2:29771/140‑1 (MQ=255)                                                                                                                                  cGGAAGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGCCAGAt  <  1:381418/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        AGACCAGCTTTCATGTGAGATAGCGGTGATTTTCCCGCTCTGGTCGGCACCGATACGCAAGTGCTGAAGGGTGGCGGGGCGGTGCGTGGTGTTATTGGGAATAGAGGGGCGGGGGAGCATCACTTTAACCGGACGTTTCACCGCTCGGGCGGCGAGGGCCGCCAGCAGCGCATCGCTTCTCAGGAACAGCTTGCCGCCAAACCCTCCGCCGATATACGGGGAGATAATACGCACATTCTCCACGGGAACTTTCAGCGTTTTTGCCAGAT  >  NC_000913/300011‑300279

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