Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,357,417	0	C	T	24.4%	160.9 / 50.9	86	R172C (CGT→TGT)	yfaD	transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base T (11/10);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GTTTGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTAGAGGCTTCACT  >  NC_000913/2357320‑2357515                                                                                                   |                                                                                                   gTTTGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGt                                                                                                  >  2:1355119/1‑100 (MQ=255) gTTTGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGt                                                                                                  >  2:412909/1‑100 (MQ=255)   ttGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGtat                                                                                                >  1:2073511/1‑100 (MQ=255)   ttGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGtat                                                                                                <  2:2162942/100‑1 (MQ=255)    tGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGtata                                                                                               <  1:1160063/100‑1 (MQ=255)     gACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAg                                                                                              >  1:1287615/1‑100 (MQ=255)      aCGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGc                                                                                             >  2:1017513/1‑100 (MQ=255)      aCGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGc                                                                                             <  1:1017515/100‑1 (MQ=255)      aCGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGc                                                                                             <  2:2030441/100‑1 (MQ=255)      aCGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGc                                                                                             >  1:1604277/1‑100 (MQ=255)       cGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGct                                                                                            >  1:1093947/1‑100 (MQ=255)       cGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGct                                                                                            >  1:2009189/1‑100 (MQ=255)        gATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGctc                                                                                           <  2:2127/100‑1 (MQ=255)        gATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGctc                                                                                           >  2:91995/1‑100 (MQ=255)           ccGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTctgc                                                                                        >  2:393111/1‑100 (MQ=255)            cGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTctgct                                                                                       >  1:1464720/1‑100 (MQ=255)            cGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTctgct                                                                                       >  2:1447747/1‑100 (MQ=255)            cGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTctgct                                                                                       >  1:2298976/1‑100 (MQ=255)             gCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTctgctg                                                                                      >  2:66621/1‑100 (MQ=255)                  ttGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACt                                                                                 <  1:1617043/100‑1 (MQ=255)                   tGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTg                                                                                >  2:2003316/1‑100 (MQ=255)                    gCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGa                                                                               >  2:402154/1‑100 (MQ=255)                     cccGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGAt                                                                              <  1:346009/100‑1 (MQ=255)                     cccGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGAt                                                                              >  1:1489734/1‑100 (MQ=255)                       cGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTc                                                                            <  1:1319389/100‑1 (MQ=255)                        ggCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCa                                                                           <  1:462991/100‑1 (MQ=255)                          cttccgatctccGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCaaa                                                                         <  2:3261/90‑1 (MQ=255)                            gTTGTACACCGAAGCTTTTCTGTTGGTGGATAG‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAg                                                                       <  2:18309/100‑1 (MQ=255)                            cgatctcacCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAg                                                                       <  2:671576/94‑1 (MQ=255)                              tGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCa                                                                     >  2:1683717/1‑100 (MQ=255)                               gTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                                    >  2:1190838/1‑100 (MQ=255)                                tACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                                   >  2:1353286/1‑100 (MQ=255)                                  cacCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCagatcggaaga                                                                 >  2:1101142/1‑90 (MQ=255)                                  cacCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAgagatcg                                                                 >  1:671577/1‑94 (MQ=255)                                   acCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTc                                                                <  2:926131/100‑1 (MQ=255)                                   acCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTc                                                                <  2:927248/100‑1 (MQ=255)                                    ccGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAgatcggaag                                                               >  1:3261/1‑92 (MQ=255)                                    ccGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTcg                                                               <  1:2126345/100‑1 (MQ=255)                                    ccGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTcg                                                               >  1:227309/1‑100 (MQ=255)                                      gAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTcgcg                                                             <  1:2323590/100‑1 (MQ=255)                                        aGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGAc                                                           <  2:763011/100‑1 (MQ=255)                                        aGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGAc                                                           <  2:1737442/100‑1 (MQ=255)                                            tttCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAATTGATTCAAAAGCATATTCGCGACcgcg                                                       <  2:2073504/100‑1 (MQ=255)                                               tacgatggtgATATCCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGAt                                                     <  1:497502/91‑1 (MQ=255)                                               cTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGAtt                                                    >  1:647173/1‑100 (MQ=255)                                               cTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGAtt                                                    <  2:1147759/100‑1 (MQ=255)                                                 gTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                                  >  1:1723444/1‑100 (MQ=255)                                                 gTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                                  >  1:2025275/1‑100 (MQ=255)                                                 gTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                                  >  1:815555/1‑100 (MQ=255)                                                  ttGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTaa                                                 <  1:2003323/100‑1 (MQ=255)                                                   atctggATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAt                                                <  1:702565/97‑1 (MQ=255)                                                    ggtggATAT‑CCCCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTCATGCAAAAGCATATTCGCGACCGCGATTTAATc                                               <  2:1534668/100‑1 (MQ=255)                                                    ggtggATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATc                                               <  2:2009182/100‑1 (MQ=255)                                                      tggATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCgg                                             <  1:402156/100‑1 (MQ=255)                                                       ggATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtagat                                            >  2:702564/1‑96 (MQ=255)                                                        gATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATGCGCGACCGCGATTTAATCGGCa                                           <  1:1447751/100‑1 (MQ=255)                                                            t‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGt                                       >  2:1497225/1‑100 (MQ=255)                                                              cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTc                                      <  1:412911/100‑1 (MQ=255)                                                              cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTc                                      <  1:1353289/100‑1 (MQ=255)                                                              cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTc                                      >  1:2113474/1‑100 (MQ=255)                                                              cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTc                                      >  2:1330288/1‑100 (MQ=255)                                                               aCCATCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCg                                     >  1:911185/1‑100 (MQ=255)                                                               aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCg                                     <  1:1355122/100‑1 (MQ=255)                                                                 cATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGAc                                   <  2:1317231/100‑1 (MQ=255)                                                                   tCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAg                                 <  2:1022262/100‑1 (MQ=255)                                                                   tCGTACCTGACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAg                                 <  2:236151/100‑1 (MQ=255)                                                                       aCCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATc                             <  1:918015/100‑1 (MQ=255)                                                                        ccTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATca                            >  2:1962020/1‑100 (MQ=255)                                                                         cTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATcac                           >  1:867531/1‑100 (MQ=255)                                                                         cTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATcac                           >  2:1293291/1‑100 (MQ=255)                                                                           gACGATGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATcacca                         <  1:365622/100‑1 (MQ=255)                                                                           gACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATcacca                         >  2:2040523/1‑100 (MQ=255)                                                                           gACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATcacca                         <  2:227307/100‑1 (MQ=255)                                                                             cGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACg                       >  2:1739392/1‑100 (MQ=255)                                                                                tGAGATCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCtt                    >  2:205336/1‑100 (MQ=255)                                                                                   gaTCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTg                 >  2:923789/1‑100 (MQ=255)                                                                                   gaTCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTg                 <  2:530542/100‑1 (MQ=255)                                                                                   gaTCATGCAACATCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGACATGGTCGACAGGATCACCACGCTTTTg                 >  1:1266488/1‑100 (MQ=255)                                                                                    aTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTgg                >  1:1154620/1‑100 (MQ=255)                                                                                        tGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTag            <  1:66623/100‑1 (MQ=255)                                                                                        tGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTag            <  2:1604272/100‑1 (MQ=255)                                                                                          cAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTagag          >  2:944399/1‑100 (MQ=255)                                                                                               tCGGTGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTAGAGGCTTc     <  2:647172/100‑1 (MQ=255)                                                                                                 ggCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACGGGATCACCACGCTTTTGGTTAGAGGCTTCAc   >  1:1471245/1‑100 (MQ=255)                                                                                                 ggCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTAGAGGCTTCAc   >  2:953397/1‑100 (MQ=255)                                                                                                  gCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTAGAGGCTTCACt  >  1:1945014/1‑100 (MQ=255)                                                                                                   |                                                                                                   GTTTGACGATCCGCAACTTGCCCGGCAGTTGTACACCGAAGCTTTTCTGTTGGTGGATAT‑CACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAATCGGCATGGTCGACAGGATCACCACGCTTTTGGTTAGAGGCTTCACT  >  NC_000913/2357320‑2357515

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