Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	711,724	C→T	A59T (GCG→ACG)	ybfE ←	LexA‑regulated protein, CopB family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	711,724	0	C	T	100.0%	226.8 / NA	70	A59T (GCG→ACG)	ybfE	LexA‑regulated protein, CopB family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/32);  total (38/32)

GAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCG  >  NC_000913/711593‑711857                                                                                                                                    |                                                                                                                                      gAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGc                                                                                                                                >  1:1459484/1‑139 (MQ=255) gAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGc                                                                                                                                >  1:375608/1‑139 (MQ=255)   aaGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCtt                                                                                                                              >  1:903407/1‑139 (MQ=255)   aaGTGGAAATTTAAACGATACCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCtt                                                                                                                              >  1:1146117/1‑139 (MQ=255)             ttAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGaca                                                                                                                    <  2:977337/139‑1 (MQ=255)             ttAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGaca                                                                                                                    <  1:539195/139‑1 (MQ=255)                   gATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                              <  2:209879/139‑1 (MQ=255)                   gATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                              <  2:634551/139‑1 (MQ=255)                     tGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAg                                                                                                            >  2:172152/1‑139 (MQ=255)                     tGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAg                                                                                                            >  2:241991/1‑139 (MQ=255)                      gCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAgg                                                                                                           <  1:791338/139‑1 (MQ=255)                       cccTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                              <  2:232529/135‑1 (MQ=255)                       cccTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                              >  1:232529/1‑135 (MQ=255)                            aCTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACg                                                                                                      <  2:568500/138‑1 (MQ=255)                            aCTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACg                                                                                                      >  1:568500/1‑138 (MQ=255)                                 gCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACttt                                                                                                >  1:1150209/1‑139 (MQ=255)                                 gCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACttt                                                                                                >  1:1353611/1‑139 (MQ=255)                                    gcgcCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTc                                                                                             >  2:1352938/1‑139 (MQ=255)                                     cgccgcCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGt                                                                                                                          >  2:57506/1‑109 (MQ=255)                                     cgccgcCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGt                                                                                                                          <  1:57506/109‑1 (MQ=255)                                      gccgccAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTcgc                                                                                           >  1:1085576/1‑139 (MQ=255)                                      gccgccAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTcgc                                                                                           >  2:713747/1‑139 (MQ=255)                                               tgctgcATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCt                                                                                  <  1:172152/139‑1 (MQ=255)                                                      ttAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                           <  2:1459484/139‑1 (MQ=255)                                                      ttAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                           <  1:1352938/139‑1 (MQ=255)                                                            tctcTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACt                                                                                                  >  2:1174273/1‑110 (MQ=255)                                                            tctcTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACt                                                                                                  <  1:1174273/110‑1 (MQ=255)                                                              tctTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTaa                                                                   >  2:970826/1‑139 (MQ=255)                                                               ctTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAAt                                                                  >  1:81448/1‑139 (MQ=255)                                                                    aTCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCa                                                             >  2:1474449/1‑139 (MQ=255)                                                                         ttCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAgg                                                                                                           <  2:3717/88‑1 (MQ=255)                                                                         ttCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAgg                                                                                                           >  1:3717/1‑88 (MQ=255)                                                                          tCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                      >  1:582417/1‑76 (MQ=255)                                                                          tCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                      <  2:582417/76‑1 (MQ=255)                                                                           cGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACg                                                               >  2:399321/1‑130 (MQ=255)                                                                           cGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACg                                                               <  1:399321/130‑1 (MQ=255)                                                                           cGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTc                                                      >  1:1323109/1‑139 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGac                                                                                                                     <  2:105699/68‑1 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGac                                                                                                                     >  1:105699/1‑68 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                      <  2:607995/67‑1 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                      >  1:607995/1‑67 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGcc                                                                                                         >  1:1064324/1‑80 (MQ=255)                                                                                   cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGcc                                                                                                         <  2:1064324/80‑1 (MQ=255)                                                                                     cATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGaa                                            <  1:749116/139‑1 (MQ=255)                                                                                        atTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGc                                         >  2:936834/1‑139 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAg                                                                                                                           >  2:113306/1‑53 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAg                                                                                                                           <  1:113306/53‑1 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTa                                                                                                    <  1:1476434/76‑1 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTa                                                                                                    >  2:1476434/1‑76 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGa                                             >  1:167656/1‑131 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGa                                             <  2:167656/131‑1 (MQ=255)                                                                                            cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGAt                                     <  2:1323109/139‑1 (MQ=255)                                                                                                aCTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAg                                                                                                            >  1:534866/1‑64 (MQ=255)                                                                                                aCTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAg                                                                                                            <  2:534866/64‑1 (MQ=255)                                                                                                aCTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGt                                 >  2:171700/1‑139 (MQ=255)                                                                                                    cGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGtttt                                                                                      >  2:964291/1‑82 (MQ=255)                                                                                                    cGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGtttt                                                                                      <  1:964291/82‑1 (MQ=255)                                                                                                      ccAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                           >  2:390523/1‑91 (MQ=255)                                                                                                      ccAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                           <  1:390523/91‑1 (MQ=255)                                                                                                           tCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGgtcgta                      <  2:903407/139‑2 (MQ=255)                                                                                                            cGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTcc                     <  2:1085576/139‑1 (MQ=255)                                                                                                              ttCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCg                   >  2:380178/1‑139 (MQ=255)                                                                                                                cAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGt                 <  1:713747/139‑1 (MQ=255)                                                                                                                    gcCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGAc             >  1:1454839/1‑139 (MQ=255)                                                                                                                        cGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTc         >  2:1149841/1‑139 (MQ=255)                                                                                                                         gACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTcgcg                                                                                          <  1:297151/57‑1 (MQ=255)                                                                                                                         gACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTcgcg                                                                                          >  2:297151/1‑57 (MQ=255)                                                                                                                               ttCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGcgcg  <  1:741730/139‑1 (MQ=255)                                                                                                                                tCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGAt                                     >  2:513916/1‑103 (MQ=255)                                                                                                                                tCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGAt                                     <  1:513916/103‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      GAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCG  >  NC_000913/711593‑711857

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