Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,947 +G 15.6% intergenic (+20/+23) RCAP_rcc01568 → / ← purF RES domain family protein/amidophosphoribosyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,699,9471.→G15.6% 37.6 ‑0.3 ‑0.177intergenic (+20/+23)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/7):  ref base (33/32):  total (38/39)
Fisher's exact test strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.53e-01

CATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGG  >  NC_014034/1699847‑1700040
                                                                                                     |                                                                                              
cATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGc                                                                                                 >  1:6300510/1‑101 (MQ=255)
cATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGc                                                                                                 >  1:4778464/1‑101 (MQ=255)
cATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGc                                                                                                 >  2:8658557/1‑101 (MQ=255)
cATCGGTTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGc                                                                                                 <  2:7246062/101‑1 (MQ=255)
 aTCTGGTGGTGTTTCGCTGGAGGCAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.t                                                                                               <  1:2985373/101‑1 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.t                                                                                               >  2:8623556/1‑101 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.t                                                                                               >  1:1653020/1‑101 (MQ=255)
  tCTGGTGCTTTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                              <  2:4768706/101‑1 (MQ=255)
  tCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                              <  1:3124630/101‑1 (MQ=255)
   cTGGTGCTGTTTCGCTGGAGCGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc                                                                                             <  1:7209621/101‑1 (MQ=255)
     ggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc.gc                                                                                          <  2:1747757/101‑1 (MQ=255)
      gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                         >  2:6605123/1‑101 (MQ=255)
      gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                         >  1:730466/1‑101 (MQ=255)
      gTGCTGTTGCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                         <  1:5076989/101‑1 (MQ=255)
       tGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                        >  2:6152652/1‑101 (MQ=255)
       gccggcttCGCGGGAGGGAGCCGCCGATCCCGGCGCAACTGCCCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                        <  1:1104981/95‑1 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGa                                                                                       <  2:3722665/101‑1 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGa                                                                                       >  2:8281716/1‑101 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGa                                                                                       >  2:4843015/1‑101 (MQ=255)
         gcgTTTTGCCGGAGTGGGCCGCTGATTTCGGCCCAACTGCCCGGCCTTCCCCCGCCTGCGGATTTTGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                      <  1:3220100/99‑1 (MQ=255)
         cTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                      <  1:5555904/101‑1 (MQ=255)
          tGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGCCCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                     <  1:1321894/101‑1 (MQ=255)
          tGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                     >  2:4133261/1‑101 (MQ=255)
          tGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                     >  2:2588736/1‑101 (MQ=255)
          tGTTTCGCTGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                     <  2:2279608/101‑1 (MQ=255)
          tGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGCTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                     <  1:6465970/101‑1 (MQ=255)
           ttttCGGTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    <  2:3217551/100‑1 (MQ=255)
           ggtgcGCTGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    <  1:1074735/97‑1 (MQ=255)
           gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    >  2:4083515/1‑101 (MQ=255)
           gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    <  1:4956025/101‑1 (MQ=255)
           gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    >  1:5856324/1‑101 (MQ=255)
           gTTTCGCGGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    <  1:6802460/101‑1 (MQ=255)
           gTTCCGCTGGGGGGGGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                    <  1:7348203/101‑1 (MQ=255)
            tttCGCTGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                   <  1:2487000/101‑1 (MQ=255)
            tttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                   >  2:5091208/1‑101 (MQ=255)
            tttCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                   <  2:125946/101‑1 (MQ=255)
            tttCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                   <  2:8630929/101‑1 (MQ=255)
             ttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGa                                                                                  >  2:2992761/1‑101 (MQ=255)
              tCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGac                                                                                 <  1:6306643/101‑2 (MQ=255)
              tCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGac                                                                                 <  1:5821519/101‑2 (MQ=255)
              tCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGac                                                                                 <  2:762969/101‑2 (MQ=255)
                gcgggagggggCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                               <  1:5575189/98‑4 (MQ=255)
                gcgggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                               <  2:3280004/98‑4 (MQ=255)
                gCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                               <  2:8675211/101‑4 (MQ=255)
                 cTGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccg                                                                              <  1:287534/101‑5 (MQ=255)
                 cTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccg                                                                              >  2:5548035/1‑97 (MQ=255)
                  tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                             >  1:1470284/1‑96 (MQ=255)
                  tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                             >  2:2709237/1‑96 (MQ=255)
                  gggggggggCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                             <  1:5049692/97‑6 (MQ=38)
                  gggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATAc                                                                             <  1:5603620/100‑1 (MQ=255)
                   ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                            >  2:5295153/1‑95 (MQ=255)
                   ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                            >  2:261125/1‑95 (MQ=255)
                   ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                            >  2:1049023/1‑95 (MQ=255)
                   ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                            >  2:2042141/1‑95 (MQ=255)
                   ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                            >  2:1638189/1‑95 (MQ=255)
                    gagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccc                                                                           >  1:4403182/1‑94 (MQ=255)
                     agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccg                                                                          >  1:394021/1‑93 (MQ=255)
                     agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccg                                                                          >  2:1557554/1‑93 (MQ=255)
                      gggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccggg                                                                         >  2:2573763/1‑92 (MQ=255)
                      gggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccgg                                                                         >  1:6369927/1‑92 (MQ=255)
                      gggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccgg                                                                         >  1:275126/1‑92 (MQ=255)
                       ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccgggg                                                                        >  2:3768735/1‑91 (MQ=255)
                       ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccgggg                                                                        >  2:2424353/1‑91 (MQ=255)
                                  tctcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCAt                                                             <  2:4909890/101‑1 (MQ=255)
                                                          gCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATc                                     <  2:8801331/101‑1 (MQ=255)
                                                                      aTTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCAc                         >  1:2720220/1‑101 (MQ=255)
                                                                          cGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGa                     >  1:2873326/1‑101 (MQ=255)
                                                                             ttGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATaa                  <  2:2720138/101‑1 (MQ=255)
                                                                             ttGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAACACCCCCTCAGCCCCTCACCGTGACGCCCTTGCCGTTCTGATCGGAAGGAGTCACCGGAAaa                  >  2:1462716/1‑101 (MQ=255)
                                                                                 ccccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg               >  2:7944172/11‑101 (MQ=25)
                                                                                  cccgccgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGc              >  2:4632147/10‑101 (MQ=25)
                                                                                   cggccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGcc             <  1:7944496/93‑1 (MQ=37)
                                                                                   ccgccgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGcc             >  2:918207/9‑101 (MQ=37)
                                                                                      ccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgg          >  1:6830421/6‑100 (MQ=37)
                                                                                      ccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCga          >  1:8491094/6‑101 (MQ=37)
                                                                                      ccgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCga          >  1:4736742/6‑101 (MQ=37)
                                                                                        gggtCGCCTCGGCGTCA.GCCGAGGCGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgaga        <  1:4606332/98‑1 (MQ=38)
                                                                                        gggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgaga        <  1:8040650/98‑1 (MQ=38)
                                                                                         gttCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa       <  2:919780/100‑1 (MQ=38)
                                                                                         ggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa       >  1:42341/3‑101 (MQ=38)
                                                                                          gtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa      <  2:1848004/100‑1 (MQ=38)
                                                                                           tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc     <  2:5895736/101‑1 (MQ=38)
                                                                                           tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc     <  2:5020368/101‑1 (MQ=38)
                                                                                            cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa    >  2:1707357/1‑101 (MQ=38)
                                                                                            cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa    >  2:3043356/1‑101 (MQ=38)
                                                                                             gCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAgg  >  2:3727875/1‑101 (MQ=255)
                                                                                                     |                                                                                              
CATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGG  >  NC_014034/1699847‑1700040

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG