Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,193,720	C→T	6.6%	R25C (CGC→TGC)	PP_5615 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,193,720	0	C	T	6.6%	176.7 / 3.0	60	R25C (CGC→TGC)	PP_5615	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/2);  total (30/30)
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.91e-01

GAACGCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAGGAGCGGCG  >  NC_002947/4193573‑4193862                                                                                                                                                    |                                                                                                                                               gAACGCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCg                                                                                                                                               >  1:128465/1‑149 (MQ=255)    cGCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGTAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCtt                                                                                                                                            <  5:327200/149‑1 (MQ=255)    cGCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCtt                                                                                                                                            <  7:322757/149‑1 (MQ=255)     gCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTg                                                                                                                                           <  7:241847/149‑1 (MQ=255)           cGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACt                                                                                                                                     >  4:107061/1‑149 (MQ=255)             tACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGc                                                                                                                                   >  8:91837/1‑149 (MQ=255)                gCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCgg                                                                                                                                >  1:371975/1‑149 (MQ=255)                  tCCGTCATGGAGGAAGTATGCTGAACGCTAAATTCAGAGGCTCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGa                                                                                                                              <  7:253898/149‑1 (MQ=255)                      tCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACATGCTTGGATACTGCCGGGAATCa                                                                                                                          >  8:32718/1‑149 (MQ=255)                      tCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGANGTATNCCGTTGCTCNACCTGCTTGGATACTGCCGGGAAt                                                                                                                            >  8:38399/1‑147 (MQ=255)                          ggAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTc                                                                                                                      >  4:9938/1‑149 (MQ=255)                          ggAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACCGCCGGGAATCAGTTc                                                                                                                      >  7:60105/1‑149 (MQ=255)                           gAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGtt                                                                                                                       <  2:371975/147‑1 (MQ=255)                           gAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGtt                                                                                                                       <  7:308660/147‑1 (MQ=255)                              gCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCa                                                                                                                  <  3:144747/149‑1 (MQ=255)                                      cTGAACGGGAAATTCAGAGGCGCTTGTAGAAGTGAGCGATGAAAATTACCGGAAGAGTCGAAGTCGAAACCGTTATAGTCGTCGTGTGCGATGTATGCCGTTGCTGCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACa                                                                                                           <  2:205271/148‑1 (MQ=255)                                       tGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAg                                                                                                         >  6:347963/1‑149 (MQ=255)                                            gcgAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCacac                                                                                                    >  7:230079/1‑149 (MQ=255)                                              gAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGc                                                                                                        >  6:213874/1‑143 (MQ=255)                                              gAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGc                                                                                                        <  5:213874/143‑1 (MQ=255)                                                aaTTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGGCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTggg                                                                                                <  8:204353/149‑1 (MQ=255)                                                     agagGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTAt                                                                                           <  2:128465/149‑1 (MQ=255)                                                      gagGCGCTTGTTGTAGTCAGCGTTGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACGCGCTTGGATTCTGCCGGGAATCAGTTGGGGACCCTACAAGCACACTGGGGCTAt                                                                                           <  3:266842/148‑1 (MQ=255)                                                        ggCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCTATGTATGCCGTTGCTCCACCTGCTTTGATACTTCCGGTAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGt                                                                                        <  7:32718/149‑1 (MQ=255)                                                        ggCGCNTGTAGAAGTCAGCGATGAAAATAANCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGNGCGATGTATGCCGTTGCTCCACCTGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGt                                                                                        <  7:38399/149‑1 (MQ=255)                                                         gcgcTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTa                                                                                       <  1:181939/149‑1 (MQ=255)                                                            cTTTTAGAAGTCATCGACGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGAGGTGGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTg                                                                                    <  8:173688/149‑1 (MQ=255)                                                                 aGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATg                                                                                          <  4:182027/138‑1 (MQ=255)                                                                 aGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATg                                                                                          >  3:182027/1‑138 (MQ=255)                                                                     gTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTc                                                                                                                      >  5:60190/1‑106 (MQ=255)                                                                     gTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTc                                                                                                                      <  6:60190/106‑1 (MQ=255)                                                                        aGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCGCACTGGGGCTATGGTACTGCCCATGATggg                                                                         >  4:7540/1‑148 (MQ=255)                                                                         gCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGa                                                                       <  8:230079/149‑1 (MQ=255)                                                                            aTGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAgcg                                                                    >  7:87580/1‑149 (MQ=255)                                                                                  aTAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTa                                                                 >  7:42982/1‑146 (MQ=255)                                                                                  aTAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTa                                                                 <  8:42982/146‑1 (MQ=255)                                                                                         ccagagTCGAACTCGAAACCGTTATAGACGTCGTGTGCGAAGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTCCAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCAc                                                       <  4:229576/147‑1 (MQ=255)                                                                                                 gaagtcgaaACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCtgtgtg                                                <  4:153557/148‑1 (MQ=255)                                                                                                     tcgaaACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGATCCACCCGATTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGAt                                           >  5:73656/1‑149 (MQ=255)                                                                                                       gaaACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATg                                                                                          >  1:213924/1‑100 (MQ=255)                                                                                                       gaaACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATg                                                                                          <  2:213924/100‑1 (MQ=255)                                                                                                        aaaCCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTgggg                                                                                               >  8:65775/1‑94 (MQ=255)                                                                                                        aaaCCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTgggg                                                                                               <  7:65775/94‑1 (MQ=255)                                                                                                        aaaCCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGc                                        >  3:337003/1‑149 (MQ=255)                                                                                                        aaaCCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGc                                        >  4:142022/1‑149 (MQ=255)                                                                                                         aaCCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGct                                       <  4:337003/149‑1 (MQ=255)                                                                                                             gTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGcttctt                                   <  7:91837/149‑1 (MQ=255)                                                                                                                   gACGTCGTCTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCa                                >  2:164310/1‑146 (MQ=255)                                                                                                                    aCGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACa                            >  8:208781/1‑149 (MQ=255)                                                                                                                      gtcgtGTTCGGTTTATGCCGTTGCTCCACCCGCTTGGATACATGCGGGAATGAGTTAGGCACCCTACAAACACACGGGGGCGTTGGTACTGCCCATGATGGGGAGCGCTATGAACATCACCGGTGTGAGGGTTTCTTCTTACAAACAtt                          >  2:102829/1‑149 (MQ=255)                                                                                                                       tcgtGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCTTGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTg                         <  1:102829/149‑1 (MQ=255)                                                                                                                       tcgtGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTg                         <  3:9938/149‑1 (MQ=255)                                                                                                                         gtgtGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTgg                        <  3:142022/148‑1 (MQ=255)                                                                                                                            tgCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTa                    >  7:305393/1‑149 (MQ=255)                                                                                                                                      gCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAg          >  4:252337/1‑149 (MQ=255)                                                                                                                                              tCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCGCCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAg             <  3:88737/138‑1 (MQ=255)                                                                                                                                              tCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCGCCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAg             >  4:88737/1‑138 (MQ=255)                                                                                                                                              tCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAGGAgcggc   >  4:249844/1‑148 (MQ=255)                                                                                                                                              tCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAGGAgcggc   >  8:356019/1‑148 (MQ=255)                                                                                                                                               ccACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAGGAgcggcg  >  1:149814/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                               GAACGCTTACCGTACGCTCCGTCCTGGAGGCAGTATGCTGAACGCGAAATTCAGAGGCGCTTGTAGAAGTCAGCGATGAAAATAACCGGAAGAGTCGAAGTCGAAACCGTTATAGACGTCGTGTGCGATGTATGCCGTTGCTCCACCCGCTTGGATACTGCCGGGAATCAGTTCGGCACCCTACAAGCACACTGGGGCTATGGTACTGCCCATGATGGGGAGCGCTATGAACTGCACCTGTGTGAGGATTGCTTCTTCCAAACATTGGCTTACCTCAAGCAGGAGCGGCG  >  NC_002947/4193573‑4193862

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