Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,950,699 A→G 3.0% intergenic (+42/‑118) RCAP_rcc02761 → / → RCAP_rcc02762 glutathione S‑transferase/4Fe‑4S ferredoxin, iron‑sulfur cluster binding protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,950,6990A→G 3.0% 10.6 ‑1.3 ‑1.0275intergenic (+42/‑118)RCAP_rcc02761/RCAP_rcc02762glutathione S‑transferase/4Fe‑4S ferredoxin, iron‑sulfur cluster binding protein
Reads supporting (aligned to +/- strand):  new base (6/2):  ref base (135/127):  total (142/133)
Fisher's exact test strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.18e-02

Only 100 of 291 total aligned reads displayed.

GGCAGATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTGAGGG  >  NC_014034/2950600‑2950799
                                                                                                    |                                                                                                        
ggCAGATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAACCGGGGGCTCTGCCCCCGGG.CCCCCGGGCTATTTGA.g                                                                                                        >  2:3814165/1‑101 (MQ=255)
 gCAGATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.gg                                                                                                       >  1:3115768/1‑101 (MQ=255)
   aGATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCa                                                                                                     <  1:1198992/101‑1 (MQ=255)
    gATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCaa                                                                                                    >  2:1430043/1‑101 (MQ=255)
    gATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCaa                                                                                                    <  2:3970311/101‑1 (MQ=255)
     aTCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.g                                                                                                  <  1:1114345/101‑1 (MQ=255)
      tCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTCA.GGCAA.Ga                                                                                                 >  1:2088247/1‑101 (MQ=255)
       cAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GAt                                                                                                >  1:3017633/1‑101 (MQ=255)
       cAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GAt                                                                                                <  1:3195013/101‑1 (MQ=255)
       cAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GAt                                                                                                <  2:642826/101‑1 (MQ=255)
        aGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATg                                                                                               <  1:1090752/101‑1 (MQ=255)
         gcggCCGGGGGTGCCGCCCTCGCCGCAAGATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGa                                                                                              <  1:2254506/98‑1 (MQ=255)
         gCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.Gagat                                                                                              >  2:3560720/1‑98 (MQ=255)
         gCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGa                                                                                              <  1:4580602/101‑1 (MQ=255)
           tGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGaaa                                                                                            >  1:909092/1‑101 (MQ=255)
             cccGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAaca                                                                                          <  2:335895/101‑1 (MQ=255)
              ccGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAacac                                                                                         >  1:814130/1‑101 (MQ=255)
               cGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACg                                                                                        <  1:4394794/101‑1 (MQ=255)
                 gggTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGcc                                                                                      <  2:3740638/101‑1 (MQ=255)
                  ggTGCCGCCCACGCCGCAATATGCGGACCTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGccc                                                                                     <  1:1798686/101‑1 (MQ=255)
                   gTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCCCTGGCCCCGGA.CCCCCGGGGTATTTTA.GGGGA.GAGGAAAAACGcccc                                                                                    >  1:3674670/1‑101 (MQ=255)
                    ggccgccCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCt                                                                                   <  2:2417159/100‑1 (MQ=255)
                       cgccCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTc                                                                                <  1:50533/101‑1 (MQ=255)
                         cccTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCAt                                                                              >  2:3189711/1‑101 (MQ=255)
                          ccTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGGGATTTTAGGGCAAGGAGAAAACAC.CCCCTTTCac                                                                              >  1:4032994/1‑100 (MQ=25)
                           cTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCt                                                                            <  2:2783058/101‑1 (MQ=255)
                            tCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCtt                                                                           <  2:3115568/101‑1 (MQ=255)
                             cgccgcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTg                                                                          <  1:3454/101‑1 (MQ=255)
                               ccgcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGcc                                                                        >  2:4514049/1‑101 (MQ=255)
                               ccgcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGcc                                                                        >  1:160444/1‑101 (MQ=255)
                                 gcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGcccc                                                                      <  2:5805/101‑1 (MQ=255)
                                 gcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGcccc                                                                      >  1:637741/1‑101 (MQ=255)
                                 gcAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGcccc                                                                      >  2:1624024/1‑101 (MQ=255)
                                  cAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCa                                                                     >  1:3164869/1‑101 (MQ=255)
                                  cAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCa                                                                     <  1:2777654/101‑1 (MQ=255)
                                    atatGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCaaa                                                                   >  2:2348276/1‑101 (MQ=255)
                                      atGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAata                                                                 >  2:3288385/1‑101 (MQ=255)
                                          ggAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATccc                                                             <  1:1017163/101‑1 (MQ=255)
                                           gAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.g                                                           >  1:3310227/1‑101 (MQ=255)
                                            aaCTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.gg                                                          >  1:2038822/1‑101 (MQ=255)
                                             aCTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGTTATTTGA.GGCAA.GATGAAACACGCCCATTTCATCTTGCCCCAAATATCCC.ggg                                                         >  2:3048645/1‑101 (MQ=255)
                                                tGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGt                                                      >  1:2747834/1‑101 (MQ=255)
                                                 gATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGACCCCCCGGGGCTTTTTG.GGCAA.AAGGAAAAACCCCCCCTCCTCCTTCCCCCCAAAACACC.gggggg                                                      >  1:1619654/1‑100 (MQ=9)
                                                  aTTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGa                                                    <  2:4635315/101‑1 (MQ=255)
                                                  aTTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGa                                                    <  2:2466163/101‑1 (MQ=255)
                                                  aTTTCTGAAGCAAGCCGGCGGCTCTGCCCCCGGG.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGa                                                    <  2:4747040/101‑1 (MQ=255)
                                                    ttCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAAt                                                  >  1:1913584/1‑101 (MQ=255)
                                                    ttCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCA.GGGGGCGAAt                                                  >  2:1157164/1‑101 (MQ=255)
                                                      cTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAAttt                                                <  1:3729744/101‑1 (MQ=255)
                                                      cTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAAttt                                                >  2:1517489/1‑101 (MQ=255)
                                                           gcaagcCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGt                                           >  2:3011852/1‑101 (MQ=255)
                                                           gcaagcCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGt                                           >  1:4310411/1‑101 (MQ=255)
                                                           gcaagcCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGt                                           >  1:417555/1‑101 (MQ=255)
                                                             aagcCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCa                                         <  2:4310164/101‑1 (MQ=255)
                                                               gcCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAgg                                       <  2:3218277/101‑1 (MQ=255)
                                                                 cGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCa                                     <  2:1590792/101‑1 (MQ=255)
                                                                 cGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCa                                     <  2:2709753/101‑1 (MQ=255)
                                                                  gggggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCaa                                    >  1:4186376/1‑101 (MQ=255)
                                                                  gggggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCaa                                    >  2:4607829/1‑101 (MQ=255)
                                                                    gggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAag                                  <  1:2642116/101‑1 (MQ=255)
                                                                    gggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAag                                  >  2:2896414/1‑101 (MQ=255)
                                                                     ggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAaga                                 <  2:2861373/101‑1 (MQ=255)
                                                                     ggCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAaga                                 <  2:1619554/101‑1 (MQ=255)
                                                                         ctGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAgggg                             >  2:1767295/1‑101 (MQ=255)
                                                                          tGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAgggg.g                           >  2:1804674/1‑101 (MQ=255)
                                                                          tGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAgggg.g                           >  1:1681552/1‑101 (MQ=255)
                                                                          cttcgcctca.CCCCCGGGATATTTGG.GGCAA.GAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCG.AAGGCCAAGAgggg.gg                          >  2:2067505/10‑100 (MQ=25)
                                                                           gCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAAAgggg.gg                          >  1:155802/1‑100 (MQ=255)
                                                                           gCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGGGAATTTGCCGTCAGGCAAAGAgggg.gg                          >  2:4065671/1‑100 (MQ=255)
                                                                            tcgcctca.CCCCCGGGATATTTGG.GGCAA.GAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCG.AAGGCCAAGAgggg.ggca                        <  2:3665917/94‑4 (MQ=25)
                                                                            cccccGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCa                         <  1:3803636/101‑1 (MQ=255)
                                                                             ccccGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAg                        <  1:2420912/101‑1 (MQ=255)
                                                                                cGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTc                     >  1:4526217/1‑101 (MQ=255)
                                                                                 ggA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTcc                    <  2:1344659/101‑1 (MQ=255)
                                                                                 ggA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTcc                    >  1:2363703/1‑101 (MQ=255)
                                                                                  gA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTccc                   >  1:3589380/1‑101 (MQ=255)
                                                                                     tccccGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTccccc                 <  2:2491365/100‑1 (MQ=255)
                                                                                      ccccGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCt                <  1:945291/101‑1 (MQ=255)
                                                                                      ccccGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCt                <  1:2959554/101‑1 (MQ=255)
                                                                                       cccGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTg               <  1:3395368/101‑1 (MQ=255)
                                                                                        ccGGGATATTTGG.GGCAA.GAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCG.AAGGCCAAGCGGGGCGCAGCGCCCCCttt              >  2:754805/1‑99 (MQ=18)
                                                                                         cGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGcc             <  2:3310025/101‑1 (MQ=255)
                                                                                          gggATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGccc            <  2:125936/101‑1 (MQ=255)
                                                                                          gggATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGccc            >  2:2756046/1‑101 (MQ=255)
                                                                                           ggATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGcccc           >  1:2890378/1‑101 (MQ=255)
                                                                                           ggATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGcccc           <  2:4579528/101‑1 (MQ=255)
                                                                                            gATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCg          <  1:1430137/101‑1 (MQ=255)
                                                                                             atatTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCgg         >  2:2635041/1‑101 (MQ=255)
                                                                                             atatTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCgg         >  2:4378238/1‑101 (MQ=255)
                                                                                              tatTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGc        <  2:3417831/101‑1 (MQ=255)
                                                                                              tatTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGc        <  2:2771842/101‑1 (MQ=255)
                                                                                              tatTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGc        >  2:4466670/1‑101 (MQ=255)
                                                                                               atTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCt       >  1:56973/1‑101 (MQ=255)
                                                                                                tttGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTg      <  2:4186130/101‑1 (MQ=255)
                                                                                                tttGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCAC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTg      >  2:4103708/1‑101 (MQ=255)
                                                                                                 ttGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTga     >  1:4616179/1‑101 (MQ=255)
                                                                                                 ttGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTga     <  2:2188631/101‑1 (MQ=255)
                                                                                                   gA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTgagg   >  2:928183/1‑101 (MQ=255)
                                                                                                   gA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTgagg   >  2:1095417/1‑101 (MQ=255)
                                                                                                   gA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTgagg   <  2:2900780/101‑1 (MQ=255)
                                                                                                    |                                                                                                        
GGCAGATCAGCTGCCCGGGGTGCCGCCCTCGCCGCAATATGCGGAACTTGATTTCTGAAGCAAGCCGGGGGCTCTGCCCCCGGA.CCCCCGGGATATTTGA.GGCAA.GATGAAACACGCCCCTTTCATCTTGCCCCAAATATCCC.GGGGGTGAATTTGCCGTCAGGCAAAGAGGGG.GCAGCTCCCCCTGCCCCGGCTGAGGG  >  NC_014034/2950600‑2950799

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG