Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,129 A→G 35.4% D826D (GAT→GAC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,1290A→G35.4% 78.0 ‑0.4 ‑0.082D826D (GAT→GAC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (18/11):  ref base (30/23):  total (48/34)
Fisher's exact test strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.48e-01

TCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCG  >  NC_014034/681034‑681229
                                                                                               |                                                                                                    
tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                                 <  1:577544/101‑1 (MQ=255)
tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                                 <  2:2130884/101‑1 (MQ=255)
tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                                 >  1:2070060/1‑101 (MQ=255)
tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGCCAGCGCGACGATGTCGCCGGt                                                                                                 >  2:3752015/1‑101 (MQ=255)
 cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                                <  2:766218/101‑2 (MQ=255)
 cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                                <  1:2361815/101‑2 (MQ=255)
  gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                               <  1:563625/101‑3 (MQ=255)
    cctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgac                                                                                             <  1:2013672/101‑5 (MQ=255)
     ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCGATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                            <  2:1512624/101‑6 (MQ=38)
     ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                            >  1:3479584/1‑96 (MQ=255)
     ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                            <  1:668132/101‑6 (MQ=255)
     ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                            >  2:447106/1‑96 (MQ=255)
      tcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacag                                                                                           <  1:1552349/101‑7 (MQ=38)
       cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                                          >  1:952811/1‑94 (MQ=38)
       cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                                          <  1:4589884/101‑8 (MQ=38)
       cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                                          >  1:1000865/1‑94 (MQ=38)
        gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                                         >  1:2311853/1‑93 (MQ=37)
         aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                                        >  2:4377718/1‑92 (MQ=37)
                   gCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcg                                                                              >  2:1208874/1‑100 (MQ=25)
                    cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGTCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                             >  2:1793744/1‑99 (MQ=25)
                    cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                             >  2:4287412/1‑99 (MQ=25)
                     tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                            >  1:2917857/1‑98 (MQ=25)
                      gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                                           >  2:1730515/1‑97 (MQ=25)
                      gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                                           >  1:3712510/1‑97 (MQ=25)
                        gcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                                         >  1:1237129/1‑95 (MQ=25)
                        cccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                                         <  1:829920/100‑7 (MQ=21)
                          cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGGCGGACAGCCCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccaccgg                                                                       >  2:3113376/1‑93 (MQ=16)
                          cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactgg                                                                       >  1:4219663/1‑93 (MQ=21)
                           gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                                      >  1:1619573/1‑92 (MQ=21)
                            tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                                     >  2:1414237/1‑91 (MQ=255)
                            tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                                     >  1:4254928/1‑91 (MQ=25)
                                                                    ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGTCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                             >  1:3921526/11‑101 (MQ=37)
                                                                    ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                             >  2:1736550/11‑101 (MQ=37)
                                                                    ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                             >  2:2361205/11‑101 (MQ=37)
                                                                     tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                            >  2:2230505/10‑101 (MQ=38)
                                                                      gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc                           <  2:3294009/93‑1 (MQ=255)
                                                                      gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc                           <  1:1556247/93‑1 (MQ=255)
                                                                      gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCCAGCCCCTTGATCGCATTGGCACCCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc                           <  1:3987090/93‑1 (MQ=37)
                                                                       cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg                          <  1:4575800/94‑1 (MQ=255)
                                                                       cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg                          >  2:4181603/8‑101 (MQ=255)
                                                                       cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg                          >  1:1892995/8‑101 (MQ=255)
                                                                       cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg                          >  1:1154046/8‑101 (MQ=255)
                                                                        gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt                         <  1:1613432/95‑1 (MQ=255)
                                                                         ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg                        <  1:2105024/96‑1 (MQ=255)
                                                                         ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCAATGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg                        >  2:4037655/6‑101 (MQ=38)
                                                                          cgtcgttaGCGCCACGATGTCATCGGTGCCGATGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc                       <  2:2776385/94‑1 (MQ=38)
                                                                          cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc                       <  1:3260145/97‑1 (MQ=255)
                                                                           gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt                      >  1:768814/4‑101 (MQ=255)
                                                                            tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                     >  1:2245014/3‑101 (MQ=255)
                                                                            tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                     <  2:3443084/99‑1 (MQ=255)
                                                                            tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                     <  1:2048137/99‑1 (MQ=255)
                                                                            tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                     <  1:193633/99‑1 (MQ=255)
                                                                            tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                     <  2:2264202/99‑1 (MQ=255)
                                                                             cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa                    <  1:3194457/100‑1 (MQ=255)
                                                                             cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa                    <  1:2330521/100‑1 (MQ=255)
                                                                              gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg                   >  1:1785196/1‑101 (MQ=255)
                                                                              ctCAGCGCCACGATGTCATCGGTGCCGAGGGCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg                   <  2:3315582/100‑1 (MQ=255)
                                                                               tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc                  >  2:245392/1‑101 (MQ=255)
                                                                               tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc                  >  1:2285044/1‑101 (MQ=255)
                                                                                cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg                 >  1:2296602/1‑101 (MQ=255)
                                                                                 aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc                <  2:2748494/101‑1 (MQ=255)
                                                                                  gcgcCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg               <  2:3760670/101‑1 (MQ=255)
                                                                                   cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg              <  1:2384733/101‑1 (MQ=255)
                                                                                    gccacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGc             <  2:657350/98‑1 (MQ=255)
                                                                                    gccacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGc             >  2:3595334/4‑101 (MQ=255)
                                                                                     ccacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGcc            <  2:1206446/99‑1 (MQ=255)
                                                                                     ccacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGcc            <  1:833544/99‑1 (MQ=255)
                                                                                      cacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCa           >  1:2652255/2‑101 (MQ=255)
                                                                                      cacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCa           >  1:4408468/2‑101 (MQ=255)
                                                                                       acgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCAt          >  2:3529150/1‑101 (MQ=255)
                                                                                       acgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGCGCTCAGCGCGGCCAt          <  1:3932219/101‑1 (MQ=255)
                                                                                         gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt        >  1:2125816/1‑101 (MQ=255)
                                                                                         gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt        >  2:2959579/1‑101 (MQ=255)
                                                                                         gaTGTCATCGGTGCCGAGGCCATTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt        >  2:818859/1‑101 (MQ=255)
                                                                                          aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg       >  2:3761644/1‑101 (MQ=255)
                                                                                          aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg       <  2:848748/101‑1 (MQ=255)
                                                                                          aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGAGCGGCCATCTg       >  2:569339/1‑101 (MQ=255)
                                                                                           tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgc      <  1:716051/101‑1 (MQ=255)
                                                                                            gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg     >  2:3662541/1‑101 (MQ=255)
                                                                                            gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg     >  1:4211027/1‑101 (MQ=255)
                                                                                            gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg     >  2:4067765/1‑101 (MQ=255)
                                                                                            gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg     >  2:4419048/1‑101 (MQ=255)
                                                                                            gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg     >  2:1790010/1‑101 (MQ=255)
                                                                                              cATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGcc   >  2:3315253/1‑101 (MQ=255)
                                                                                              cATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGcc   >  2:1094543/1‑101 (MQ=255)
                                                                                               aTCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCg  <  1:4121699/101‑1 (MQ=255)
                                                                                               aTCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCg  <  1:2410835/101‑1 (MQ=255)
                                                                                               aTCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCg  <  2:892700/101‑1 (MQ=255)
                                                                                               |                                                                                                    
TCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCG  >  NC_014034/681034‑681229

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