Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,441 C→G 9.9% A722A (GCG→GCC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4410C→G 9.9% 30.8 ‑4.2 ‑0.0162A722A (GCG→GCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/13):  ref base (77/69):  total (80/82)
Fisher's exact test strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01

Only 100 of 189 total aligned reads displayed.

TCTGCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGG  >  NC_014034/681341‑681541
                                                                                                    |                                                                                                     
tCTGCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATcc                                                                                                       <  2:4375309/101‑1 (MQ=255)
  tGCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgc                                                                                                     <  1:2560773/101‑1 (MQ=255)
   gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                    >  2:2124931/1‑101 (MQ=255)
   gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                    >  1:2943067/1‑101 (MQ=255)
    ccGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgccg                                                                                                   <  1:4937373/101‑1 (MQ=255)
    ccGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgccg                                                                                                   <  1:613770/101‑1 (MQ=255)
     cgcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCcgccgc                                                                                                  <  1:3701221/101‑1 (MQ=255)
        gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                               >  2:18645/1‑101 (MQ=255)
        gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                               >  2:2888385/1‑101 (MQ=255)
         tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                              >  1:1452562/1‑101 (MQ=255)
           gACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAg                                                                                            >  1:3982941/1‑101 (MQ=255)
           gACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAg                                                                                            >  2:2352019/1‑101 (MQ=255)
             cAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgcg                                                                                          <  1:4784404/101‑1 (MQ=255)
              aGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgcag                                                                                         >  1:3115245/1‑99 (MQ=255)
               gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                        <  1:4393216/101‑1 (MQ=255)
                 gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                      <  1:2042673/101‑1 (MQ=255)
                 gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                      <  2:1171606/101‑1 (MQ=255)
                 ccctcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                      <  2:4501568/97‑1 (MQ=255)
                       aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                >  2:1129468/1‑101 (MQ=255)
                        tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                               <  2:2445812/101‑1 (MQ=255)
                         cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                              >  1:514676/1‑101 (MQ=255)
                         cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                              >  2:3050349/1‑101 (MQ=255)
                           gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                            <  2:3141744/101‑1 (MQ=255)
                           gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                            >  2:2928914/1‑101 (MQ=255)
                            cgcgcCTGAAGGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCCGATCCGCCGCCTCGAGCGGCGCGAGCTGCGTGc                                                                           >  1:4363226/1‑101 (MQ=255)
                            cgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGc                                                                           <  2:637081/101‑1 (MQ=255)
                             gcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                          <  1:2843335/101‑1 (MQ=255)
                                 cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                      <  2:2032147/101‑1 (MQ=255)
                                 cTGAAAGCCCGACCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                      >  1:1955285/1‑101 (MQ=255)
                                   gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                    <  2:1845928/101‑1 (MQ=255)
                                    aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                   >  2:472597/1‑101 (MQ=255)
                                     aaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcg                                                                  >  1:4505818/1‑101 (MQ=255)
                                      aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                 >  2:4487618/1‑101 (MQ=255)
                                      aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                 <  2:1368828/101‑1 (MQ=255)
                                      aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                 <  2:3048653/101‑1 (MQ=255)
                                          cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                             >  2:1532296/1‑101 (MQ=255)
                                            gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                           <  2:4790044/101‑1 (MQ=255)
                                            gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                           <  1:690547/101‑1 (MQ=255)
                                             ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                          <  2:4616803/101‑1 (MQ=255)
                                               aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                        >  1:3997497/1‑101 (MQ=255)
                                                tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgc                                                       >  2:2026532/1‑101 (MQ=255)
                                                tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgc                                                       <  2:2505273/101‑1 (MQ=255)
                                                  gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                     <  1:2085021/101‑1 (MQ=255)
                                                   ttgCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGTGCGCCGCGATCTGCCTGCTGGTGGGCGCGCCGATCTgcgcg                                                    >  2:932109/3‑100 (MQ=38)
                                                    ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                   >  1:4486761/1‑101 (MQ=255)
                                                       cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                >  2:866258/1‑101 (MQ=255)
                                                       cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                <  1:3872949/101‑1 (MQ=255)
                                                        gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                               <  1:547224/101‑1 (MQ=255)
                                                        gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                               <  2:2828017/101‑1 (MQ=255)
                                                          gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCa                                             <  2:3982639/101‑1 (MQ=255)
                                                           cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                            <  2:4143183/101‑1 (MQ=255)
                                                           cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                            >  1:3356525/1‑101 (MQ=255)
                                                           cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                            <  1:871152/101‑1 (MQ=255)
                                                            ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGa                                           >  1:1621578/1‑101 (MQ=255)
                                                                cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgc.g                                      >  1:4137795/1‑101 (MQ=255)
                                                                  cGGAGAGCTCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                    <  1:4246863/101‑1 (MQ=255)
                                                                  cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                    <  2:4839140/101‑1 (MQ=255)
                                                                  cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                    <  1:371856/101‑1 (MQ=255)
                                                                   ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gccg                                   >  1:127957/1‑101 (MQ=255)
                                                                    gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gccgc                                  >  1:4873147/1‑101 (MQ=255)
                                                                     agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.gccgcc                                 <  1:1019758/101‑1 (MQ=255)
                                                                      gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCt                                <  2:341714/101‑1 (MQ=255)
                                                                       agCGCCGTCGCGCTGATCTGCGCTAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTg                               <  2:1829619/101‑1 (MQ=255)
                                                                         cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTgcg                             >  1:1941853/1‑101 (MQ=255)
                                                                          gcCGTCGCGCTGATCTGCGCCAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTgcgc                            >  1:4490819/1‑101 (MQ=255)
                                                                           ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGcc                           <  1:4075698/101‑1 (MQ=255)
                                                                            cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcg                          >  2:4075235/1‑101 (MQ=255)
                                                                             gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgt                         >  2:48974/1‑101 (MQ=255)
                                                                             gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgt                         >  1:3677838/1‑101 (MQ=255)
                                                                             gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCATACGCcgt                         >  2:4007452/1‑101 (MQ=255)
                                                                               cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcg                       >  2:2406248/1‑101 (MQ=255)
                                                                                  gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCa                    <  1:3279146/101‑1 (MQ=255)
                                                                                   gggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAg                   <  1:166051/91‑1 (MQ=2)
                                                                                   cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAg                   >  2:3578378/1‑101 (MQ=255)
                                                                                   cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGACTGCGCCGTCGCCAg                   >  1:4445062/1‑101 (MQ=255)
                                                                                    ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAgc                  <  2:2056722/92‑1 (MQ=2)
                                                                                    ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAgc                  <  1:1093843/92‑1 (MQ=2)
                                                                                     gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAgcg                 >  2:4562526/1‑101 (MQ=255)
                                                                                        ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCa              <  1:5008976/96‑1 (MQ=2)
                                                                                         cgcgcCAGATCCCCCGCCCCGAGCGCCGCCATCTTCCTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCAc             <  2:5045812/100‑1 (MQ=25)
                                                                                         cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCAc             <  1:901341/97‑1 (MQ=2)
                                                                                          gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg            >  2:3791042/1‑101 (MQ=255)
                                                                                          gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg            <  2:4261747/101‑1 (MQ=255)
                                                                                          gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCTTGCTCGACAGGGCCAGAAGCTGCGTCGAGCCCAG.CGCGGCCGCCTGCGCCGTCGTCAGCCCCACg            >  2:1636690/4‑101 (MQ=2)
                                                                                          gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTGGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg            <  2:3281452/98‑1 (MQ=2)
                                                                                          gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg            <  1:3925425/98‑1 (MQ=2)
                                                                                          gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg            <  2:2598648/98‑1 (MQ=2)
                                                                                           caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGa           <  2:1836725/99‑1 (MQ=9)
                                                                                           caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGACGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGGGCCACGa           >  2:5021035/3‑101 (MQ=9)
                                                                                            gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt          <  1:18655/101‑1 (MQ=255)
                                                                                             ccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg         <  2:2788658/101‑1 (MQ=9)
                                                                                             ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg         >  1:3937935/1‑101 (MQ=255)
                                                                                              cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGt        <  1:3513890/101‑1 (MQ=255)
                                                                                                 aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt     >  1:2413151/1‑101 (MQ=255)
                                                                                                 aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt     <  2:3819267/101‑1 (MQ=255)
                                                                                                 aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAAGC.GCAGCCTGCGCCGTCGTCAGCGCCACGATGTCAt     >  2:2497699/1‑101 (MQ=255)
                                                                                                  tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc    <  1:4607355/101‑1 (MQ=255)
                                                                                                   ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg   <  1:4181143/101‑1 (MQ=255)
                                                                                                   ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg   <  1:3046268/101‑1 (MQ=255)
                                                                                                    cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg  >  2:599939/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
TCTGCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGG  >  NC_014034/681341‑681541

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