Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,494	0	T	G	53.7%	56.1 / 84.2	54	intergenic (‑512/+91)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (3/22);  new base G (29/0);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 2.95e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-09

TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTA  >  NC_000913/987349‑987634                                                                                                                                                  |                                                                                                                                                  ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                                                                                   >  1:622694/1‑146 (MQ=255) ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                                                                                               <  1:306766/150‑1 (MQ=255)    aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                                                                                           <  1:918740/151‑1 (MQ=255)          caaGAATGCCAAAAAATAGGATGAAACTTGTGCTATTAGCTAAGCAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATAGCGGAGTTgggaaaaaagt                                                                                                                                       >  1:528509/3‑139 (MQ=255)           cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAA‑TT‑gtt                                                                                                                                  >  2:500389/1‑149 (MQ=255)           cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAA‑TT‑GCt                                                                                                                                  >  2:341097/1‑151 (MQ=255)            aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCtt                                                                                                                                 <  1:245363/151‑1 (MQ=255)             gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTg                                                                                                                                <  1:465159/151‑1 (MQ=255)               aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAA‑TT‑GCtt                                                                                                                                 >  1:767475/1‑148 (MQ=255)                   cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑a                                                                                                                         >  2:444476/1‑150 (MQ=255)                     aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGA‑AT‑GCTTG‑ACtt                                                                                                                           >  2:567319/1‑147 (MQ=255)                        aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTCAAAAAAA‑TT‑GTTTT‑ACTT‑aaaaaa                                                                                                                    >  2:411598/1‑150 (MQ=255)                        aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGAGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑aacaaa                                                                                                                    >  1:539187/1‑150 (MQ=255)                                   aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGa                                                                                                        <  2:433343/151‑1 (MQ=255)                                   aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGA‑TT‑TCTTG‑TCTT‑TACAAAAAAATAAAtt                                                                                                          >  2:772727/1‑149 (MQ=255)                                    aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGAc                                                                                                       <  2:70294/151‑1 (MQ=255)                                      tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAATATA‑TT‑TCTTG‑TCTT‑TACAAAAAGaaaaaattaca                                                                                                      >  2:812229/1‑143 (MQ=255)                                       gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGa                                                                                                    <  1:925453/151‑1 (MQ=255)                                       gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAA‑TT‑GCTTG‑ACTT‑TACAAAAAAATaa                                                                                                             >  1:330924/1‑142 (MQ=255)                                         gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGaaa                                                                                                  <  1:611694/151‑1 (MQ=255)                                         gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGaaa                                                                                                  <  2:837071/151‑1 (MQ=255)                                         gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGaa                                                                                                   <  1:281056/150‑1 (MQ=255)                                            aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGa                                                                                               <  2:677319/151‑1 (MQ=255)                                              tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGATAAAGA‑TT‑GTTTT‑GGTT‑AAAAAAAAGATAAATTTGCTGAAAAg                                                                                                >  2:106947/1‑148 (MQ=255)                                               aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATaa                                                                                            <  1:811314/151‑1 (MQ=255)                                               aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATa                                                                                             <  2:778454/150‑1 (MQ=255)                                                gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATaaa                                                                                           <  1:207609/151‑1 (MQ=255)                                                            cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAA‑AT‑GCTTT‑CATT‑TAAAAAAAAAAAAATTTAAAGaaaaaaaaaaaaaaaa                                                                                     >  1:109712/1‑143 (MQ=255)                                                             tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAATAGA‑TT‑GCTTG‑TCTT‑TAAAAAAAAATAACATTACAGAAAAAATaaaaaaaaaacccac                                                                               >  1:18411/1‑142 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGA‑AT‑TCTTT‑TCTT‑ATAAAAAAAACAACATTATAGAAAAAAAAAAAAAaaagac                                                                                  >  1:352236/1‑144 (MQ=255)                                                               aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTCACAACGG‑TT‑GGATT‑ACTT‑TAAAAACAGAAAAAAATACAGGAACAAAAAAAAAAGAG‑ACAACggg                                                                           >  2:275160/1‑151 (MQ=255)                                                                 ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGgc                                                                          <  2:347417/150‑1 (MQ=255)                                                                 ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGA‑TT‑TCTTG‑ACTT‑TACAAAAAGATAAATTTACAGAAAAAATAAAAAAAAAA‑CCAGCggg                                                                           >  2:686152/1‑149 (MQ=255)                                                                  cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAA‑TT‑TCTTT‑ACTT‑TAAAAAAAAAAAAATTTACAGaaaaaaaaaaaaaaaaaacaacgg                                                                             >  1:553147/1‑137 (MQ=255)                                                                     aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAA‑TTATTTTA‑AATT‑AAAAAAAAAAAAAC‑TGACAGAAAAAAAAAAAAAAAAG‑ACAGCCGgcggcgg                                                                     >  1:580076/1‑147 (MQ=255)                                                                      aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAA‑TT‑GGGTG‑GCTT‑AACAAAAAAATAACTATACAGAAAAAAAAAAAAAAAAGACCAGCggcgccgg                                                                      >  2:385486/1‑144 (MQ=255)                                                                       tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGAGGAAATAAA‑TT‑GTTTT‑GCAT‑TACAAACAAACAAACTCACTGAAAAGATAAAAAAAAAA‑ACAGCCCGCGCTGGc                                                                    >  1:755446/1‑150 (MQ=255)                                                                           aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAA‑TT‑TTTTT‑ACTTAAAAAAAAAAAAAAATTGA‑AGAAAAAAAAAAAAAAAAA‑ACAACGGgcgggggggctt                                                                >  2:886358/1‑141 (MQ=255)                                                                                ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAA‑TT‑GCTTG‑GCTT‑AAAAAACAAATAAATTTACAGAAAAGAAAAAAAAAAAA‑CCAGCGGGGGCGGGCTCTTTaa                                                             >  2:754872/1‑148 (MQ=255)                                                                                  aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTT‑TTTTT‑ATTT‑AAAAAAAAAAAAAATT‑AAAGAAAAAAAAAAAAAAAAA‑ACAGCGGGGGCTGGCGCtt                                                                >  1:269994/1‑143 (MQ=255)                                                                                      gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                    <  1:677427/151‑1 (MQ=255)                                                                                              gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                            <  1:263455/151‑1 (MQ=255)                                                                                                  aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                        <  2:362981/151‑1 (MQ=255)                                                                                                  aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCGAAAAAA‑AT‑GCTTT‑TCTT‑AACAAAAAAAAAACTTTGCAGAAAAAAAAAAAAAAAAA‑GCAACGGGCGGCGGCTCCTGAAGAAGCCAAAATaaaaaaaggg                                        >  1:497124/1‑145 (MQ=255)                                                                                                   aaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAA‑TT‑TTTTGTAATA‑AAAAAAAAAAAAATTTGA‑AGAAAAGAAAAAAAAAAAG‑ACAGAGGGGGCTGGGTCTTTAAAAAGCGAAAATaaaaacaaggg                                       >  1:520781/1‑144 (MQ=255)                                                                                                    aTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAA‑TT‑GCTTT‑ACTT‑AACAAACAAAAAAATTTTCAGAAAAAAAAAAAAAAAAA‑CCAACGGGGGCTGGCTCCTGAAGAAGCCAAAATaaaaagtgg                                         >  2:549650/1‑143 (MQ=255)                                                                                                      gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                    <  2:269574/151‑1 (MQ=255)                                                                                                       ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAA‑TT‑GCTTT‑ACTT‑TACAAAAAAAAAACATTACAAAAAAAAAAAAAAAAAAA‑ACAACGGGCGCGGGCCCTTTAAGAAGCTAAATTaaaaaaaggggtta                                    >  1:778512/1‑140 (MQ=255)                                                                                                            aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                              <  1:772706/151‑1 (MQ=255)                                                                                                            aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAA‑TT‑TGTTT‑ACTT‑TACAAAAAGATAAATTTTCAGAAAAAAAAAAAAAAAAA‑ACAGCGGGGGCTGGGGCTTTAAGACGCTAAAATAAAAACTGGCGGTACGCgg                               >  1:491415/1‑150 (MQ=255)                                                                                                                   atTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAA‑TT‑TTTTG‑TATT‑TAAAAACAAAAAACTTTTAAAAAAAAATAAAAAAAAAA‑ACAGCGGGGGGCGGCCCTTTGAGAAGCCAAAATAAAAAGAGGCGGTGCGCGGAGGACGTg                       >  2:780259/1‑151 (MQ=255)                                                                                                                      ccGTAACAGGATGATCGTTATGGCGGGGGTGCGAAAAAA‑AT‑TTTTG‑GATT‑AACCAAAAGATATTTTTAGAGAAAAAATAAAAAAAAAG‑ACCGCGGGTGCTGGCTCTTTGCGACGCCAAATTAGAAACCGGCGTTTcgcgcgggacggggg                     >  2:287740/1‑139 (MQ=255)                                                                                                                        gTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGA‑AT‑GTTTG‑GCTT‑TAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAG‑ACAACGGGGGCCGGGCCTTGAAGAAGCTAAAATAGAAAATGGCGGTACGCGGAGTACGTggg                     >  2:403806/1‑148 (MQ=255)                                                                                                                                   atCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc       <  1:171049/151‑1 (MQ=255)                                                                                                                                         aTGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgta  <  1:76773/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGA‑TT‑GCTTG‑ACTT‑AACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAG‑CCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTA  >  NC_000913/987349‑987634

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