New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	= 1253890	93 (0.220)	72 (0.180)	4/178	38.2	0.29	intergenic (+48/+43)	RCAP_rcc01181/hisH	hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
	?	NC_014034	2590190 =	257 (0.560)					intergenic (‑89/‑167)	glnE/RCAP_rcc02415	glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

CGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG...................................................................................................  >  NC_014034/1253782‑1253890 ....................................................................................................GGGGGTGCGGGGGGCAGCGCCCCCGCCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTT  >  NC_014034/2590190‑2590288                                                                                                                                                                                                                           GGTGCTGGGCAAGTGCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGG                                                                                                    <  2:3666647/101‑1 (MQ=0)          GGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGG                                                                                                    <  2:6300064/101‑1 (MQ=1)          GGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGG                                                                                                    <  1:6957159/101‑1 (MQ=1)          GGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGG                                                                                                    <  1:34238/101‑1 (MQ=1)          GGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGATCTCATCTTGCCCTAAATACCTCGGGGGTGCGG                                                                                                    <  1:7193819/101‑1 (MQ=1)           GTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCCAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGG                                                                                                   <  1:8287905/101‑1 (MQ=1)            TGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGG                                                                                                  >  1:2447492/1‑101 (MQ=1)            TGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGG                                                                                                  >  1:3188291/1‑101 (MQ=1)            TGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGG                                                                                                  <  2:8711436/101‑1 (MQ=1)            TGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGG                                                                                                  <  2:6057410/101‑1 (MQ=1)             GCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGG                                                                                                 >  2:3409708/1‑101 (MQ=1)             GCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGG                                                                                                 <  2:4055508/101‑1 (MQ=1)              CTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGc                                                                                                <  2:7796760/101‑2 (MQ=1)              CTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGc                                                                                                <  1:4417576/101‑2 (MQ=1)               TGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcg                                                                                               >  1:4186448/1‑99 (MQ=1)               TGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcg                                                                                               >  2:4733060/1‑99 (MQ=1)                GGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgg                                                                                              <  1:6735646/101‑4 (MQ=1)                GGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgg                                                                                              >  2:5855873/1‑98 (MQ=1)                GGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgg                                                                                              >  2:6707128/1‑98 (MQ=1)                GGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgg                                                                                              >  1:1781962/1‑98 (MQ=1)                 gaaAAGCGCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgga                                                                                             >  2:5476609/4‑97 (MQ=0)                 GGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgga                                                                                             <  2:6302869/101‑5 (MQ=1)                 GGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgga                                                                                             >  1:149369/1‑97 (MQ=1)                 GGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcgga                                                                                             <  2:8793931/101‑5 (MQ=1)                 GGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGGGCGGGGGcgga                                                                                             >  1:7338919/1‑97 (MQ=1)                  GCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggag                                                                                            >  2:3033536/1‑96 (MQ=1)                  GCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggag                                                                                            <  2:6864967/101‑6 (MQ=1)                  GCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggag                                                                                            >  2:1055552/1‑96 (MQ=1)                   CAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagc                                                                                           <  2:1441839/101‑7 (MQ=1)                   CAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagc                                                                                           >  1:3855006/1‑95 (MQ=1)                   CAAGCTCGGCTTCCAGACCGGCGGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagc                                                                                           <  1:875040/101‑7 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          >  2:5947550/1‑94 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          >  2:288260/1‑94 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          <  2:3754412/101‑8 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          <  2:4968898/101‑8 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          <  1:7826247/101‑8 (MQ=1)                    AAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcc                                                                                          >  2:6728619/1‑94 (MQ=1)                     AGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccc                                                                                         <  1:2681763/101‑9 (MQ=1)                     AGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccc                                                                                         <  1:883868/101‑9 (MQ=1)                     AGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccc                                                                                         <  1:3920448/101‑9 (MQ=1)                     AGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccc                                                                                         <  1:7497415/101‑9 (MQ=1)                     AGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccc                                                                                         <  1:5692483/101‑9 (MQ=1)                      GCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcccc                                                                                        <  1:8769849/101‑10 (MQ=1)                      GCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcccc                                                                                        <  1:4094477/101‑10 (MQ=1)                      GCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagcccc                                                                                        <  1:5937093/101‑10 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       <  1:2610486/101‑11 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       <  1:4450793/101‑11 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       >  1:4801322/1‑91 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       >  2:179280/1‑91 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       <  1:7795931/101‑11 (MQ=1)                       CTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGcggagccccc                                                                                       <  2:6681810/101‑11 (MQ=1)                         CGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGC                                                                                    >  1:4752015/1‑101 (MQ=11)                         CGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGC                                                                                    >  2:6725337/1‑101 (MQ=11)                         CGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGC                                                                                    >  1:7840961/1‑101 (MQ=11)                           GCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCG                                                                                  >  2:9002162/1‑101 (MQ=32)                            CTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGG                                                                                 >  2:3427964/1‑101 (MQ=32)                             TTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGG                                                                                >  2:1331803/1‑101 (MQ=32)                             TTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGG                                                                                >  1:3841833/1‑101 (MQ=32)                               CCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTC                                                                              >  2:6048674/1‑101 (MQ=32)                               CCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTC                                                                              >  1:7662563/1‑101 (MQ=32)                                CAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCG                                                                             >  2:5588614/1‑101 (MQ=32)                                 AGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGC                                                                            >  2:2360179/1‑101 (MQ=32)                                 AGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGC                                                                            >  2:8666979/1‑101 (MQ=32)                                 AGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGC                                                                            >  1:7042996/1‑101 (MQ=32)                                  GACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCC                                                                           >  2:3738698/1‑101 (MQ=32)                                  GACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCC                                                                           >  2:6902188/1‑101 (MQ=32)                                  GACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCC                                                                           >  1:6489713/1‑101 (MQ=32)                                    CCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTC                                                                         >  2:5956447/1‑101 (MQ=32)                                    CCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTC                                                                         >  2:6361437/1‑101 (MQ=32)                                      GGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGG                                                                       >  2:8409318/1‑101 (MQ=32)                                      GGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGG                                                                       >  1:5366094/1‑101 (MQ=32)                                      GGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGG                                                                       >  1:8605054/1‑101 (MQ=32)                                      GGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGG                                                                       >  2:1605107/1‑101 (MQ=32)                                       GCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGC                                                                      >  1:1132787/1‑101 (MQ=32)                                       GCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGC                                                                      >  2:5213820/1‑101 (MQ=32)                                             CCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGCAAGC                                                                >  1:1163531/1‑101 (MQ=17)                                                 CGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGCAAGCCGAG                                                            >  1:8353035/1‑101 (MQ=17)                                                    GGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGC.GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGCAAGCCGAGGCG                                                         >  1:5023332/1‑101 (MQ=17)                                                                                                                                                                                                                  CGCCGCGATGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG...................................................................................................  >  NC_014034/1253782‑1253890 ....................................................................................................GGGGGTGCGGGGGGCAGCGCCCCCGCCGGGTCGCCTCGGCGCAAGCCGAGGCGAAATCTCTCCCCGACCCCTTGCGGCGCAATGTAAAAATCTTTCACATCCCCTCTT  >  NC_014034/2590190‑2590288 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG