Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,616,195	A→G	E156G (GAA→GGA)	radA →	DNA repair protein RadA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,616,195	0	A	G	96.9%	221.2 / ‑1.0	65	E156G (GAA→GGA)	radA	DNA repair protein RadA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (26/37);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGG  >  NZ_CP009273/4616065‑4616329                                                                                                                                   |                                                                                                                                       cTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCa                                                                                                                                <  1:960010/139‑1 (MQ=255)     aaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                            <  2:1248776/139‑1 (MQ=255)     aaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                            <  2:599793/139‑1 (MQ=255)      aaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGaa                                                                                                                           <  2:231933/139‑1 (MQ=255)      aaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGaa                                                                                                                           >  2:285696/1‑139 (MQ=255)        cGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGaa                                                                                                                           >  2:996380/1‑137 (MQ=255)         gCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAg                                                                                                                        <  1:380889/139‑1 (MQ=255)           tgtgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGaa                                                                                                                           >  2:991968/1‑134 (MQ=255)             tgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTGACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATc                                                                                                                     <  2:856356/138‑1 (MQ=255)             tgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCt                                                                                                                    >  1:139774/1‑139 (MQ=255)                       cccAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGc                                                                                                          >  1:577866/1‑139 (MQ=255)                       cccAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGc                                                                                                          >  1:214233/1‑139 (MQ=255)                            cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCgaag                                                                                                     <  2:999240/139‑1 (MQ=255)                                  aaaaCGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAagag                                                                                                <  1:948671/138‑1 (MQ=255)                                      cGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAAcc                                                                                           >  2:159651/1‑139 (MQ=255)                                        cTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGa                                                                                         <  1:1238328/139‑1 (MQ=255)                                         tGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCg                                                                                          >  1:1165346/1‑137 (MQ=255)                                                   cGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGt                                                                               >  2:899711/1‑138 (MQ=255)                                                   cGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGt                                                                               >  1:1164009/1‑138 (MQ=255)                                                     gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAt                                                                            >  2:611490/1‑139 (MQ=255)                                                     gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAt                                                                            >  1:476858/1‑139 (MQ=255)                                                     gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAt                                                                            >  2:208015/1‑139 (MQ=255)                                                     gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAt                                                                            >  2:1177667/1‑139 (MQ=255)                                                     gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTCACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGTGCAACCGATGCTGATGGTAAt                                                                            >  2:345896/1‑139 (MQ=255)                                                         agagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAc                                                                        <  1:1177667/139‑1 (MQ=255)                                                         agagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAc                                                                        <  1:899711/139‑1 (MQ=255)                                                          gagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACt                                                                       <  2:300684/139‑1 (MQ=255)                                                               gCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGAt                                                                   <  1:116892/138‑1 (MQ=255)                                                               gCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATc                                                                  <  2:1164009/139‑1 (MQ=255)                                                                cTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATcc                                                                 <  2:214233/139‑1 (MQ=255)                                                                    aaCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCgaag                                                                                                     >  1:1001226/1‑99 (MQ=255)                                                                    aaCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCgaag                                                                                                     <  2:1001226/99‑1 (MQ=255)                                                                     aCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTg                                                            <  1:159651/139‑1 (MQ=255)                                                                          tGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCa                                                       >  2:306293/1‑139 (MQ=255)                                                                               aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  2:139774/139‑1 (MQ=255)                                                                               aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  2:149959/139‑1 (MQ=255)                                                                               aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  1:285696/139‑1 (MQ=255)                                                                               aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  2:476858/139‑1 (MQ=255)                                                                               aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  1:208015/139‑1 (MQ=255)                                                                                tGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGc                                                 >  1:975174/1‑139 (MQ=255)                                                                                  cgcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATTCAGGTGATGCATATggag                                                >  2:497049/1‑136 (MQ=255)                                                                                  cgcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAATCGTTCCAGGTGATGCATATGGCgg                                               <  2:626276/139‑1 (MQ=255)                                                                                    cgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  1:191393/134‑1 (MQ=255)                                                                                    cgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  >  2:191393/1‑134 (MQ=255)                                                                                    cgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGAt                                             >  2:203597/1‑139 (MQ=255)                                                                                            gCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGt                                      <  1:203597/138‑1 (MQ=255)                                                                                            gCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGt                                      <  2:1085937/138‑1 (MQ=255)                                                                                            gCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGt                                      <  2:320355/138‑1 (MQ=255)                                                                                             ccTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCa                                    <  1:497049/139‑1 (MQ=255)                                                                                             ccTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCa                                    <  1:306293/139‑1 (MQ=255)                                                                                                 ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAc                                                                        >  1:732336/1‑99 (MQ=255)                                                                                                 ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAc                                                                        <  2:732336/99‑1 (MQ=255)                                                                                                 ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCg                                 <  2:165888/138‑1 (MQ=255)                                                                                                 ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGc                                >  1:592461/1‑139 (MQ=255)                                                                                                       ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                            <  2:559410/137‑1 (MQ=255)                                                                                                       ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                            >  1:804943/1‑137 (MQ=255)                                                                                                          aCTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCg                       <  1:882423/139‑1 (MQ=255)                                                                                                               cAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGgcg                  <  2:577866/139‑1 (MQ=255)                                                                                                                      aaTATGTTGGCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGt             <  2:496137/137‑1 (MQ=255)                                                                                                                        tatGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  >  2:586880/1‑98 (MQ=255)                                                                                                                        tatGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                  <  1:586880/98‑1 (MQ=255)                                                                                                                             tgtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAAc    >  1:59311/1‑139 (MQ=255)                                                                                                                              gtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGa                                                                    <  1:24198/74‑1 (MQ=255)                                                                                                                              gtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGa                                                                    >  2:24198/1‑74 (MQ=255)                                                                                                                               tCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAAcgg  <  1:413218/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       CTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGG  >  NZ_CP009273/4616065‑4616329

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