Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	46.3%	64.2 / 16.3	56	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base C (0/25);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 1.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868                                                                                                                                                      |                                                                                                                                                  ccGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTa                                                                                                                                                   >  1:880174/1‑150 (MQ=255)     gcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCttt                                                                                                                                              <  2:626804/151‑1 (MQ=255)     gcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCttt                                                                                                                                              <  1:416065/151‑1 (MQ=255)      cgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTc                                                                                                                                             <  2:548462/151‑1 (MQ=255)      cgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTc                                                                                                                                             <  2:455878/151‑1 (MQ=255)       gcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTcc                                                                                                                                            >  2:113794/1‑151 (MQ=255)               aGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacg                                                                                                                                     >  1:551721/1‑150 (MQ=255)                 tGGCTTCCGGCAACTGATAAATTTTCATCGAATAGACGCGACCACGGCGGGAGAAGCACAGAATATGGTCGTGGGTGTTCGCCCCCAGCAGTCGGTCGATAAAGCCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCacga                                                                                                                                    <  2:923668/149‑1 (MQ=255)                  ggtttCCGGCAACTAAAAAATTTTCATCAAATAGCCGCGCCCCCGGCGGGAAAAGCACAAAATATGGTCGGGAGTGTCCCCCACCAGCAGTCGGTCAATAAAGTCTTCTTTTTTAATCCGTGCGGCAGTTTTACCTTTCCCGCCacgacg                                                                                                                                  <  2:869511/147‑1 (MQ=255)                    cTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTg                                                                                                                               >  1:528997/1‑151 (MQ=255)                     ttCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTgc                                                                                                                              >  2:379667/1‑151 (MQ=255)                       ccGGCAACTGATAAACTTTCATCAAATAGACGCGCCCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTACCTTTCCCGCCACGACGCTgc                                                                                                                              <  1:544398/149‑1 (MQ=255)                         ggCAACTGAAAAACTTTCATCGAATAGACGCGCCCACGGCTGGAGAAGCACAGAATATGGTCGGGGGTGTCCCCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATACTTGGGGCAGTTTTCCCTTTCCCGCCACGACGCt                                                                                                                                <  2:873382/145‑1 (MQ=255)                            aaCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGt                                                                                                                       >  2:489412/1‑151 (MQ=255)                                   aacttttCATCGAATAGCCGCGCCCCGGGCTGGAGAAGCACAGAATATGGTCGTGTGTGTTCGCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTTTTTAATACGTGCGGCAAATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGa                                                                                                                <  2:677057/147‑1 (MQ=255)                                   aaaCTTTCATCGAATAGCCGCGCCCCGGGCGGGGAAAGCAAAAAATATGGTCGTGATTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTTTTCTTAAATCCGTGCGGCAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGa                                                                                                                <  2:643756/151‑1 (MQ=255)                                        ttCATCGAATAGACGCGCCCACGGCTGGAAAAGCACAGAATATGGTCGTGATTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAg                                                                                                             <  1:448185/149‑1 (MQ=255)                                                aTAGACGCGCCCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAc                                                                                                   <  1:542784/151‑1 (MQ=255)                                                aTAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAc                                                                                                   >  1:237154/1‑151 (MQ=255)                                                  aGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGGCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                 >  1:305938/1‑151 (MQ=255)                                                    aCGCGACCCCGGCTGGAGAAGCACAGAATAGGGTCGTGAGTTTTCGCCACCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATCCGTGCGGAAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACt                                                                                                  <  2:937332/148‑1 (MQ=255)                                                      gcgcCCCCGGCTGGAGAAGCACAGAATATGTTCGTGAGTTTTCGCCCCCAGCAGTCGGTGGATAAAGTCTTCTTTTTTAATTCGGGCGGCATATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATGCt                                                                                                  <  2:631053/146‑1 (MQ=255)                                                       cgcccACGGCTGGAGAAGCACAGAATATGGTCGTGAGTTTCCGCCCCCAGCAGTCGTTCAATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                              <  2:74957/146‑1 (MQ=255)                                                        gACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            >  2:425317/1‑150 (MQ=255)                                                          ccACGGCGGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  1:580300/148‑1 (MQ=255)                                                              ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCCCCACCAGCAGGCGGTCGATAAAGTTTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg                                                                                     <  1:316699/151‑1 (MQ=255)                                                                cTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                   >  1:525371/1‑151 (MQ=255)                                                                 tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  >  1:644707/1‑151 (MQ=255)                                                                  ggAGAAGCAAAAAATTTGGTTGTGGGTGTTCGCCACCAGCGGTCGGTCAATAAAGTTTTCTTTTTAAATCCGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                   <  1:407602/149‑1 (MQ=255)                                                                      aaccacaGAATATGGTCGTGAGTGTCCGCCGCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATGCGTGCGGAAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                   <  2:395312/142‑1 (MQ=255)                                                                      aaGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                                >  1:660254/1‑148 (MQ=255)                                                                      aaGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                                >  2:280532/1‑148 (MQ=255)                                                                      aaGCACAAAATATGGCCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTATTCCGTGGGGCAGATTTTCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgaga                                                                               <  2:218554/149‑1 (MQ=255)                                                                         caaaGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTAAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                            <  1:432445/146‑1 (MQ=255)                                                                          acaAAATATGGTCGTGAGTGTTCGCTCCCAGCGGTCGGTCGAAAAAGTTTTCTTTTTTAATCCGGGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                            <  2:857987/148‑1 (MQ=255)                                                                           caGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTc                                                                        >  1:118468/1‑151 (MQ=255)                                                                                      cGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc                                                             >  2:723962/1‑151 (MQ=255)                                                                                        tGAGTGTCCCCCCCCAGCATTCGGTCGATAAATTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACAt                                                              <  2:643257/148‑1 (MQ=255)                                                                                            tgtTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                       >  2:269659/1‑151 (MQ=255)                                                                                              ttCGCCCCCAGCAGTCGGTCAAAAAAGTCTTCTTTTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                       <  1:786367/149‑1 (MQ=255)                                                                                              tcggCCACCAGCGGTGGGTCAATAAGTTCTTTTTCTTTACTAGGTGGGGAAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgg                                                      <  2:800120/147‑1 (MQ=255)                                                                                              tccGCCCCAAGCAGTTGGTCGATAAACTTTTTTTTTTAAAGGCGTGCGAAAATTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                       <  1:561749/147‑1 (MQ=255)                                                                                                cGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCATATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                   <  2:269386/151‑1 (MQ=255)                                                                                                 gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                   >  2:127848/1‑150 (MQ=255)                                                                                                  ccaccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  1:795012/147‑1 (MQ=255)                                                                                                              ggtcgATAAAGTCTTCTTTTTTAATACGTGCGGCAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                       <  2:716255/149‑1 (MQ=255)                                                                                                                 cgAAAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGt                                   <  2:613503/150‑1 (MQ=255)                                                                                                                  gaaaaaGTCTTCTTCTTTAATACGTGCGGCGGTTTTTCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                    <  1:562591/145‑1 (MQ=255)                                                                                                                     aaaTTCTTCTTCTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt                               <  1:904530/150‑1 (MQ=255)                                                                                                                     aaaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                              <  1:179626/151‑1 (MQ=255)                                                                                                                           ttcttcttTAATCCGTGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt                           <  2:160585/148‑1 (MQ=255)                                                                                                                                 ttttatccgtcaggCAGATTTCCTTTTTCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                      <  1:50951/135‑1 (MQ=255)                                                                                                                                   tAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc                >  2:345143/1‑151 (MQ=255)                                                                                                                                         gggCGGCATTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTg            <  2:830064/147‑1 (MQ=255)                                                                                                                                         gggCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGt             <  2:551656/146‑1 (MQ=255)                                                                                                                                                aGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt    <  1:511772/150‑1 (MQ=255)                                                                                                                                                  aTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac  <  1:526509/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868

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