Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,622,021 +G 1.2% coding (104/717 nt) RCAP_rcc03420 → LrgB family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,622,0211.→G 1.2% 3.7 ‑3.3 ‑0.0323coding (104/717 nt)RCAP_rcc03420LrgB family protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (139/180):  total (143/180)
Fisher's exact test strand distribution p-value = 3.75e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

Only 100 of 393 total aligned reads displayed.

AATGACGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTATTTCGAG  >  NC_014034/3621921‑3622121
                                                                                                     |                                                                                                    
aatgaCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCt                                                                                                       <  2:6005306/101‑1 (MQ=255)
 atgaCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.g                                                                                                     <  1:4982395/101‑1 (MQ=255)
   gaCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGc                                                                                                   <  1:6227242/101‑1 (MQ=255)
    aCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.gggc                                                                                                  >  2:1205585/1‑99 (MQ=255)
    aCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.gggc                                                                                                  >  1:4710685/1‑99 (MQ=255)
    aCGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGcc                                                                                                  <  2:1590238/101‑1 (MQ=255)
     cgcgGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGGGGCGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGgcc                                                                                                 >  1:839223/1‑101 (MQ=255)
     cgcgGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGgcc                                                                                                 >  2:1136511/1‑101 (MQ=255)
       cggcgCTGGTCTGGCCGTATCTCGCCGCGCGGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCgg                                                                                               <  2:831800/101‑1 (MQ=255)
        ggcgCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCTGGGCCCgg                                                                                               >  1:5571576/1‑101 (MQ=255)
        ggcgCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGgc                                                                                              <  2:3665845/101‑1 (MQ=255)
           gcTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCa                                                                                           <  2:5355657/101‑1 (MQ=255)
             tGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCCGGGCAAt                                                                                         >  2:2019705/1‑101 (MQ=255)
             tGATCTGGTAGTATCCCGCATCGCAGCCGCTGTTGTGGCTGACGACGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAAt                                                                                         <  2:5629066/101‑1 (MQ=255)
              gATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCagac                                                                                        <  2:6730998/101‑4 (MQ=255)
                  tGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCAATCCCgc                                                                                    >  2:5298250/1‑101 (MQ=255)
                    gTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGgccg                                                                                  <  1:2323852/101‑1 (MQ=255)
                     tCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGCCAGTCTGTATCGGCT.GGGCCGGGGCCATCCCgccgc                                                                                 >  1:4400047/1‑101 (MQ=255)
                      ccTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGgccgcc                                                                                <  2:5446039/101‑1 (MQ=255)
                       cTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGGGCAATCCCgcccccc                                                                               >  1:4542711/1‑100 (MQ=255)
                        tATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCCGCGCCATCCCGGCGCCaa                                                                              >  2:3019324/1‑101 (MQ=255)
                         aTCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAAt                                                                             <  1:3499619/101‑1 (MQ=255)
                          tctcGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCTGGGCCCGGC.CAATCCGgcccgccatc                                                                            >  2:1498844/1‑94 (MQ=255)
                          tctcGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCCATCCGggcgcccaac                                                                            >  2:3742655/1‑97 (MQ=255)
                          tctcGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATc                                                                            <  2:5571286/101‑1 (MQ=255)
                            tcgccacgcAGCCGCGGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCCGGGCGATCCCGCCGCCCATCCg                                                                          >  1:6871165/1‑101 (MQ=255)
                             cgcctcgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCgg                                                                         <  2:6660312/101‑1 (MQ=255)
                               ccacgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTg                                                                       <  1:1396530/98‑1 (MQ=255)
                                  cgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCGAACCGGGCCCCCACCCggcggag                                                                    >  2:6601814/1‑96 (MQ=255)
                                  cgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATg                                                                    <  1:2253536/101‑1 (MQ=255)
                                  cgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATg                                                                    <  2:2354769/101‑1 (MQ=255)
                                  cgcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATg                                                                    <  2:3810523/101‑1 (MQ=255)
                                   gcAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCAATCCGGCCGCCCATCCGGGGATGc                                                                   >  1:1048960/1‑101 (MQ=255)
                                    cAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCt                                                                  <  1:4140003/101‑1 (MQ=255)
                                      gCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCAGCGCCAACCGGGCGCCCAGCCCGTGTTGCCgg                                                                >  2:1903446/1‑101 (MQ=255)
                                       ccGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGGGCAATCCCGCCGCCAATCCCGGGATGCTggg                                                               >  2:7058604/1‑100 (MQ=255)
                                           tgttgtGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCCATCCGGGCGCCCATCCGGTGATGCTGGCgggg                                                           >  2:2555532/1‑99 (MQ=255)
                                            gttgtGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCCGCGCAATCCCGCCGCCCAACCGGGGGTGCTGGCggggg                                                          >  1:5506443/1‑98 (MQ=255)
                                             ttgtGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGGGATGCTGGGgggggt                                                         >  1:4249993/1‑97 (MQ=255)
                                              tgtgGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCCATACGGCCGCCAATCCGGTGATGCTGGGggggggt                                                        >  2:5251072/1‑96 (MQ=255)
                                              tgtgGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCCATCCGGGGATGCTGGCggggggt                                                        >  1:5277906/1‑96 (MQ=255)
                                              tgtgGCCGACGGGGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGAtt                                                        <  1:952590/101‑1 (MQ=255)
                                              tggggCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGAtt                                                        <  2:4975270/98‑1 (MQ=255)
                                                tgGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGGCCGGCGCAATCCGGCCGCCAATCCGGGGATGCTGGCGGTGAtttt                                                      >  2:6899172/1‑101 (MQ=255)
                                                 ggCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGAttttt                                                     <  2:2873364/101‑1 (MQ=255)
                                                 ggCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGAttttt                                                     <  2:5144415/101‑1 (MQ=255)
                                                 ggCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGAttttt                                                     <  2:4769485/101‑1 (MQ=255)
                                                  gCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTc                                                    <  2:2669398/101‑1 (MQ=255)
                                                  gCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTc                                                    <  2:6493705/101‑1 (MQ=255)
                                                   cTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGAGTTTTCt                                                   >  2:1770921/1‑101 (MQ=255)
                                                    tGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTg                                                  <  1:6043126/101‑1 (MQ=255)
                                                      aCGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGc                                                <  1:1625836/101‑1 (MQ=255)
                                                          cGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGtg                                            <  2:6249056/101‑1 (MQ=255)
                                                            aCGCTGTTTGCCTATTGGGCGGGGGACATTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGtgct                                          <  1:2944897/101‑1 (MQ=255)
                                                             cccttttggCCTAGGGGGCGGGGGACAGTCTGTATCGGCT.GCCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGtgctg                                         <  2:2299549/93‑1 (MQ=255)
                                                             cGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGtgctg                                         <  2:4881662/101‑1 (MQ=255)
                                                             cGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGtgctg                                         <  1:2531933/101‑1 (MQ=255)
                                                                tGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCtt                                      <  2:768208/101‑1 (MQ=255)
                                                                tGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGCGCTGCtt                                      <  1:394806/101‑1 (MQ=255)
                                                                 gTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTg                                     <  1:20855/101‑1 (MQ=255)
                                                                 gTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTg                                     >  2:2158604/1‑101 (MQ=255)
                                                                  tttGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAAACCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCtttg                                    >  2:2107718/1‑99 (MQ=255)
                                                                   ttGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGc                                   >  1:7071128/1‑101 (MQ=255)
                                                                   ttGCCTATGGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGc                                   <  2:1986493/101‑1 (MQ=255)
                                                                    tGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCTGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCa                                  <  1:6529028/101‑1 (MQ=255)
                                                                    tGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGCa                                  >  2:1670667/1‑101 (MQ=255)
                                                                      ccTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATc                                <  1:64688/101‑1 (MQ=255)
                                                                       cTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGCATCa                               >  1:5930210/1‑101 (MQ=255)
                                                                       cTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGCATCa                               >  1:6628481/1‑101 (MQ=255)
                                                                         agtGGGCGGGGGACAGTCTGTATCGGCT.GCCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCAcc                             <  1:3815681/99‑1 (MQ=255)
                                                                          ttGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCg                            <  2:5506171/101‑1 (MQ=255)
                                                                            gggCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGc                          <  1:5839586/101‑1 (MQ=255)
                                                                            gggCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGc                          <  1:6687368/101‑1 (MQ=255)
                                                                            gggCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGGGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGc                          >  1:3275432/1‑101 (MQ=255)
                                                                             ggCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGAGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCa                         >  1:3376288/1‑101 (MQ=255)
                                                                              gCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCAc                        <  1:5203211/101‑1 (MQ=255)
                                                                              gCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCAc                        <  2:3993486/101‑1 (MQ=255)
                                                                              gCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCAc                        <  1:1996111/101‑1 (MQ=255)
                                                                               cGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACg                       <  1:1221727/101‑1 (MQ=255)
                                                                                gggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGGCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGGTTTGCATCCCAGGCACGc                      >  2:1360701/1‑101 (MQ=255)
                                                                                gggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTTTGGCCGTGCCGCTTTGCATCACCGGCCCGc                      >  1:5698972/1‑101 (MQ=255)
                                                                                gggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTTGCCGTGCTGCTTGGCATCACCGGCAAGc                      >  2:3357545/1‑101 (MQ=255)
                                                                                 ggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTTGCATCACCGGCACGcc                     >  2:3131818/1‑101 (MQ=255)
                                                                                 ggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGcc                     <  1:6135431/101‑1 (MQ=255)
                                                                                 ggggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGcc                     >  2:3887551/1‑101 (MQ=255)
                                                                                  gggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGccc                    <  1:3681382/101‑1 (MQ=255)
                                                                                  gggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCCGCCGCCAATCCGGGGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCATGccc                    >  2:4528102/1‑101 (MQ=255)
                                                                                   ggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCt                   <  1:3074804/101‑1 (MQ=255)
                                                                                   ggACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCt                   <  1:2793961/101‑1 (MQ=255)
                                                                                    gACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTa                  <  2:5859522/101‑1 (MQ=255)
                                                                                       aGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCg               <  1:1709297/101‑1 (MQ=255)
                                                                                         tCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCCCCCggg             >  1:1027167/1‑101 (MQ=255)
                                                                                            gTATCGGCT.GGCCCGGCGCAATCCGGCCGCCCATCCGGTGATGCTGGCGGGGATTTTTCTGGCCGTGCTGCTTGGGCTCACCGGCACGCCCCACCGCGAgg          >  2:4237457/1‑100 (MQ=255)
                                                                                             tATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTa         >  1:4622586/1‑101 (MQ=255)
                                                                                              aTCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCCACCcgggggat        >  1:6313155/1‑96 (MQ=255)
                                                                                              aTCGGCG.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGCTGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTAt        <  1:5096009/101‑1 (MQ=255)
                                                                                               tCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTAtt       <  1:6979668/101‑1 (MQ=255)
                                                                                                 ggCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTATTTc     <  2:4634917/101‑1 (MQ=255)
                                                                                                  gCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTATTTCg    >  1:6775504/1‑101 (MQ=255)
                                                                                                    t.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTATTTCGAg  >  2:4260151/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
AATGACGCGGCGCTGATCTGGTCCTATCTCGCCTCGCAGCCGCTGTTGTGGCTGACGGCGACGCTGTTTGCCTATTGGGCGGGGGACAGTCTGTATCGGCT.GGCCCGGCGCAATCCGGCCGCCAATCCGGTGATGCTGGCGGTGATTTTTCTGGCCGTGCTGCTTGGCATCACCGGCACGCCCTACCGGGAGTATTTCGAG  >  NC_014034/3621921‑3622121

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 41 ≤  ATCG