Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,550 C→A 8.1% G353W (GGG→TGG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034682,5500C→A 8.1% 31.3 ‑0.2 ‑0.6185G353W (GGG→TGG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (6/9):  ref base (77/93):  total (83/102)
Fisher's exact test strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.77e-01

Only 100 of 198 total aligned reads displayed.

GCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAA  >  NC_014034/682450‑682650
                                                                                                    |                                                                                                     
gcCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGcccc                                                                                                       <  2:2850464/101‑1 (MQ=255)
 ccGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCg                                                                                                      >  1:4210937/1‑101 (MQ=255)
 ccGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCg                                                                                                      <  1:623598/101‑1 (MQ=255)
  cGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCgg                                                                                                     >  2:897385/1‑101 (MQ=255)
  cGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCgg                                                                                                     >  1:3817317/1‑101 (MQ=255)
  cGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCgg                                                                                                     <  2:1689232/101‑1 (MQ=255)
   gAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGt                                                                                                    <  1:4673577/101‑1 (MQ=255)
    aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                                   <  2:4125258/101‑1 (MQ=255)
    aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                                   <  1:197716/101‑1 (MQ=255)
    aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                                   <  1:2179143/101‑1 (MQ=255)
       tCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAgc                                                                                                <  1:951539/101‑1 (MQ=255)
        cAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAgcg                                                                                               >  2:3052743/1‑101 (MQ=255)
        cAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAgcg                                                                                               <  1:1329168/101‑1 (MQ=255)
          gcgcCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGcgcc                                                                                             <  2:4324308/101‑1 (MQ=255)
             ccTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCgccgcc                                                                                          >  2:3921250/1‑101 (MQ=255)
             ccTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCgccgcc                                                                                          >  1:2072036/1‑101 (MQ=255)
               tCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAg                                                                                        >  1:184456/1‑101 (MQ=255)
               tCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCCCCGCCAg                                                                                        >  1:1224285/1‑101 (MQ=255)
                cggcggCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGt                                                                                       >  2:629815/1‑101 (MQ=255)
                cggcggCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGt                                                                                       >  2:2053254/1‑101 (MQ=255)
                 ggcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGtt                                                                                      >  2:3115832/1‑101 (MQ=255)
                  gcggcCTCCTCGGACAGCTCCGCGATCTGCNGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                     <  1:1386888/101‑1 (MQ=255)
                  gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                     <  2:4819430/101‑1 (MQ=255)
                  gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                     <  2:2740670/101‑1 (MQ=255)
                     gcCTCCGCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGAt                                                                                  <  1:798382/101‑1 (MQ=255)
                       ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGATCCCCCCGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCg                                                                                <  2:1224168/101‑1 (MQ=255)
                       ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCg                                                                                <  1:556014/101‑1 (MQ=255)
                           tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                           >  1:1190011/1‑101 (MQ=255)
                            cGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTc                                                                          >  2:1110278/1‑101 (MQ=255)
                             ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                         >  1:3147225/1‑101 (MQ=255)
                             ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                         >  1:4394295/1‑101 (MQ=255)
                              gACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAg                                                                        <  2:1392703/101‑1 (MQ=255)
                                 aGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGccc                                                                     <  1:2995205/101‑1 (MQ=255)
                                  gcgcCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCg                                                                    <  1:1317446/101‑1 (MQ=255)
                                    gccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGc                                                                  >  2:2639325/1‑101 (MQ=255)
                                     ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                 <  1:1825308/101‑1 (MQ=255)
                                      cgcgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCa                                                                <  1:3295986/101‑1 (MQ=255)
                                      cgcgATCTGCCGCGTCCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCa                                                                <  2:2881080/101‑1 (MQ=255)
                                        cgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATg                                                              <  2:4025890/101‑1 (MQ=255)
                                          aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                            <  1:1604025/101‑1 (MQ=255)
                                          aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                            <  1:2564807/101‑1 (MQ=255)
                                           tCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTcc                                                           <  1:460444/101‑1 (MQ=255)
                                             tGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.GGCTCAGCCCGTCCATGTCCTc                                                         >  1:1589040/1‑101 (MQ=255)
                                              gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                        <  2:936180/101‑1 (MQ=255)
                                              gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                        <  1:4161895/101‑1 (MQ=255)
                                                 gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                                     <  2:1432816/101‑1 (MQ=255)
                                                   gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                   >  2:5000619/1‑101 (MQ=255)
                                                   gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                   <  1:2659120/101‑1 (MQ=255)
                                                   gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                   >  2:2163579/1‑101 (MQ=255)
                                                   gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                   >  2:3336365/1‑101 (MQ=255)
                                                   gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                   <  1:3336628/101‑1 (MQ=255)
                                                    tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                  <  2:184440/101‑1 (MQ=255)
                                                    tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                  >  1:1021938/1‑101 (MQ=255)
                                                     gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                 >  2:1345404/1‑101 (MQ=255)
                                                     gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                 >  2:473556/1‑101 (MQ=255)
                                                     gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                 >  1:713118/1‑101 (MQ=255)
                                                     gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                 <  1:2765948/101‑1 (MQ=255)
                                                      cccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAg                                                <  1:2615838/101‑1 (MQ=255)
                                                         aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                             >  1:2839268/1‑101 (MQ=255)
                                                         aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                             <  1:32098/101‑1 (MQ=255)
                                                          tcgccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGcccc                                            <  1:3737253/95‑1 (MQ=37)
                                                           accccccAGCTGTGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCTCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCt                                           <  1:2483253/99‑1 (MQ=255)
                                                           aaCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCt                                           <  2:3633457/101‑1 (MQ=255)
                                                           aaCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCt                                           <  1:3658307/101‑1 (MQ=255)
                                                             ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                         <  1:3083497/98‑2 (MQ=37)
                                                              cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                                        >  2:3771201/3‑99 (MQ=37)
                                                              ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                        <  2:4642083/101‑1 (MQ=255)
                                                              ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                        <  2:4902726/101‑1 (MQ=255)
                                                               gccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAg                                       <  2:3819934/100‑1 (MQ=25)
                                                               cccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAg                                       <  2:2485230/101‑1 (MQ=255)
                                                                ccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                      <  2:2992515/101‑1 (MQ=255)
                                                                ccACCTGGGCCGTGCTGCGCCCCGCCATCTGCGCCCCGGGCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                      <  1:1421511/101‑1 (MQ=255)
                                                                 cAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcg                                     >  1:853404/1‑101 (MQ=25)
                                                                  aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                    >  1:2389598/1‑101 (MQ=255)
                                                                  aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcgc                                    <  2:2466930/101‑1 (MQ=25)
                                                                     tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                 >  2:4359047/1‑101 (MQ=255)
                                                                     tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                 >  2:4221088/1‑101 (MQ=255)
                                                                      gggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggcg                                <  2:738270/101‑1 (MQ=255)
                                                                      ggaccGCGCTTAGCGGCGTCCCCTACGACCCGGTGAGCGACGCCAGTTCGATGG.TTCTCACCCCGGCCATGTCTTTGGTCGAAAGCCCCTTCAGCgcggcg                                <  1:1062642/98‑1 (MQ=255)
                                                                        gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                              <  2:2368605/101‑6 (MQ=21)
                                                                         ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                             >  2:2823665/1‑101 (MQ=255)
                                                                         ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                             <  2:2592166/101‑1 (MQ=255)
                                                                          cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccactgg                            <  2:1053054/101‑8 (MQ=25)
                                                                           gTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGcc                           <  2:3169182/101‑1 (MQ=255)
                                                                            tGCTCAGCGCCGCCATCTTCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                          >  2:3738516/1‑101 (MQ=255)
                                                                               tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                       <  1:1774271/101‑1 (MQ=255)
                                                                               tCAGCGCCGCCATCTGCCCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                       <  2:1588908/101‑1 (MQ=255)
                                                                                  gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc                    <  2:4579640/101‑1 (MQ=255)
                                                                                  gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc                    <  2:3870579/101‑1 (MQ=255)
                                                                                     ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                 >  1:4183555/1‑101 (MQ=255)
                                                                                     ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                 <  1:3424451/101‑1 (MQ=255)
                                                                                      cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaa                >  2:1482470/1‑101 (MQ=255)
                                                                                       gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa               <  2:4210620/101‑1 (MQ=255)
                                                                                           tCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGcc           <  2:4470436/101‑1 (MQ=255)
                                                                                            cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg          >  1:4852117/1‑101 (MQ=255)
                                                                                            cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg          >  1:3456272/1‑101 (MQ=255)
                                                                                             tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc         <  1:3589035/101‑1 (MQ=255)
                                                                                                 ccccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCg     <  1:475025/101‑1 (MQ=255)
                                                                                                   ccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGa   >  1:4354544/1‑101 (MQ=255)
                                                                                                   ccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGa   <  2:543005/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
GCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAA  >  NC_014034/682450‑682650

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG