New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 74479176 (0.170)73 (0.160) 5/198 39.0 0.40 intergenic (+176/+36) idi2/bchO isopentenyl‑diphosphate delta‑isomerase/magnesium‑chelatase BchO
?NC_014034 = 2330174 143 (0.310)intergenic (+149/+163) RCAP_rcc02161/mrcA hypothetical protein/penicillin‑binding protein 1A
Rejected: Position hash score below cutoff.

CAAGATGAAATGAATGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGATCAG  >  NC_014034/744692‑744889
                                                                                                   |                                                                                                  
cAAGATGAAATGAATGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCACAATACGCCGGGGGGTGCGGGGGGCGGAGccccc                                                                                                   <  1:8205333/101‑1 (MQ=255)
cAAGATGAAATGAATGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGccccc                                                                                                   <  1:5853030/101‑1 (MQ=255)
 aaGATGAAATGAATGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGcccccc                                                                                                  <  2:7396048/101‑1 (MQ=255)
      gaaatgaaTGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCggg                                                                                             >  2:2085745/1‑101 (MQ=255)
          tgaatgaCTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGc                                                                                         >  2:6380914/1‑101 (MQ=255)
            aatgaCTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGccc                                                                                       >  2:4460499/1‑100 (MQ=255)
               gaCTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCGCgg                                                                                    >  2:7407418/1‑101 (MQ=255)
                   ttCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGGCt                                                                                >  2:1514966/1‑101 (MQ=255)
                    tCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCCGGcccg                                                                               >  1:6773094/1‑99 (MQ=255)
                     gaccttGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGc                                                                              >  2:4667884/6‑101 (MQ=255)
                       gTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGccc                                                                            >  1:8112208/1‑101 (MQ=255)
                        tCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCCAGGGCCTGCCCg                                                                           >  2:5341697/1‑101 (MQ=255)
                         cTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCCGGCCCGcccca                                                                          >  2:3060257/1‑99 (MQ=255)
                             ggcggTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCCAGGCg                                                                      >  2:5128411/1‑101 (MQ=255)
                                                     gCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCCGCCCGAGGCGAAATCTCCCCCCCCGAGCCACTcc                                              >  2:1168703/1‑101 (MQ=18)
                                                         ttCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGCCCGGGGCGAAAACTCACCCCTCGAGCCAGGCCAGaa                                          >  2:2111519/1‑101 (MQ=21)
                                                           ccTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTTACCCCTCGAGCCGGCCCAggaca                                        >  2:8928160/1‑97 (MQ=18)
                                                           ccTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCCGCCCGAGGCGAAATCTCACCCCCCGAGCCAGTCCAGAatc                                        >  1:8984600/1‑101 (MQ=255)
                                                           ccTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGGCCCCGGCGGAATCCCACCCCTCGAGCCAGTCCAGAatc                                        >  1:2020886/1‑101 (MQ=21)
                                                           ccTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCCCCGGCCCGCCCGGGGCCAAATCCCACCCCCCGCGCCCGCCCCGAatc                                        >  1:361955/1‑101 (MQ=14)
                                                              tGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCCCCCCCCGCGCCCGTCCAGAatcacc                                     >  2:7402720/1‑99 (MQ=21)
                                                                 cAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCCGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAAACATCCCt                                  >  1:5836933/1‑101 (MQ=255)
                                                                   aaaTACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCCCGAGCCAGTCCAGAATCATCCCTGa                                >  2:3958028/1‑101 (MQ=255)
                                                                    aaTACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGTGGGGCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGaa                               >  1:8253354/1‑101 (MQ=255)
                                                                     aTACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCCGGCCCGCCCGAGGCGAAATCTCACCCCTCGAGCCCGTCCCGGATCATCCCTGgaa                              >  2:5642983/1‑101 (MQ=25)
                                                                       aCGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGc                            >  1:3791112/1‑101 (MQ=255)
                                                                       aCGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCCGAAAGc                            >  2:5142798/1‑101 (MQ=255)
                                                                         gCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCCCGGGGCGGCCCCAGGGGAAAACTCACCCCCCGAGCCAGTTCAGCATCAGCCCCGaaattcc                          >  2:7535877/1‑97 (MQ=17)
                                                                            ggggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCCGTCCAGAATCATCCCTGAAAGCCCATc                       >  2:8124001/1‑101 (MQ=255)
                                                                            ggggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCCAGCCAGTCCAGAATCATCCCTGAAAGCCCATc                       >  1:8565512/1‑101 (MQ=255)
                                                                             gggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCCGCCCAGAATCATCCCTGAAAGCCCCTcc                      >  2:8146669/1‑101 (MQ=255)
                                                                             gggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCCAGCCAGTCCAGAATCATCCCTGAAAGccccccc                      >  2:8528997/1‑97 (MQ=255)
                                                                             gggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCCCGAGCCCGTCCAGAATCATCCCTGAAAGCCCATcc                      >  1:3322825/1‑101 (MQ=255)
                                                                             gggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCCCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATcc                      >  2:3492123/1‑101 (MQ=255)
                                                                              ggggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCcg                     >  1:2380372/1‑101 (MQ=255)
                                                                               gggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATcccc                    >  2:2576048/1‑99 (MQ=255)
                                                                               gggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCcgc                    >  1:8635070/1‑101 (MQ=255)
                                                                               gggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCcgc                    >  1:2476275/1‑101 (MQ=255)
                                                                               gggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCcgc                    >  2:5717759/1‑101 (MQ=255)
                                                                               gggTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCCGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCcgc                    >  2:2655454/1‑101 (MQ=255)
                                                                                 gTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCgg                  >  2:2238676/1‑101 (MQ=255)
                                                                                 gTGCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGCCCGAGGCGAAATCTCAGCCCTCCAGCCAGTCCAGAATCATCCCTGAAAGCCCACCCGCgg                  >  1:2538715/1‑101 (MQ=38)
                                                                                 gTGCGGGGGGCGGAGCCCCCCGCGGGCCGCCTCGGGCCTGCCCGAGGCGCAATCTCACCCCTCGCGCCAGCCCACGATGCTCCCTGAAAGCCCATCcccgc                  >  2:9173179/1‑100 (MQ=25)
                                                                                 gTGCAGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCgg                  >  2:6362770/1‑101 (MQ=255)
                                                                                  tGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGGCTGCCCGGGGCGAAATCTCACCCCTCGAGCCCGTCCAGAATCATCCCTGAAAGGCCATCCGCGGc                 >  2:1046761/1‑101 (MQ=37)
                                                                                  tGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCCGAAAGCCCATCCGCGGc                 >  1:2276904/1‑101 (MQ=255)
                                                                                   gCGGGGGGCGGAGCCCCCCGCGGGTCGGCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCa                >  2:8137368/1‑101 (MQ=255)
                                                                                   gCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCa                >  1:9123440/1‑101 (MQ=255)
                                                                                   gCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCCGCa                >  2:2022012/1‑101 (MQ=255)
                                                                                   gCGGGGGGCGGAGCCCCCCGCGGGGCGCCTCGGGCCTGCCCGCGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCa                >  1:754691/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGGCTGCCCGAGGCGAAATCTCACCCCTCGCGCGGGTCCAGAACCATCCCTGAAAGCCCATCCGGGGCAAc              >  2:5564567/1‑101 (MQ=25)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  1:2123448/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  2:5147446/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  1:1275112/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  1:7116802/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  1:7795537/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAc              >  2:1092770/1‑101 (MQ=255)
                                                                                     ggggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAAACCCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCCAc              >  1:7676787/1‑101 (MQ=255)
                                                                                      gggggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGACCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAAcc             >  1:8253842/1‑101 (MQ=255)
                                                                                         ggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCg          >  1:4215601/1‑101 (MQ=255)
                                                                                         ggCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCg          >  2:8311122/1‑101 (MQ=255)
                                                                                          gCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGt         >  2:2292803/1‑101 (MQ=255)
                                                                                          gCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGt         >  1:241286/1‑101 (MQ=255)
                                                                                           cGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTg        >  2:851923/1‑101 (MQ=255)
                                                                                           cGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTg        >  1:8434188/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGGGAAATCTCACCCCTCGAGCCAGACCAGAATCATCCCTGAAAGCCCATCCGCGGCaaaccctggg       >  1:2211246/1‑93 (MQ=38)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCCGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCCCGTgg       >  1:2289792/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  1:1116284/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  1:2427221/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  2:6882453/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  1:287619/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  2:6353634/1‑101 (MQ=255)
                                                                                            ggAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTgg       >  1:3299662/1‑101 (MQ=255)
                                                                                             gAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGa      >  1:3461889/1‑101 (MQ=255)
                                                                                             gAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGa      >  1:1781942/1‑101 (MQ=255)
                                                                                             gAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGa      >  2:2367529/1‑101 (MQ=255)
                                                                                                ccccccGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGATCa   >  2:7897841/1‑101 (MQ=255)
                                                                                                 cccccGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGATCAg  >  1:5076585/1‑101 (MQ=255)
                                                                                                   |                                                                                                  
CAAGATGAAATGAATGACTTTCCGTCTGCGGCGGTACAGCCTGTCCCACCCCTGCTTTTCCTTGGCAAAATACGCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCGGGTCGCCTCGGGCCTGCCCGAGGCGAAATCTCACCCCTCGAGCCAGTCCAGAATCATCCCTGAAAGCCCATCCGCGGCAACCTCGTGGATCAG  >  NC_014034/744692‑744889

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 28 ≤  ATCG  < 38 ≤  ATCG