Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,666 G→C 11.2% L647L (CTC→GCC RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6660G→C11.2% 85.1 ‑1.5 ‑0.1320L647L (CTC→CTGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (21/15):  ref base (131/153):  total (152/168)
Fisher's exact test strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.28e-01

Only 100 of 334 total aligned reads displayed.

TGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCG  >  NC_014034/681566‑681766
                                                                                                    |                                                                                                    
tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                                                                                                      <  1:4662729/101‑1 (MQ=255)
    cAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcg                                                                                                  >  2:5577221/1‑101 (MQ=255)
      gcgcCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcg                                                                                                <  1:4918593/101‑1 (MQ=255)
       cgccgcCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcgc                                                                                               <  1:2623537/101‑1 (MQ=255)
          cgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGa                                                                                            >  2:3254359/1‑101 (MQ=255)
             cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                         <  2:1314277/101‑1 (MQ=255)
               ccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGa                                                                                       >  2:4803640/1‑101 (MQ=255)
                 gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                                     <  1:421261/101‑1 (MQ=255)
                 gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                                     >  2:1441944/1‑101 (MQ=255)
                  tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                                    <  1:2224891/101‑1 (MQ=255)
                   cGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCgg                                                                                   >  2:6563992/1‑101 (MQ=255)
                    gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                                                                  <  2:6452798/101‑1 (MQ=255)
                    gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                                                                  >  2:1824375/1‑101 (MQ=255)
                      ggtggcaaccgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGc                                                                                <  2:4008345/92‑1 (MQ=2)
                       gtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGcc                                                                               >  2:6344278/9‑100 (MQ=2)
                       gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCt                                                                               <  1:6635571/101‑1 (MQ=255)
                        aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc                                                                              <  1:4714928/101‑1 (MQ=255)
                        aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc                                                                              >  2:1723369/1‑101 (MQ=255)
                         ggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                             >  2:2527213/7‑101 (MQ=2)
                         ggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                             <  1:4470702/95‑1 (MQ=2)
                         ggcaaccgccggGATCTGCGCCGCGCTCAGCGCCGCGATCGGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                             <  2:4940359/95‑1 (MQ=2)
                         aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                             <  1:3483707/101‑1 (MQ=255)
                         aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                             <  2:6168609/101‑1 (MQ=255)
                           ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGc                                                                           >  2:5633643/1‑101 (MQ=255)
                           caatcgccgcGATCTGCGCCGCGCTCAGCCCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                           <  1:3875050/97‑2 (MQ=2)
                            cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                          <  1:4812057/101‑1 (MQ=255)
                            aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                          >  2:3366478/4‑99 (MQ=2)
                            aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                          <  1:3753234/98‑3 (MQ=2)
                             atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                                         <  1:1664634/99‑4 (MQ=2)
                             aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                         <  1:5853909/101‑1 (MQ=255)
                             aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                         >  2:4874798/1‑101 (MQ=255)
                               cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                                                       <  2:2113950/101‑1 (MQ=255)
                               cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                                                       >  2:3675793/1‑101 (MQ=255)
                                gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                                      >  2:5127375/1‑95 (MQ=2)
                                gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAGCCTggc                                                                      >  1:2893768/1‑99 (MQ=2)
                                 ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                                     <  1:7104338/101‑1 (MQ=255)
                                 ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                                     <  2:731995/101‑1 (MQ=255)
                                  cgcgATCTGCGCGGCGCTGAGCGCGGCGATCTGCCCGGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggcca                                                                    <  2:1146962/101‑9 (MQ=2)
                                     gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                 >  2:2380581/1‑101 (MQ=255)
                                     gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                 >  2:5905907/1‑101 (MQ=255)
                                      aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                >  1:10054/1‑101 (MQ=255)
                                      aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                >  2:1915670/1‑101 (MQ=255)
                                      aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                >  2:4980/1‑101 (MQ=255)
                                      aTCTGCGCCGTGCTCAGCGCCGCCGCCGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                <  1:1909691/101‑1 (MQ=255)
                                      aTCTCCGCCTTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                <  2:1432009/101‑1 (MQ=255)
                                       tCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGc                                                               >  1:3100605/1‑101 (MQ=255)
                                         tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                                                             <  2:6811279/101‑1 (MQ=255)
                                         tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                                                             >  2:2835930/1‑101 (MQ=255)
                                          gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                            <  1:5112569/101‑1 (MQ=255)
                                          gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                            <  1:981795/101‑1 (MQ=255)
                                          gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                            >  1:1319941/1‑101 (MQ=255)
                                          gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                            <  1:1110902/101‑1 (MQ=255)
                                              cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCAGATcggaagagcgt                                                        >  2:575541/1‑91 (MQ=255)
                                               gTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCt                                                       <  2:3386152/101‑1 (MQ=255)
                                                 gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                                     >  1:2366850/1‑101 (MQ=255)
                                                  cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                                    >  2:4136435/1‑101 (MQ=255)
                                                    cAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGc                                                  <  2:3846553/101‑1 (MQ=255)
                                                       cgccgccgcCTGCCCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                               <  2:3594422/101‑1 (MQ=255)
                                                       cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                               <  2:3085901/101‑1 (MQ=255)
                                                          cgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGaaa                                            <  2:2879745/101‑1 (MQ=255)
                                                            ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                          <  1:4544833/101‑1 (MQ=255)
                                                              gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                        <  1:5906210/101‑1 (MQ=255)
                                                              cccTGCCCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                        <  2:5449923/100‑1 (MQ=255)
                                                               ccggCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgcc                                       <  1:6884782/98‑1 (MQ=255)
                                                                 tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                     >  2:4561230/1‑101 (MQ=255)
                                                                 tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                     >  2:6047679/1‑101 (MQ=255)
                                                                 tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                     <  1:1723433/101‑1 (MQ=255)
                                                                      cGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGc                                >  2:861678/1‑101 (MQ=255)
                                                                      cGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGc                                <  1:5597120/101‑1 (MQ=255)
                                                                        aGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTTAGGCCCGAAACCGCCGCGCTGccc                              <  1:3649564/101‑1 (MQ=255)
                                                                         gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                             <  2:4319250/101‑1 (MQ=255)
                                                                          ccAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAg                            >  1:6769038/1‑101 (MQ=255)
                                                                           cAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGc                           <  1:1441995/101‑1 (MQ=255)
                                                                             aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc                         <  1:1824439/101‑1 (MQ=255)
                                                                             aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc                         >  2:5477530/1‑101 (MQ=255)
                                                                              gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                        <  2:6245714/101‑1 (MQ=255)
                                                                               cccGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCcg                       >  2:6754521/1‑101 (MQ=255)
                                                                                ccGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCcgc                      >  1:24551/1‑101 (MQ=255)
                                                                                  ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa                    >  1:5227815/1‑101 (MQ=255)
                                                                                   gCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGAt                   >  1:2103962/1‑101 (MQ=255)
                                                                                     cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt                 >  1:6037560/1‑101 (MQ=255)
                                                                                       ccGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgc               <  1:4336338/101‑1 (MQ=255)
                                                                                        cGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcg              <  2:67185/101‑1 (MQ=255)
                                                                                        cGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcg              >  2:1646848/1‑101 (MQ=255)
                                                                                          gaaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGCCCGGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc            <  2:2041667/100‑1 (MQ=255)
                                                                                          cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGGCGCGCTAAGGCCCGAAACCGCCGGCCTGCCCAGCCCCGCGCTCTgcccg            >  1:5702557/1‑98 (MQ=255)
                                                                                           aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg           >  2:3011580/1‑101 (MQ=255)
                                                                                             ggCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTg         >  2:2730670/1‑101 (MQ=255)
                                                                                               cTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGcc       >  2:2252478/1‑101 (MQ=255)
                                                                                               cTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGcc       >  1:652753/1‑101 (MQ=255)
                                                                                               cTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGcc       <  2:4476271/101‑1 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      <  2:347675/101‑1 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      <  1:217664/101‑1 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      <  1:5017887/101‑1 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      >  2:5281199/1‑101 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      <  2:5043282/101‑1 (MQ=255)
                                                                                                tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc      >  1:6346823/1‑101 (MQ=255)
                                                                                                 ggAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa     >  1:4140720/1‑101 (MQ=255)
                                                                                                 ggAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa     <  1:5309630/101‑1 (MQ=255)
                                                                                                    gagCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCg  <  1:3254475/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCG  >  NC_014034/681566‑681766

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