Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	468,115	A→C	Q235P (CAA→CCA)	cof →	thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	468,115	0	A	C	100.0%	222.2 / NA	59	Q235P (CAA→CCA)	cof	thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/29);  total (30/29)

AATAAAGGCGCTGCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCAACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCG  >  NC_000913/467973‑468262                                                                                                                                               |                                                                                                                                                    aaTAAAGGCGCTGCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTgcgcg                                                                                                                                             <  2:257777/151‑1 (MQ=255)  aTAAAGGCGCTGCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTgcgcgc                                                                                                                                            >  1:585299/1‑151 (MQ=255)             gCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGc                                                                                                                                 >  2:398947/1‑151 (MQ=255)             gCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGAAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGc                                                                                                                                 >  1:610065/1‑151 (MQ=255)               aTTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCAt                                                                                                                               >  2:829557/1‑151 (MQ=255)               aTTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCAt                                                                                                                               <  2:674175/151‑1 (MQ=255)                  gACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCAttt                                                                                                                             >  1:883820/1‑150 (MQ=255)                   aCGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTAc                                                                                                                           >  2:619956/1‑151 (MQ=255)                          tGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGAt                                                                                                                    >  2:136714/1‑151 (MQ=255)                            aCCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTg                                                                                                                  >  1:799664/1‑151 (MQ=255)                             cccAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTgg                                                                                                                 <  2:135199/151‑1 (MQ=255)                              ccAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTgg                                                                                                                 >  1:86889/1‑150 (MQ=255)                                 aCATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACa                                                                                                              <  2:592029/150‑1 (MQ=255)                                    tttAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGc                                                                                                          >  2:562793/1‑151 (MQ=255)                                     ttAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGc                                                                                                          <  2:550130/150‑1 (MQ=255)                                      tAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGcc                                                                                                         <  2:88082/150‑1 (MQ=255)                                        ggTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGaa                                                                                                      <  1:61611/151‑1 (MQ=255)                                            tATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATc                                                                                                   <  1:347500/150‑1 (MQ=255)                                              tCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAgg                                                                                                >  2:407394/1‑151 (MQ=255)                                                gTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCt                                                                                              >  1:861075/1‑151 (MQ=255)                                                   gcgcGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGt                                                                                            <  1:692089/150‑1 (MQ=255)                                                     gcgATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGtctc                                                                                         >  2:792559/1‑151 (MQ=255)                                                                gCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGa                                                                              <  1:891661/151‑1 (MQ=255)                                                                   tttGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGc                                                                           >  1:598453/1‑151 (MQ=255)                                                                     tGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCAc                                                                         >  2:380394/1‑151 (MQ=255)                                                                       gTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCActg                                                                       >  2:856665/1‑151 (MQ=255)                                                                        tGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCActgg                                                                      >  2:581125/1‑151 (MQ=255)                                                                         gATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCActggc                                                                     >  1:709536/1‑151 (MQ=255)                                                                                   aCGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATcc                                                           <  2:435700/151‑1 (MQ=255)                                                                                   aCGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATcc                                                           <  1:118355/151‑1 (MQ=255)                                                                                      aTCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCcaca                                                        <  2:773924/151‑1 (MQ=255)                                                                                        cgcgAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACAt                                                       <  1:195770/150‑1 (MQ=255)                                                                                          cgAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCt                                                     >  2:432932/1‑150 (MQ=255)                                                                                               tGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCtt                                                >  1:526258/1‑150 (MQ=255)                                                                                                  tAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTAtt                                             <  1:843244/150‑1 (MQ=255)                                                                                                    gTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTccc                                          <  1:859665/151‑1 (MQ=255)                                                                                                      cAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTccccc                                        >  1:877282/1‑151 (MQ=255)                                                                                                      cAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTcccc                                         <  2:927000/150‑1 (MQ=255)                                                                                                        gCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCg                                       >  1:889500/1‑150 (MQ=255)                                                                                                        gCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGa                                      <  1:526577/151‑1 (MQ=255)                                                                                                          gTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAAt                                    <  1:62589/151‑1 (MQ=255)                                                                                                           tCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAAt                                    >  2:665492/1‑150 (MQ=255)                                                                                                            cGGTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATa                                   <  1:213275/150‑1 (MQ=255)                                                                                                              gTAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACg                                <  2:83326/151‑1 (MQ=255)                                                                                                               tAGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACg                                >  1:180982/1‑150 (MQ=255)                                                                                                                aGCGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACg                                >  2:377086/1‑149 (MQ=255)                                                                                                                  cGGATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGAt                            <  2:220823/151‑1 (MQ=255)                                                                                                                    gATTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATcc                          >  1:543663/1‑151 (MQ=255)                                                                                                                     aTTTATTATGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATccc                         <  2:25120/151‑1 (MQ=255)                                                                                                                       ttattaTGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCt                        <  1:881413/150‑1 (MQ=255)                                                                                                                          ttatGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTcc                     >  2:1608/1‑150 (MQ=255)                                                                                                                           tatGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCa                    <  2:263428/150‑1 (MQ=255)                                                                                                                            atGGGCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGc                  >  2:632367/1‑151 (MQ=255)                                                                                                                              gggCAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCAc                <  2:140863/151‑1 (MQ=255)                                                                                                                                 cAATGCGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCgg             >  1:31192/1‑151 (MQ=255)                                                                                                                                      cGATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAAt        >  1:273052/1‑151 (MQ=255)                                                                                                                                       gATGCCGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAAtt       <  2:598472/151‑1 (MQ=255)                                                                                                                                            cGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCg  <  2:521641/151‑1 (MQ=255)                                                                                                                                            cGCCACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCATTCCAGCACCGGGCAATTGCCCg  <  2:467428/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    AATAAAGGCGCTGCATTGACGGTGCTGACCCAACATTTAGGTTTATCGTTGCGCGATTGCATGGCCTTTGGTGATGCGATGAACGATCGCGAAATGTTAGTCAGCGTCGGTAGCGGATTTATTATGGGCAATGCGATGCCGCAACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCG  >  NC_000913/467973‑468262

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