Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	924,324:1	+A	6.8%	coding (872/2169 nt)	PP_0804 ←	protein secretion ABC transporter permease/ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	924,324	1	.	A	6.8%	229.8 / 10.8	59	coding (872/2169 nt)	PP_0804	protein secretion ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (29/26);  new base A (2/2);  total (31/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 = 7.82e-01

ACCGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCGG  >  NC_002947/924181‑924462                                                                                                                                                  |                                                                                                                                           aCCGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGc                                                                                                                                      <  2:963399/150‑1 (MQ=255)   cGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTg                                                                                                                                    >  2:1086455/1‑150 (MQ=255)     tGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACATCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGa                                                                                                                                  <  1:280843/150‑1 (MQ=255)        cAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCGAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGt                                                                                                                               <  1:242292/150‑1 (MQ=255)        cAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGt                                                                                                                               >  1:156072/1‑150 (MQ=255)         aGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTg                                                                                                                              >  1:7463/1‑150 (MQ=255)           gcgcCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTAAGACGCGGAGCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGAt                                                                                                                            >  2:947626/1‑150 (MQ=255)               cTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAgg                                                                                                                        >  1:25455/1‑150 (MQ=255)                     gcaagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCgg                                                                                                                  >  1:489874/1‑150 (MQ=255)                      caagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGa                                                                                                                 <  1:303412/150‑1 (MQ=255)                       aagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGTCCAGCACCGCGA‑CGGTGCTggcc                                                                                                                                 >  1:1097800/1‑131 (MQ=255)                       aagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGTCCAGCACCG‑CGACGGTGCTggcc                                                                                                                                 <  2:1097800/133‑3 (MQ=255)                            aCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACt                                                                                                           >  2:729179/1‑150 (MQ=255)                                   cagcaCCATGGCCAGCAGCAAGACCTCCGGAATGAACAGCAGTACGCCGCCGATCAGCCAGAACAACGAACA‑GAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTc                                                                                                    <  1:209455/150‑1 (MQ=255)                                    agcaCCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTcg                                                                                                   >  2:445248/1‑150 (MQ=255)                                       aCCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCgcgc                                                                                                <  2:25455/150‑1 (MQ=255)                                         cATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCaa                                                                                              >  1:637403/1‑150 (MQ=255)                                          aTGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAAc                                                                                             <  2:637403/150‑1 (MQ=255)                                               cagcagCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCga                                                                                        >  1:949909/1‑150 (MQ=255)                                                  cagcGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGAt                                                                                     <  1:1111209/150‑1 (MQ=255)                                                     cgcgACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGaa                                                                                  <  1:877582/150‑1 (MQ=255)                                                            aGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCg                                                                           >  1:238864/1‑150 (MQ=255)                                                             gCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAAc                                                                                 <  2:1052680/143‑1 (MQ=255)                                                             gCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGAAGGCAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAAc                                                                                 >  1:1052680/1‑143 (MQ=255)                                                              cGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGt                                                                         <  1:729179/150‑1 (MQ=255)                                                              cGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGt                                                                         >  1:242313/1‑150 (MQ=255)                                                                ggATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTgg                                                                       >  1:111459/1‑150 (MQ=255)                                                                 gATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGa                                                                      <  1:352407/150‑1 (MQ=255)                                                                  aTGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGAc                                                                     <  1:356296/150‑1 (MQ=255)                                                                      acacCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTg                                                                                                                                    >  1:610198/1‑83 (MQ=255)                                                                      acacCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTg                                                                                                                                    <  2:610198/83‑1 (MQ=255)                                                                      acacCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTgg                                                                 <  2:292492/150‑1 (MQ=255)                                                                         ccAGTACGCCGCCGATCAGCCACAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCTAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCCGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGcc                                                              <  1:360404/150‑1 (MQ=255)                                                                             tACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCAc                                                          >  2:539243/1‑150 (MQ=255)                                                                                gccgccGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGa                                                       <  1:539243/150‑1 (MQ=255)                                                                                           gCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcaggc                                            >  1:238392/1‑150 (MQ=255)                                                                                               gAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcaggcgcag                                        <  2:242313/150‑1 (MQ=255)                                                                                                aaC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGcgcaggcgcagg                                       <  2:949909/150‑1 (MQ=255)                                                                                                            gaacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTggg                                                                <  2:1079175/114‑1 (MQ=255)                                                                                                            gaacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTggg                                                                >  1:1079175/1‑114 (MQ=255)                                                                                                             aacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAg                           <  2:156072/150‑1 (MQ=255)                                                                                                                   agGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAgg                     >  2:158198/1‑150 (MQ=255)                                                                                                                        gaaCGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCt                >  2:406246/1‑150 (MQ=255)                                                                                                                           cGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAc             >  1:1007144/1‑150 (MQ=255)                                                                                                                            ggCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAAAGAGCCGGTGGa                                                                      <  2:697691/92‑1 (MQ=255)                                                                                                                            ggCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAAAGAGCCGGTGGa                                                                      >  1:697691/1‑92 (MQ=255)                                                                                                                                  tCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGt                                                                         >  2:126389/1‑83 (MQ=255)                                                                                                                                  tCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGt                                                                         <  1:126389/83‑1 (MQ=255)                                                                                                                                  tCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCa                       <  2:562505/133‑1 (MQ=255)                                                                                                                                  tCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCa                       >  1:562505/1‑133 (MQ=255)                                                                                                                                    gtccAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGc                                     >  1:833090/3‑117 (MQ=255)                                                                                                                                    gtccAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGc                                     <  2:833090/115‑1 (MQ=255)                                                                                                                                    gtccAGCACCGCGA‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTc    >  2:1140587/3‑150 (MQ=255)                                                                                                                                    gtccAGCACCG‑CGACGGTGCTGGCCGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcag                                              >  1:942782/3‑108 (MQ=255)                                                                                                                                    gtccAGCACCG‑CGACGGTGCTGGCCGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcag                                              <  2:942782/106‑1 (MQ=255)                                                                                                                                     gCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCg   >  2:623078/1‑150 (MQ=255)                                                                                                                                      ccAGCACCG‑CGACGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCgg  >  1:550967/1‑150 (MQ=255)                                                                                                                                           aCCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAccc           >  2:1120131/1‑136 (MQ=255)                                                                                                                                           aCCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAccc           <  1:1120131/136‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           ACCGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAAC‑ACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCG‑CGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCGG  >  NC_002947/924181‑924462

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