Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,582 C→G 3.8% T675T (ACG→ACC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5820C→G 3.8% 6.0 ‑0.4 ‑0.0159T675T (ACG→ACC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (98/55):  total (103/56)
Fisher's exact test strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.52e-01

Only 100 of 173 total aligned reads displayed.

GATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTG  >  NC_014034/681483‑681680
                                                                                                   |                                                                                                  
gagcTGGGCCGAGGCCAGACGCCCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAACCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                                                                                                   <  2:2964576/98‑1 (MQ=25)
gATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                                                                                                   >  1:2919565/1‑101 (MQ=255)
 aTCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt                                                                                                  >  2:2137350/1‑101 (MQ=255)
 aTCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt                                                                                                  >  1:1812013/1‑101 (MQ=255)
 aTCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt                                                                                                  >  2:4725939/1‑101 (MQ=255)
  tCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTc                                                                                                 >  2:1060705/1‑101 (MQ=255)
    tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCCCCGTCGc                                                                                               >  1:2002993/1‑101 (MQ=255)
    tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGc                                                                                               >  1:4283666/1‑101 (MQ=255)
       gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa                                                                                            <  2:3870689/101‑1 (MQ=255)
        ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                                                           <  1:1854705/101‑1 (MQ=255)
          gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc                                                                                         >  1:1537200/1‑101 (MQ=255)
          gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc                                                                                         >  1:2763104/1‑101 (MQ=255)
           aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGTCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                                                        <  1:4304949/101‑1 (MQ=255)
           aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                                                        <  2:670230/101‑1 (MQ=255)
              ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                                     <  2:573376/101‑1 (MQ=255)
              ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                                     >  1:1931582/1‑101 (MQ=255)
                 gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                                                  >  1:4064672/1‑101 (MQ=255)
                  ccgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccg                                                                                 <  1:3975072/100‑1 (MQ=255)
                  aCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGCATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccg                                                                                 <  2:4362909/101‑1 (MQ=255)
                    gcgcCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCcgcg                                                                               <  2:1941716/101‑1 (MQ=255)
                     cgccgcCTGCGCCGTCGTCAGTGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgcaa                                                                              >  1:3723404/1‑99 (MQ=25)
                      gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                                             >  1:143154/1‑101 (MQ=255)
                      gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                                             <  2:2924538/101‑1 (MQ=255)
                      gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                                             <  1:3578648/101‑1 (MQ=255)
                      gcccccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgcaag                                                                             <  1:2498758/101‑4 (MQ=25)
                       ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                                            <  1:2825639/101‑1 (MQ=255)
                         gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTg                                                                          >  2:2539822/1‑101 (MQ=37)
                         gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                                          <  1:4070964/101‑1 (MQ=255)
                         gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                                          <  1:1402954/101‑1 (MQ=255)
                          ccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTgc                                                                         >  1:3818739/1‑101 (MQ=37)
                          ccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgc                                                                         >  2:1751623/1‑101 (MQ=255)
                          ccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgc                                                                         >  2:3338086/1‑101 (MQ=255)
                            tGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgagc                                                                       >  2:970050/1‑101 (MQ=255)
                              ctcggtggtcATCGCCACGTTGTACTGGGTGCGAAGCCCCCGGTTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCg                                                                     <  2:4444743/94‑1 (MQ=21)
                               gcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGGCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGt                                                                    >  1:4271667/1‑101 (MQ=255)
                                ccGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTg                                                                   >  1:1952956/1‑101 (MQ=255)
                                 cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                                                  <  1:272077/101‑1 (MQ=255)
                                 cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                                                  >  2:1595729/1‑101 (MQ=255)
                                   tcgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTc                                                                <  1:1511117/101‑1 (MQ=255)
                                    cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                                               >  2:3719622/1‑101 (MQ=255)
                                    cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                                               >  1:967349/1‑101 (MQ=255)
                                      tcAGCGCCACGATTTCATCGGTGCAAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgc                                                             <  2:2899497/101‑1 (MQ=255)
                                      tcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgc                                                             >  2:4622481/1‑101 (MQ=255)
                                       cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                            >  1:1469820/1‑101 (MQ=255)
                                         gcgcCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                                          >  2:1323768/1‑101 (MQ=255)
                                          caccACGATGTCATCGGTGCCAAAGCCGCGTATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCGGCGCCGTGCTCAGcgccg                                                         <  2:3856709/99‑1 (MQ=255)
                                            ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                       >  2:3543454/1‑101 (MQ=255)
                                            ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                       <  1:1169443/101‑1 (MQ=255)
                                                gATGTGAGCGGTGCCAAGGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                   <  2:3326284/101‑1 (MQ=255)
                                                gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                   >  2:3265453/1‑101 (MQ=255)
                                                  tGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGc                                                 <  1:2406422/101‑1 (MQ=255)
                                                    tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                               <  1:600018/101‑1 (MQ=255)
                                                    tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                               >  1:2413334/1‑101 (MQ=255)
                                                    tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                               >  2:1243477/1‑101 (MQ=255)
                                                     cATCGGTGCCAAAGCCGCAGATCGCCGCCGTGCCCAGCACCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGACGCCTGCAcc                                              <  1:771317/101‑1 (MQ=255)
                                                      aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                             >  1:4146562/1‑101 (MQ=255)
                                                      aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                             >  1:1688751/1‑101 (MQ=255)
                                                       tCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGa                                            >  1:1866995/1‑101 (MQ=255)
                                                         ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                          >  1:2008928/1‑101 (MQ=255)
                                                          gTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGcc                                         >  1:950041/1‑101 (MQ=255)
                                                           tGCTAAAGCCGCGGGTCGCGGCCGTGCCCCGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCCACGCCCGCACCGAGccc                                        >  2:3227701/1‑100 (MQ=37)
                                                              cAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGc                                     >  2:851872/1‑101 (MQ=255)
                                                              cAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGc                                     <  1:4562847/101‑1 (MQ=255)
                                                              cAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGc                                     <  2:1931445/101‑1 (MQ=255)
                                                               aaaGCCGCGGATCGCCGCCGTGCCCAGCTCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGcc                                    >  1:4574784/1‑101 (MQ=255)
                                                               aaaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGcc                                    >  2:2081491/1‑101 (MQ=255)
                                                                aaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGccc                                   >  2:1646207/1‑101 (MQ=255)
                                                                aaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGccc                                   <  2:1621446/101‑1 (MQ=255)
                                                                 aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg                                  >  2:4981791/1‑101 (MQ=255)
                                                                  gCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCgg                                 >  1:3812830/1‑101 (MQ=255)
                                                                    cgcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc                               <  2:1537070/101‑1 (MQ=255)
                                                                      cgGACCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAc                             <  2:1296089/101‑1 (MQ=255)
                                                                       ggATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAcc                            >  1:2121663/1‑101 (MQ=255)
                                                                       ggATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAcc                            <  1:813354/101‑1 (MQ=255)
                                                                        gATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCg                           >  2:188535/1‑101 (MQ=255)
                                                                        gATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCg                           >  1:1625579/1‑101 (MQ=255)
                                                                         accgccgcCGTCCCCAGCGCCGCCGCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGc                          <  2:2810861/99‑1 (MQ=255)
                                                                         aTCGCCGCCGTGCCGAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGc                          >  2:3124419/1‑101 (MQ=255)
                                                                         aTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAcccc                          >  2:3455190/1‑99 (MQ=255)
                                                                         aTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGc                          >  1:4011009/1‑101 (MQ=255)
                                                                           cgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa                        >  1:3438270/1‑101 (MQ=255)
                                                                            gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCTCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                       >  1:3518039/1‑101 (MQ=255)
                                                                            gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                       >  2:1070564/1‑101 (MQ=255)
                                                                            cccgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                       <  2:4750929/100‑1 (MQ=255)
                                                                             ccgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAgg                      >  1:303004/1‑101 (MQ=255)
                                                                               gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGcg                    >  2:4238439/1‑100 (MQ=255)
                                                                                 cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                  <  1:859597/101‑1 (MQ=255)
                                                                                 cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                  >  1:3864831/1‑101 (MQ=255)
                                                                                   tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                <  2:4505498/101‑1 (MQ=255)
                                                                                   tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                >  1:2677099/1‑101 (MQ=255)
                                                                                   ggCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                <  2:3937643/100‑1 (MQ=255)
                                                                                    gCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgaga               >  2:513929/1‑101 (MQ=255)
                                                                                    gCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgaga               >  1:4891768/1‑101 (MQ=255)
                                                                                      ccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc             >  1:2340437/1‑101 (MQ=255)
                                                                                        aGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgc           >  1:1242343/1‑101 (MQ=255)
                                                                                            ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg       <  2:2453129/101‑1 (MQ=255)
                                                                                              gccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc     >  1:3124973/1‑101 (MQ=255)
                                                                                               ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc    <  1:4514104/101‑1 (MQ=255)
                                                                                               ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc    >  1:4820271/1‑101 (MQ=255)
                                                                                                 cccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg  <  1:5019358/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTG  >  NC_014034/681483‑681680

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