Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,812,778	C→G	7.0%	intergenic (+72/‑38)	rplT → / → pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,812,778	0	C	G	7.0%	177.0 / 3.8	57	intergenic (+72/‑38)	rplT/pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base G (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGAG  >  NC_002947/2812635‑2812926                                                                                                                                                 |                                                                                                                                                     gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCaa                                                                                                                                                    <  7:334485/146‑1 (MQ=255)    gATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGctct                                                                                                                                              >  8:263128/1‑149 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                                                         <  3:107115/148‑1 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTc                                                                                                                                            <  5:259651/145‑1 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCAATAGGGGAAGAGCTTTGAAGGCTCTTc                                                                                                                                            <  3:213467/145‑1 (MQ=255)                aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                  <  7:224593/149‑1 (MQ=255)                aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                  >  5:171032/1‑149 (MQ=255)                   gAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCATTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTa                                                                                                                               >  2:251968/1‑149 (MQ=255)                     aaaaGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATc                                                                                                                             >  6:190072/1‑149 (MQ=255)                      aaaGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                                                            >  8:90341/1‑149 (MQ=255)                               tttGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTg                                                                                                                    >  8:10841/1‑148 (MQ=255)                               tttGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGt                                                                                                                   <  8:211115/149‑1 (MQ=255)                                  gctgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACa                                                                                                                >  7:98774/1‑149 (MQ=255)                                   ctgcGATTGTCGAGAAAGCTAAAGCCATTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACAt                                                                                                               <  1:251968/149‑1 (MQ=255)                                            tCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt                                                                                                      >  7:187413/1‑149 (MQ=255)                                             cGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTg                                                                                                     >  8:130407/1‑149 (MQ=255)                                                 aaaGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATg                                                                                                 >  1:333689/1‑149 (MQ=255)                                                    gCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATgcgc                                                                                              >  5:222508/1‑149 (MQ=255)                                                      tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACTCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCt                                                                                             >  1:61596/1‑149 (MQ=255)                                                         aGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGAt                                                                                                  <  4:161831/140‑1 (MQ=255)                                                         aGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGAt                                                                                                  >  3:161831/1‑140 (MQ=255)                                                           ccGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc                                                                                       <  7:10841/149‑1 (MQ=255)                                                             gTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTccc                                                                                     <  7:90341/149‑1 (MQ=255)                                                              ttCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCa                                                                                    <  7:263128/149‑1 (MQ=255)                                                                     ttAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATTGAAAACCTGGATGCGCTGGTCTCCCAAGccc                                                                               >  6:126258/1‑147 (MQ=255)                                                                      tAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCACTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGccc                                                                               >  4:87919/1‑146 (MQ=255)                                                                        aGTAGCCGCGCCAATGAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                                                                           <  4:220911/148‑1 (MQ=255)                                                                        aGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                                                                           <  6:222508/148‑1 (MQ=255)                                                                            cccGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGTAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGcc                                                                       <  5:49176/148‑1 (MQ=255)                                                                             ccGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGt                                                                     >  3:305079/1‑149 (MQ=255)                                                                              cgcgACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCa                                                                                                                                                     >  7:14816/1‑68 (MQ=255)                                                                              cgcgACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCa                                                                                                                                                     <  8:14816/68‑1 (MQ=255)                                                                                           ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                        <  8:187413/149‑1 (MQ=255)                                                                                           ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa                                                        <  8:98774/149‑1 (MQ=255)                                                                                                       gcgcTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGcc                                            >  7:295268/1‑149 (MQ=255)                                                                                                             ccGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTAc                                                                                                                 >  6:16858/1‑74 (MQ=255)                                                                                                             ccGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTAc                                                                                                                 <  5:16858/74‑1 (MQ=255)                                                                                                                   aTTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATc                                >  6:42412/1‑149 (MQ=255)                                                                                                                      gAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCgg                             >  8:175218/1‑149 (MQ=255)                                                                                                                       aaCGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCgg                             <  2:242654/148‑1 (MQ=255)                                                                                                                       aaCGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCgg                             <  7:130407/148‑1 (MQ=255)                                                                                                                             aTCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAtt                      <  8:358450/149‑1 (MQ=255)                                                                                                                             aTCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAtt                      <  6:171032/149‑1 (MQ=255)                                                                                                                                gATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtc                   >  7:340673/1‑149 (MQ=255)                                                                                                                                  tAGGGGAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa                        >  1:330576/1‑142 (MQ=255)                                                                                                                                  tAGGGGAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa                        <  2:330576/142‑1 (MQ=255)                                                                                                                                   aGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCg                >  4:274387/1‑149 (MQ=255)                                                                                                                                     gggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCgg               >  4:88354/1‑148 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCt                                          <  8:24572/121‑1 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCt                                          >  7:30184/1‑121 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCt                                          <  8:30184/121‑1 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCt                                          >  7:24572/1‑121 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAt                       <  8:204356/140‑1 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAAt                       >  7:204356/1‑140 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTACGTATCTGGAGTCTGTACATGGAAAACATGGATGCGCTGTTATCCCAACCCCTAGAGGCCGTGCACCGCCCGCAAGACAACAACGCCCTGTAACAGATCCGGGATAATAATCACgg               >  2:160871/1‑148 (MQ=255)                                                                                                                                              cTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGc     <  5:126258/149‑1 (MQ=255)                                                                                                                                              cTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGc     <  5:190072/149‑1 (MQ=255)                                                                                                                                                 cAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGAg  >  4:53975/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCAC‑CCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGAG  >  NC_002947/2812635‑2812926

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: ‑