Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,330,175 +C 76.6% intergenic (+150/+162) RCAP_rcc02161 → / ← mrcA hypothetical protein/penicillin‑binding protein 1A

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,330,1757.→C76.6% 58.9 0.0 ‑0.177intergenic (+150/+162)RCAP_rcc02161/mrcAhypothetical protein/penicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  new base (15/44):  ref base (4/14):  total (19/58)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.13e-01

GGCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......CGG.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCG  >  NC_014034/2330076‑2330245
                                                                                                        |                                                                                    
ggCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGgc....gg                                                                                      <  2:4574039/101‑2 (MQ=255)
  cTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggg                                                                             <  2:2924482/101‑1 (MQ=255)
   ttGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Gggg                                                                            >  2:6030360/1‑100 (MQ=255)
    tGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggg                                                                           >  2:2859929/1‑99 (MQ=255)
    tGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggc                                                                           <  2:6132951/101‑3 (MQ=255)
     ggCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......GGG.......g.cc                                                                  >  1:3210075/1‑101 (MQ=255)
      gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctc                                                                         <  2:5630159/101‑5 (MQ=255)
      gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctc                                                                         <  1:2416082/101‑5 (MQ=255)
      gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctc                                                                         <  2:2916443/101‑5 (MQ=255)
       ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctcc                                                                        <  2:2396785/101‑6 (MQ=255)
       ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctcc                                                                        <  1:1471946/101‑6 (MQ=255)
       ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctcc                                                                        <  1:126968/101‑6 (MQ=255)
       ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctcc                                                                        <  1:6796250/101‑6 (MQ=255)
        cTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctccg                                                                       <  1:5671554/101‑7 (MQ=255)
           cGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctccgccc                                                                    <  1:3598627/101‑10 (MQ=38)
            gCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctccgcccc                                                                   <  1:1686314/101‑11 (MQ=37)
            gCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctccgcccc                                                                   <  2:1500400/101‑11 (MQ=37)
            gCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......Ggggctccgcccc                                                                   <  2:5872949/101‑11 (MQ=37)
                      gcAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......GGG.......G.CGCCCCccccccccccccc                                                 >  1:3757946/1‑93 (MQ=37)
                                   gagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGcc                                           <  1:5200623/101‑1 (MQ=18)
                                   gagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGcc                                           <  1:3486040/101‑1 (MQ=18)
                                   gagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGcc                                           <  2:6435403/101‑1 (MQ=18)
                                   gagCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGGATTTGCGcc                                           >  2:6868732/1‑101 (MQ=18)
                                    agCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCa                                          <  2:4797578/101‑1 (MQ=18)
                                      cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAg                                        <  1:2860038/101‑1 (MQ=18)
                                      cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAg                                        <  1:1095973/101‑1 (MQ=18)
                                      cgccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAg                                        <  1:6421514/101‑1 (MQ=18)
                                       gccgcGCGGATCGGCCGATGGCGGAGAGTAATTTCGCCGCGCAGGCGCGGCGACCCGGTC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGa                                       <  2:3397096/101‑1 (MQ=9)
                                       gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGa                                       <  2:1123883/101‑1 (MQ=18)
                                       gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGa                                       <  2:5152467/101‑1 (MQ=18)
                                       gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGa                                       <  1:6755410/101‑1 (MQ=18)
                                       gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTGGCGCCAAGa                                       >  1:3081317/1‑101 (MQ=18)
                                       cccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGa                                       <  2:6776246/100‑1 (MQ=18)
                                        ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaa                                      <  2:6570006/101‑1 (MQ=18)
                                        ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaa                                      >  2:6725769/1‑101 (MQ=18)
                                         tctaCGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                     <  2:2914220/98‑1 (MQ=17)
                                         cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                     <  1:4275050/101‑1 (MQ=18)
                                         cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                     >  1:403950/1‑101 (MQ=18)
                                         cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                     <  1:2943822/101‑1 (MQ=18)
                                         cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                     <  1:2856323/101‑1 (MQ=18)
                                         cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGGATTTGCGCCAAGaaa                                     >  2:3257390/1‑101 (MQ=18)
                                          gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaaa                                    <  2:4957290/101‑1 (MQ=18)
                                          gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaaa                                    <  2:4150104/101‑1 (MQ=18)
                                          gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGaaaa                                    <  1:848701/101‑1 (MQ=18)
                                             cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                 <  2:1433623/101‑1 (MQ=18)
                                             cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGCATTTGCGCCAAGAAAAAGc                                 >  1:2685712/1‑101 (MQ=17)
                                              ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGcc                                <  2:4590684/101‑1 (MQ=18)
                                              ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGcc                                >  1:710528/1‑101 (MQ=17)
                                               gATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCg                               <  2:2910530/101‑1 (MQ=18)
                                                aTCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgg                              <  1:1312536/101‑1 (MQ=18)
                                                 tCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCggg                             <  1:6271252/101‑1 (MQ=18)
                                                    gCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcg                          <  1:1350339/101‑1 (MQ=17)
                                                     ccGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......GGGGCTCCCCC.CCCGCCCCCCCGGCGGATTTGCGCCAAGAAAAAACCCGGgcgc                         >  2:6912677/1‑101 (MQ=16)
                                                       gATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgcgc                       <  1:3217558/101‑1 (MQ=17)
                                                       gATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCCCACCCCCCGCGCATTTTCGCCCAGAAAAAACCGGGgcgcgc                       >  2:2011358/1‑101 (MQ=9)
                                                        aTGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCa                      >  1:3472288/1‑101 (MQ=17)
                                                         tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                     <  2:3130997/101‑1 (MQ=17)
                                                         tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCT.CCG.......CGCCCGCACCCCCGGCGTCTTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                     >  1:7114724/1‑101 (MQ=16)
                                                            cGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCGCCGC.......C.CCCCCCCCCCCGGCGGATTTGCGCCCAGAAAAAAACGGGGGGCGCATGCg                  >  1:5328718/1‑101 (MQ=9)
                                                            cGGAGAGGAATCTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCg                  <  1:4088172/101‑1 (MQ=17)
                                                             ggAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGCGGGGG.......CGGCGC....C.CCCGCACCCCCCGCGCATTTGCGGCAAGAAAAAGCCGGGGGGCCCATGCgg                 >  1:4588170/1‑101 (MQ=16)
                                                             ggAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCgg                 >  2:1453252/1‑101 (MQ=17)
                                                              gagaGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGGGCCGC.......C.CCCCCCCCCCCCCCGCATTTGCGGCCAGAAAAAACCGGGGCGCGCCTGCggg                >  1:1017519/1‑101 (MQ=9)
                                                                 agGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcc             <  1:4266993/101‑1 (MQ=17)
                                                                  ggAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......CGGGCTCCCCC.CCCGCACCCCCCCCGCACTTTCGCCCAGCAAAAGCCGGGGGGCGCCAGCGGGGccc            >  1:570590/1‑101 (MQ=9)
                                                                   gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGGGCCGC.......C.CCCCCACCCCCCGCGGATTTTCGCCCAGAAAAAAACCGGGCGCGCATGCggggggcc           >  2:2405363/1‑97 (MQ=9)
                                                                   gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc           <  1:1199662/101‑1 (MQ=17)
                                                                   gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc           <  1:2163672/101‑1 (MQ=17)
                                                                   gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc           <  2:4587931/101‑1 (MQ=17)
                                                                   gAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc           <  1:2354857/101‑1 (MQ=17)
                                                                    aaTTTCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCg          <  2:1212021/101‑1 (MQ=17)
                                                                      tttCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCACGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGa        >  1:488416/1‑101 (MQ=17)
                                                                       ttCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGAt       <  2:5138646/101‑1 (MQ=18)
                                                                        tCGCCGCGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATc      <  2:4758247/101‑1 (MQ=18)
                                                                         cgccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCa     <  2:4489916/101‑1 (MQ=18)
                                                                         cgccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCa     <  2:2300350/101‑1 (MQ=18)
                                                                           ccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc   >  1:2991479/1‑101 (MQ=18)
                                                                           ccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc   <  2:1017484/101‑1 (MQ=18)
                                                                           ccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc   <  2:4708265/101‑1 (MQ=18)
                                                                           ccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc   <  2:4634479/101‑1 (MQ=18)
                                                                           ccgcGCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc   <  1:584221/101‑1 (MQ=18)
                                                                            cgcgCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACcg  <  2:1320046/101‑1 (MQ=18)
                                                                            cgcgCAGGCGCGGCGACCCGGGC....GGGGGCTCCGC.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACcg  <  1:6910732/101‑1 (MQ=18)
                                                                                                        |                                                                                    
GGCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC....G.......CGG.......C.CCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCG  >  NC_014034/2330076‑2330245

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG