Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,866,697	A→C	12.0%	intergenic (+20/+37)	PP_5140 → / ← thyA	MerR family transcriptional regulator/thymidylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,866,697	0	A	C	12.0%	155.1 / 10.3	57	intergenic (+20/+37)	PP_5140/thyA	MerR family transcriptional regulator/thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (4/3);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGGC  >  NC_002947/5866551‑5866844                                                                                                                                                   |                                                                                                                                                     ggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAAc                                                                                                                                                    <  5:217703/149‑1 (MQ=255)          gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                           <  1:197763/149‑1 (MQ=255)             gcgcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGTGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt                                                                                                                                       <  7:212245/149‑1 (MQ=255)               gcTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCa                                                                                                                                     <  1:28376/149‑1 (MQ=255)                 tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgcg                                                                                                                                               >  6:178413/1‑138 (MQ=255)                 tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgcg                                                                                                                                               <  5:178413/138‑1 (MQ=255)                 tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt                                                                                                                                       >  7:245417/1‑145 (MQ=255)                 tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt                                                                                                                                       <  8:245417/145‑1 (MQ=255)                           aGGTGAGTGAATTCGCGTTCTTGCTGCAACTGGAGACACAGGGGGCGGTTTCGGTACCGGATACCGAACATTCCGATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGTTGAACCCGCGAAGA‑GGTCAGTTCTAACAATa                                                                                                                         <  6:11475/149‑1 (MQ=255)                                         gcgATCTTGCAGCAACTGGAGCCAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGGTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGAc                                                                                                           <  3:35037/149‑1 (MQ=255)                                            aTCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGc                                                                                                        >  7:67311/1‑149 (MQ=255)                                                 gcagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCg                                                                                                   >  2:58059/1‑149 (MQ=255)                                                  cagcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCgg                                                                                                  <  2:196181/149‑1 (MQ=255)                                                   agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt                                                                                                                                       >  3:208619/1‑111 (MQ=255)                                                   agcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGt                                                                                                                                       <  4:208619/111‑1 (MQ=255)                                                      aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACTGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt                                                                                              <  3:160639/149‑1 (MQ=255)                                                      aaCTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt                                                                                              >  7:46203/1‑149 (MQ=255)                                                         tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa                                                                                                                                             <  7:236445/100‑1 (MQ=255)                                                         tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa                                                                                                                                             >  8:236445/1‑100 (MQ=255)                                                         tGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGt                                                                                            >  5:37907/1‑148 (MQ=255)                                                                acaCGGGGCTGTCTCTGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCAGTGAGCCAAAAGCTTAGCGGGTGACCCCGCGAAGA‑GGTCAGTTCTAACAATATATATCTATCAGACCGCCATCGGCGCTGTCATCGgcg                                                                                    <  6:149451/147‑1 (MQ=255)                                                                aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAAGAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATAGgcg                                                                                    >  8:96996/1‑149 (MQ=255)                                                                aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACAGCCATCGGCGCTGTCATCGgcg                                                                                    >  5:227818/1‑149 (MQ=255)                                                                       gCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTg                                                                             <  6:37907/149‑1 (MQ=255)                                                                           tATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCa                                                                         <  8:46203/149‑1 (MQ=255)                                                                                 tACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACTGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCt                                                                   <  3:235935/149‑1 (MQ=255)                                                                                     ggATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca                                                               <  8:233054/149‑1 (MQ=255)                                                                                         aCCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGcc                                                                    <  5:198224/140‑1 (MQ=255)                                                                                         aCCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGcc                                                                    >  6:198224/1‑140 (MQ=255)                                                                                           cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagccag                                                          >  6:49558/1‑148 (MQ=255)                                                                                           cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagccag                                                          >  5:136548/1‑148 (MQ=255)                                                                                            gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGTCCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaa                                                        <  1:70997/149‑1 (MQ=255)                                                                                            gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCAAACAATa                                                                                                                         <  5:205999/84‑1 (MQ=255)                                                                                            gAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCAAACAATa                                                                                                                         >  6:205999/1‑84 (MQ=255)                                                                                             aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGaaa                                                       <  6:227818/149‑1 (MQ=255)                                                                                                aTTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCa                                                    >  7:241913/1‑149 (MQ=255)                                                                                                  tCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACt                                                  >  6:159922/1‑149 (MQ=255)                                                                                                   cccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTc                                                 >  7:30685/1‑149 (MQ=255)                                                                                                    ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTc                                                 <  8:241913/148‑1 (MQ=255)                                                                                                    ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg                                                >  6:268046/1‑149 (MQ=255)                                                                                                    ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg                                                >  8:203856/1‑149 (MQ=255)                                                                                                    ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTCCTCATAGCCTTCTAGCCAGAAATCACTCg                                                <  4:179712/149‑1 (MQ=255)                                                                                                    ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTACTCATAGCCTTCTAGCCAGAAATCACTCg                                                >  3:179712/1‑149 (MQ=255)                                                                                                     cATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCgg                                               >  2:84889/1‑149 (MQ=255)                                                                                                                     aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGc                               >  6:61562/1‑149 (MQ=255)                                                                                                                          cATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCATCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTc                          >  7:142493/1‑149 (MQ=255)                                                                                                                                 ccAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa                   >  6:8374/1‑149 (MQ=255)                                                                                                                                 ccAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa                   <  8:30685/149‑1 (MQ=255)                                                                                                                                   aaaaGCTTCGCGGGTGAACCCGCGAAGAGGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCa                   >  1:152875/1‑148 (MQ=255)                                                                                                                                      ctcTTCGCGGGTTCACCCGCGAAGC‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACtt              >  1:157266/3‑149 (MQ=255)                                                                                                                                        cTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc            <  3:148897/149‑1 (MQ=255)                                                                                                                                        cTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc            <  7:203856/149‑1 (MQ=255)                                                                                                                                         ttCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc            >  8:146764/1‑148 (MQ=255)                                                                                                                                          tCGCGGGTTCACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc            >  2:196323/1‑147 (MQ=255)                                                                                                                                           cgcgGGTTCACCCGCGAAGC‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa          <  2:157266/148‑1 (MQ=255)                                                                                                                                           cgcgGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa          <  5:268046/148‑1 (MQ=255)                                                                                                                                                gTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCg    <  5:49558/149‑1 (MQ=255)                                                                                                                                                   aaCCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGGc  <  7:146764/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGA‑GGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGGC  >  NC_002947/5866551‑5866844

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