Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,330,165	A→G	E278E (GAA→GAG)	sohB →	S49 peptidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,330,165	0	A	G	100.0%	216.7 / NA	61	E278E (GAA→GAG)	sohB	S49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/33);  total (28/33)

TTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAACACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGA  >  NC_000913/1330034‑1330291                                                                                                                                    |                                                                                                                               ttGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTa                                                                                                                        <  2:508244/140‑1 (MQ=255)      gggTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAgcac                                                                                                                             <  1:869684/130‑1 (MQ=255)      gggTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAgcac                                                                                                                             >  2:869684/1‑130 (MQ=255)       ggTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGaca                                                                                                                  >  1:562834/1‑140 (MQ=255)        gTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGacaa                                                                                                                 >  2:91276/1‑140 (MQ=255)          gAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGacaaca                                                                                                               <  1:412318/140‑1 (MQ=255)            aaaTACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAgg                                                                                                             <  1:203995/140‑1 (MQ=255)            aaaTACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAgg                                                                                                             <  2:789393/140‑1 (MQ=255)            aaaTACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAgg                                                                                                             <  2:485893/140‑1 (MQ=255)                           gCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGcc                                                                                              >  1:99614/1‑140 (MQ=255)                            cGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGcc                                                                                              <  1:951837/139‑1 (MQ=255)                            cGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGcc                                                                                              >  2:951837/1‑139 (MQ=255)                               gagaAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGt                                                                                          >  2:725168/1‑140 (MQ=255)                                     ttCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTgg                                                                                           <  2:696466/133‑1 (MQ=255)                                     ttCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTgg                                                                                           >  1:696466/1‑133 (MQ=255)                                     ttCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTtgatga                                                                                    >  1:892932/1‑140 (MQ=255)                                         gcgAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTtgatg                                                                                     >  2:156720/1‑135 (MQ=255)                                         gcgAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTtgatg                                                                                     <  1:156720/135‑1 (MQ=255)                                         gcgAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATc                                                                                >  1:781321/1‑140 (MQ=255)                                              gagTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTgg                                                                                           <  2:802556/124‑1 (MQ=255)                                              gagTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTgg                                                                                           >  1:802556/1‑124 (MQ=255)                                              gagTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAacac                                                                           <  1:565386/140‑1 (MQ=255)                                                   gAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTAAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGtg                                                                      >  2:870002/1‑140 (MQ=255)                                                    aaCGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGtga                                                                     <  1:770644/140‑1 (MQ=255)                                                     aCGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGtgat                                                                    <  1:823987/140‑1 (MQ=255)                                                           cGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGtt                                                              >  2:500136/1‑140 (MQ=255)                                                                   ttgttTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAg                                                      >  1:251504/1‑140 (MQ=255)                                                                           aGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATgg                                              >  2:166329/1‑140 (MQ=255)                                                                            gATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGa                                             <  1:166329/140‑1 (MQ=255)                                                                            gATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGa                                             <  1:178772/140‑1 (MQ=255)                                                                             aTTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGaa                                            >  1:551983/1‑140 (MQ=255)                                                                               tttGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAgg                                          <  2:740066/140‑1 (MQ=255)                                                                                ttGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGa                                             >  1:641193/1‑136 (MQ=255)                                                                                ttGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGa                                             <  2:641193/136‑1 (MQ=255)                                                                                ttGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGc                                         >  1:694473/1‑140 (MQ=255)                                                                                 tgtgAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGcc                                        <  2:201293/140‑1 (MQ=255)                                                                                 tgtgAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGcc                                        <  2:551983/140‑1 (MQ=255)                                                                                        cGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGt                                 <  1:725168/140‑1 (MQ=255)                                                                                          tATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTgg                               <  2:820486/140‑1 (MQ=255)                                                                                          tATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTgg                               <  1:682676/140‑1 (MQ=255)                                                                                           aTGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGt                              <  2:694473/140‑1 (MQ=255)                                                                                           aTGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGt                              <  1:500136/140‑1 (MQ=255)                                                                                           aTGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGt                              <  1:902847/140‑1 (MQ=255)                                                                                            tGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTc                             >  1:274954/1‑140 (MQ=255)                                                                                            tGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTc                             >  2:868855/1‑140 (MQ=255)                                                                                               gTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAAt                          <  2:562834/140‑1 (MQ=255)                                                                                               gTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAAt                          <  1:833038/140‑1 (MQ=255)                                                                                                 ccTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGt                        <  2:99614/140‑1 (MQ=255)                                                                                                  cTTCTCTGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTa                       >  1:502339/1‑140 (MQ=255)                                                                                                     ctctGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGc                    >  1:450547/1‑140 (MQ=255)                                                                                                      tctGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCt                   <  2:781321/140‑1 (MQ=255)                                                                                                        tGGATATTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCtat                 >  2:690347/1‑140 (MQ=255)                                                                                                           atatTGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATg              <  1:91276/140‑1 (MQ=255)                                                                                                               tGAACAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGc          >  2:7409/1‑140 (MQ=255)                                                                                                                 aaCAGGTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGt        >  1:366268/1‑140 (MQ=255)                                                                                                                     ggTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAAc    <  1:690347/140‑1 (MQ=255)                                                                                                                     ggTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAAc    <  2:274954/140‑1 (MQ=255)                                                                                                                      gcggCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACg   >  2:454488/3‑140 (MQ=255)                                                                                                                      gTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACg   >  2:692772/1‑140 (MQ=255)                                                                                                                      gTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACg   >  2:85075/1‑140 (MQ=255)                                                                                                                      gTGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACg   >  1:960845/1‑140 (MQ=255)                                                                                                                       tGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGa  <  2:502339/140‑1 (MQ=255)                                                                                                                       tGGCAACGGGTGAGCACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGa  <  2:450547/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                               TTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAGATTTTGTGAAGCGTATGCGTCCTTCTCTGGATATTGAACAGGTGGCAACGGGTGAACACTGGTACGGACAACAGGCGGTAGAGAAAGGCCTGGTTGATGAAATCAACACCAGTGATGAAGTTATTCTTAGCCTGATGGAAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGA  >  NC_000913/1330034‑1330291

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