Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,039,297	A→C	8.8%	intergenic (‑261/+33)	PP_5290 ← / ← pyrE	transposase/orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,039,297	0	A	C	8.8%	157.2 / 6.1	56	intergenic (‑261/+33)	PP_5290/pyrE	transposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base C (3/2);  total (29/27)
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.86e-01

CTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGA  >  NC_002947/6039168‑6039445                                                                                                                                   |                                                                                                                                                     cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAtt                                                                                                                                  <  1:447329/150‑1 (MQ=255)      gTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                             <  1:417996/150‑1 (MQ=255)      gTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTAAACGCGAAAGGGCAAGATTAGGGGCGAG‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                             <  1:692568/150‑1 (MQ=18)          gTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGc                                                                                                                                            <  2:451569/131‑1 (MQ=18)          gTACTGCGCGAACAGTCAGAAGCGAGGTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGc                                                                                                                                            >  1:451569/1‑131 (MQ=18)              tGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                     >  2:945892/1‑150 (MQ=255)                cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATc                                                                                                                    <  2:857623/150‑1 (MQ=255)                  cgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAAGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCATGAGCCAGGATCAGa                                                                                                                 >  1:334783/1‑150 (MQ=255)                         taaGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTa                                                                                                          >  2:1556740/3‑150 (MQ=255)                                 gAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGccc                                                                                                  >  1:439451/1‑150 (MQ=255)                                   gTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCgg                                                                                                >  2:333057/1‑150 (MQ=255)                                      aTGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAg                                                                                             <  1:385033/150‑1 (MQ=255)                                          cccTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                         <  1:1002488/150‑1 (MQ=255)                                           cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTc                                                                                        >  1:1543299/1‑150 (MQ=255)                                           cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTc                                                                                        >  2:979118/1‑150 (MQ=255)                                             tcctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                                      >  2:460856/1‑150 (MQ=255)                                                  ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTg                                                                                                                              <  2:1047875/105‑1 (MQ=255)                                                  ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTg                                                                                                                              >  1:1047875/1‑105 (MQ=255)                                                      cTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTGTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCGCGGTAGGCCTCGACAGCCGGCa                                                                             >  1:21826/1‑150 (MQ=255)                                                           aGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCCGAAACGGCGCTTG‑AAAGGAAAGCGAGCTTCTGAATAAGGTCGCCTTTTTCATTCCTGGGCAAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGc                                                                        >  1:31022/1‑150 (MQ=255)                                                           aGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGc                                                                        >  1:1220626/1‑150 (MQ=255)                                                             tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                   <  1:1007732/89‑1 (MQ=255)                                                             tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                   >  2:1007732/1‑89 (MQ=255)                                                              ccTACCATACAAGCGGGTTAAGCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                                     <  2:1613712/150‑1 (MQ=255)                                                                tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGaagaag                                                                                                                                                 >  2:15128/1‑72 (MQ=255)                                                                tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGaagaag                                                                                                                                                 <  1:15128/72‑1 (MQ=255)                                                                 aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                      <  1:480659/82‑1 (MQ=255)                                                                 aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                      >  2:480659/1‑82 (MQ=255)                                                                        aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgtg                                                           >  1:1540781/1‑150 (MQ=255)                                                                         aGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                        >  1:274440/1‑72 (MQ=21)                                                                         aGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                        <  2:274440/72‑1 (MQ=25)                                                                           cGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGAt                                                        <  1:945892/150‑1 (MQ=255)                                                                            gggTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                       <  2:850438/150‑1 (MQ=255)                                                                                  aCCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                 <  1:460856/150‑1 (MQ=255)                                                                                   cccGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGcc                                                >  2:46447/1‑150 (MQ=255)                                                                                      gcgAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAgg                                             >  2:45433/1‑150 (MQ=255)                                                                                      gcgAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAgg                                             >  2:1339557/1‑150 (MQ=255)                                                                                        gAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGaa                                           <  2:1543299/150‑1 (MQ=255)                                                                                        gAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGaa                                           <  1:358521/150‑1 (MQ=255)                                                                                           aGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTAATGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCGGTTTCAGCTGTGGATCGTCAGGCAGGAACTc                                        <  2:495349/150‑1 (MQ=255)                                                                                              gCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                     >  2:1536347/1‑150 (MQ=255)                                                                                                     cAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                              <  1:333057/150‑1 (MQ=255)                                                                                                           gcAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGAATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGc                        <  2:1220626/150‑1 (MQ=255)                                                                                                             gaatgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCg                       >  2:1355422/6‑150 (MQ=38)                                                                                                             ac‑aAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGa                      >  1:643444/1‑150 (MQ=255)                                                                                                             ac‑aAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTACAGCCCCTGGATcggaga                      >  2:1168859/1‑145 (MQ=255)                                                                                                                aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                    <  1:1556740/150‑1 (MQ=255)                                                                                                                aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                    <  1:979118/150‑1 (MQ=255)                                                                                                                     aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg               <  2:111641/150‑1 (MQ=255)                                                                                                                     aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg               <  2:439451/150‑1 (MQ=255)                                                                                                                       ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTTGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                                 >  1:715359/1‑98 (MQ=255)                                                                                                                       ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                                 <  2:715359/98‑1 (MQ=255)                                                                                                                       ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt             <  1:1536347/150‑1 (MQ=255)                                                                                                                       ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt             >  1:782761/1‑150 (MQ=255)                                                                                                                              ctctGAAGAAGGTCTCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc              <  2:1326883/142‑1 (MQ=255)                                                                                                                              ctctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc              >  1:1326883/1‑142 (MQ=255)                                                                                                                               tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg     >  1:193256/1‑150 (MQ=255)                                                                                                                               tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg     >  1:1500451/1‑150 (MQ=255)                                                                                                                                  gaagaaGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa  <  1:1339557/150‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                     CTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGA  >  NC_002947/6039168‑6039445

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