Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	390,782	C→T	6.9%	V352V (GTC→GTT)	soxA →	sarcosine oxidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	390,782	0	C	T	6.9%	170.9 / 2.4	58	V352V (GTC→GTT)	soxA	sarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGGGCAACATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGA  >  NC_002947/390636‑390921                                                                                                                                                    |                                                                                                                                               gTGCCGGGCAACATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGc                                                                                                                                              <  2:94108/149‑1 (MQ=255)         cAACATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGcc                                                                                                                                       >  1:250280/1‑148 (MQ=255)          aaCATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGa                                                                                                                                     >  3:66168/1‑149 (MQ=255)          aaCATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGa                                                                                                                                     >  4:164202/1‑149 (MQ=255)           aCATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGAc                                                                                                                                    <  7:94020/149‑1 (MQ=255)               gCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGAcgccc                                                                                                                                >  7:209728/1‑149 (MQ=255)                   gCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGAcgcccgc                                                                                                                              <  4:212826/147‑1 (MQ=255)                   gCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACgccagcca                                                                                                                            >  2:185116/1‑149 (MQ=255)                          cTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAAccc                                                                                                                      <  4:107111/148‑1 (MQ=255)                          cTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCa                                                                                                                     <  8:77184/149‑1 (MQ=255)                          cTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCa                                                                                                                     >  7:77184/1‑149 (MQ=255)                           tGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCAAGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCa                                                                                                                     >  8:165883/1‑148 (MQ=255)                             tATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTACGGGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACg                                                                                                                   >  1:66915/1‑149 (MQ=255)                                ccACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGc                                                                                                               >  2:128868/1‑149 (MQ=255)                                  aCCTATGTTCGCCGCTATGGCGTGGCGCCGGGACGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTc                                                                                                             >  1:169363/1‑149 (MQ=255)                                    cTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGc                                                                                                           <  1:128868/149‑1 (MQ=255)                                      aTGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTg                                                                                                         >  8:210299/1‑149 (MQ=255)                                         ttCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGAGCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACaa                                                                                                                         >  1:190411/1‑130 (MQ=255)                                         ttCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACaa                                                                                                                         <  2:190411/130‑1 (MQ=255)                                            gccgcTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTgaa                                                                                                   <  4:66168/149‑1 (MQ=255)                                                   tggcgtggcgCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAgcgc                                                                                            <  3:164202/149‑1 (MQ=255)                                                     gcgtggcgCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGcgcg                                                                                           <  8:88998/148‑1 (MQ=255)                                                        tggcgCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGc                                                                                        >  2:197206/1‑148 (MQ=255)                                                          gcgcCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCaa                                                                                     >  7:61843/1‑149 (MQ=255)                                                           cgcCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCaaa                                                                                    <  1:232190/149‑1 (MQ=255)                                                                  tCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTc                                                                             <  7:210299/149‑1 (MQ=255)                                                                  tCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTc                                                                             <  2:250280/149‑1 (MQ=255)                                                                    tcAATTTGGTAC‑TGTCGAACAACAACGAACACGCTTATCGCGCAGCGCTGGACTGGCACGACGAAGGCCTGCAAGTGGTTGACATCGACGACGACCGACACAAACCACGTGGCTCGCTGGTTGAAGAATCGCGTGCCAAAGGCATTCg                                                                            >  5:61906/2‑148 (MQ=255)                                                                    gCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTTGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCg                                                                            >  5:65688/1‑148 (MQ=255)                                                                     cAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTgg                                                                                                        <  1:146507/119‑1 (MQ=255)                                                                     cAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTgg                                                                                                        >  2:146507/1‑119 (MQ=255)                                                                     cAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCa                                                                          >  8:17060/1‑149 (MQ=255)                                                                      aaGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCa                                                                          >  2:179157/1‑148 (MQ=255)                                                                       aGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATa                                                                        <  1:255020/149‑1 (MQ=255)                                                                         ttGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACACCCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACt                                                                      >  4:186409/1‑149 (MQ=255)                                                                           ggTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTc                                                                     <  4:107103/148‑1 (MQ=255)                                                                           ggTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTc                                                                     >  3:107103/1‑148 (MQ=255)                                                                                 tGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTc                                                                >  6:203290/1‑148 (MQ=255)                                                                                   tCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCa                                                              <  8:209728/148‑1 (MQ=255)                                                                                        caacaaCGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCg                                                     >  3:187054/1‑149 (MQ=255)                                                                                          acaacGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGt                                                                                                                  <  5:93928/89‑1 (MQ=255)                                                                                          acaacGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGt                                                                                                                  >  6:93928/1‑89 (MQ=255)                                                                                          acaacGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTg                                                   <  8:61843/149‑1 (MQ=255)                                                                                              cGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATc                                                >  2:54594/1‑148 (MQ=255)                                                                                                aCCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAg                                             <  5:203290/149‑1 (MQ=255)                                                                                                aCCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAg                                             <  2:230796/149‑1 (MQ=255)                                                                                                 ccACGCTTATCGCGCCGCGCTGGACTTGCACGACGCAGGCCTGCAAGTGGTTGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAgg                                            <  6:65688/149‑1 (MQ=255)                                                                                                 ccACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAgg                                            <  2:169363/149‑1 (MQ=255)                                                                                                 acACGCTTATCGGGACGCTCTGGACTGTCAAGACGCAGGCCTGAAAGTGGTTGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTTATCGAgg                                            <  6:61906/148‑1 (MQ=255)                                                                                                      cTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCaaa                                       <  7:165883/149‑1 (MQ=255)                                                                                                             cgccgcgcTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTgg                                                                                                                >  4:131987/1‑72 (MQ=255)                                                                                                                  cgcTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCAc                           <  7:17060/149‑1 (MQ=255)                                                                                                                          tGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCgg                   >  1:95566/1‑149 (MQ=255)                                                                                                                               cgacgCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGccc              <  1:54594/149‑1 (MQ=255)                                                                                                                                   gCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCATGGGGTGCGTGATCGAGGCCAAAGGCAGCAAGCACTTCACCGGCGCCCg             >  1:89868/1‑149 (MQ=255)                                                                                                                                   gCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCgtgt          <  2:66915/149‑1 (MQ=255)                                                                                                                                           gCAAGTTGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGcga  >  8:149220/1‑149 (MQ=255)                                                                                                                                            cAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGcga  <  1:197206/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GTGCCGGGCAACATGCTGGCCGGTGCTGTATCCACCTATGTTCGCCGCTATGGCGTGGCGCCGGGTCGCAAGTTGGTAC‑TGTCGACCAACAACGACCACGCTTATCGCGCCGCGCTGGACTGGCACGACGCAGGCCTGCAAGTGGTCGCCATCGCCGACGCCCGCCACAACCCACGTGGCTCGCTGGTTGAAGAAGCGCGTGCCAAAGGCATTCGCATACTCACCTCCA‑‑GCG‑CCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGA  >  NC_002947/390636‑390921

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