Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,303	T→C	M201T (ATG→ACG)	cmoA →	carboxy‑SAM synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,303	0	T	C	100.0%	216.2 / NA	59	M201T (ATG→ACG)	cmoA	carboxy‑SAM synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/29);  total (30/29)

GGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAG  >  NC_000913/1953158‑1953431                                                                                                                                                  |                                                                                                                                 gggCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTgg                                                                                                                             >  1:797689/1‑151 (MQ=255)      gAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGaaa                                                                                                                          <  1:428405/149‑1 (MQ=255)      gAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAAc                                                                                                                        >  2:353228/1‑151 (MQ=255)           ccGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGAt                                                                                                                   >  1:755238/1‑151 (MQ=255)            cggcggCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATg                                                                                                                  <  1:142695/151‑1 (MQ=255)               cggcgCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTg                                                                                                               >  2:319351/1‑151 (MQ=255)                     gcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGAt                                                                                                         >  1:483077/1‑151 (MQ=255)                        ggTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTcc                                                                                                      <  1:194877/151‑1 (MQ=255)                        ggTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTc                                                                                                       <  1:319377/150‑1 (MQ=255)                         gttcTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACAGTTACAGCGAACTGGATATCAGCCAGAAGCGTAGCACGCTGGAAAACGTGATGCTGACCGATTcc                                                                                                      <  2:707581/147‑1 (MQ=255)                               tCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGaa                                                                                                >  2:100923/1‑150 (MQ=255)                                  gAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAAccc                                                                                            <  1:116248/151‑1 (MQ=255)                                   aaaaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGTTGACCGATTCCGTGGAAACCCa                                                                                           <  2:267274/151‑1 (MQ=255)                                   aaaaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCa                                                                                           <  1:533791/151‑1 (MQ=255)                                             tttCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGc                                                                                 >  2:430041/1‑151 (MQ=255)                                               tCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCt                                                                               >  2:839761/1‑151 (MQ=255)                                                  aaGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGc                                                                             >  2:555589/1‑150 (MQ=255)                                                  aaGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCa                                                                            <  2:902682/151‑1 (MQ=255)                                                    gATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATa                                                                          <  1:481213/151‑1 (MQ=255)                                                    gATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATa                                                                          <  2:110508/151‑1 (MQ=255)                                                         cAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCACAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGcc                                                                     >  2:796894/1‑151 (MQ=255)                                                           aaGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCgg                                                                   <  1:104222/151‑1 (MQ=255)                                                             gTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGtt                                                                 <  1:572510/151‑1 (MQ=255)                                                              ttGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACGGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGAAAACCCATAAAGCACGCCGGCATAACGCCGGttt                                                                >  2:885498/1‑151 (MQ=255)                                                              ttGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGttt                                                                >  1:274998/1‑151 (MQ=255)                                                              ttGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGttt                                                                >  2:35932/1‑151 (MQ=255)                                                                ggTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGtttt                                                               <  2:102966/150‑1 (MQ=255)                                                                  tGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCTGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTg                                                              >  1:79375/1‑149 (MQ=255)                                                                  tGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTg                                                              >  1:65793/1‑149 (MQ=255)                                                                   gAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGc                                                           >  1:33256/1‑151 (MQ=255)                                                                    aaCTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGc                                                           >  1:872454/1‑150 (MQ=255)                                                                    aaCTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCa                                                          <  2:755188/151‑1 (MQ=255)                                                                    aaCTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACACCGGTTTTGAGc                                                           >  1:257050/1‑150 (MQ=255)                                                                       tgttgtTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCAt                                                         >  2:721843/1‑149 (MQ=255)                                                                        gttgttCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATa                                                        <  2:878862/149‑1 (MQ=255)                                                                        gttgttCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGt                                                      >  1:673537/1‑151 (MQ=255)                                                                          tgttCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTg                                                     >  2:223361/1‑150 (MQ=255)                                                                                       ccacGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCa                                        <  2:311046/150‑1 (MQ=255)                                                                                        cacGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGt                                      <  2:19495/151‑1 (MQ=255)                                                                                         acGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTg                                     <  2:751653/151‑1 (MQ=255)                                                                                         acGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTg                                     >  2:503307/1‑151 (MQ=255)                                                                                           gACTTTAAACGTGCCAACGGTTACAGCGAAATGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGAGGTTCCAGTGCt                                   >  1:840949/1‑151 (MQ=255)                                                                                             cTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCttt                                 >  2:576717/1‑151 (MQ=255)                                                                                                 aaaCGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAAc                              <  1:582828/150‑1 (MQ=255)                                                                                                   aCGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACttt                           <  1:442381/151‑1 (MQ=255)                                                                                                      tGCCAACGGTTACAGCGAACTGGAGATCAGCAAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTgg                         <  1:743845/150‑1 (MQ=255)                                                                                                            cGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTTCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACtg                  <  2:149520/151‑1 (MQ=255)                                                                                                               ttACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACtggtg               >  2:748655/1‑151 (MQ=255)                                                                                                                  cAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACttt                           >  1:318843/1‑136 (MQ=255)                                                                                                                  cAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACttt                           <  2:318817/136‑1 (MQ=255)                                                                                                                  cAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCa            <  2:223244/151‑1 (MQ=255)                                                                                                                   aGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCAt           >  1:91778/1‑151 (MQ=255)                                                                                                                   aGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCAt           <  2:274972/151‑1 (MQ=255)                                                                                                                     cGAACTGGAGATCAGCAAGAAGCGCAGCACGCTGGAAACAGTGATGCTGACCGATTCCGTGGACACCCATAAAGCACGCCTGCATAACGCCGGGTTTGAGCATAGGGAGCTGTGGATCCAGTGCTTTAACTTTGGTTCACTGGTGGCAtt          >  2:247260/1‑150 (MQ=255)                                                                                                                        aCTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTaaa       >  2:606320/1‑150 (MQ=255)                                                                                                                         cTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTaaa       >  2:723754/1‑149 (MQ=255)                                                                                                                           ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGc    <  1:667514/150‑1 (MQ=255)                                                                                                                           ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCa   <  1:828707/151‑1 (MQ=255)                                                                                                                            gagaTCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCag  <  1:2244/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                 GGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAG  >  NC_000913/1953158‑1953431

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