Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 678,771 G→A 8.5% I1612I (ATC→ATT)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034678,7710G→A 8.5% 40.3 ‑0.2 ‑0.6233I1612I (ATC→ATT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (11/9):  ref base (111/102):  total (122/111)
Fisher's exact test strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.49e-01

Only 100 of 254 total aligned reads displayed.

CGCGGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCG  >  NC_014034/678671‑678871
                                                                                                    |                                                                                                    
cgcgGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCcgcg                                                                                                      <  1:3018171/101‑1 (MQ=255)
 gcgGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGa                                                                                                     >  1:67518/1‑101 (MQ=255)
   ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                   <  1:4266390/101‑1 (MQ=255)
    gATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCt                                                                                                  >  2:535449/1‑101 (MQ=255)
    gATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCt                                                                                                  >  1:1515712/1‑101 (MQ=255)
    gATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCt                                                                                                  <  1:3207642/101‑1 (MQ=255)
    gATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCt                                                                                                  >  1:1879346/1‑101 (MQ=255)
      tCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgc                                                                                                <  2:2359722/101‑1 (MQ=255)
      tCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgc                                                                                                <  2:480570/101‑1 (MQ=255)
        tGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgcgc                                                                                              >  2:1187694/1‑101 (MQ=255)
        tGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgcgc                                                                                              >  1:2198117/1‑101 (MQ=255)
        tGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgcgc                                                                                              >  1:1700419/1‑101 (MQ=255)
          cgcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCg                                                                                            >  2:519064/1‑101 (MQ=255)
          cgcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCg                                                                                            <  2:135624/101‑1 (MQ=255)
           gcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGt                                                                                           >  1:1370575/1‑101 (MQ=255)
            ccGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCAGCGATCTGCGCCGTg                                                                                          <  1:3635520/101‑1 (MQ=255)
                gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                      <  2:2406785/101‑1 (MQ=255)
                gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                      <  2:1771301/101‑1 (MQ=255)
                gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                      <  1:4539934/101‑1 (MQ=255)
                 cTCAGGCTTGGCATCAGCGCCGGCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAg                                                                                     <  1:3510550/101‑1 (MQ=255)
                  tCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgc                                                                                    <  1:3102086/101‑1 (MQ=255)
                    aGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                                                                  >  1:2683086/1‑101 (MQ=255)
                     ggCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                                                                 <  1:1759855/101‑1 (MQ=255)
                      gCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                                                                                >  2:1097078/1‑101 (MQ=255)
                       cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                               <  2:3337095/101‑1 (MQ=255)
                       cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                               >  2:2140305/1‑101 (MQ=255)
                          ggCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGAt                                                                            >  1:879537/1‑101 (MQ=255)
                          ggCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGAt                                                                            <  2:869383/101‑1 (MQ=255)
                             aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                         <  1:899555/101‑1 (MQ=255)
                             aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                         <  1:1147247/101‑1 (MQ=255)
                              tCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGc                                                                        >  2:458744/1‑101 (MQ=255)
                                 gcgcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTc                                                                     <  1:235137/101‑1 (MQ=255)
                                   gcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTcgc                                                                   >  1:1067975/1‑101 (MQ=255)
                                        cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                              <  1:3466893/101‑1 (MQ=255)
                                        cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                              >  1:933827/1‑101 (MQ=255)
                                        cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                              >  2:3742495/1‑101 (MQ=255)
                                        cccGACCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                              <  1:519112/101‑1 (MQ=255)
                                         ccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAg                                                             <  2:2197967/101‑1 (MQ=255)
                                         ccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAg                                                             >  2:3829685/1‑101 (MQ=255)
                                           gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                           >  1:355940/1‑101 (MQ=255)
                                            gtCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                          <  1:2388166/100‑1 (MQ=255)
                                             tCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGt                                                         >  2:2821398/1‑101 (MQ=255)
                                             tCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGt                                                         >  2:3252072/1‑101 (MQ=255)
                                             tCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGt                                                         <  1:1568126/101‑1 (MQ=255)
                                                ccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgcc                                                      >  1:3997205/1‑101 (MQ=255)
                                                 cGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccg                                                     <  2:1879220/101‑1 (MQ=255)
                                                 cGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccg                                                     >  2:3988100/1‑101 (MQ=255)
                                                 cGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccg                                                     <  2:158425/101‑1 (MQ=255)
                                                 cGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccg                                                     <  1:2892154/101‑1 (MQ=255)
                                                  gCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgc                                                    >  2:3674400/1‑101 (MQ=255)
                                                   ccATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgcc                                                   >  2:1954315/1‑101 (MQ=255)
                                                      tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                <  2:2075514/101‑1 (MQ=255)
                                                      tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                <  1:292847/101‑1 (MQ=255)
                                                      tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                >  2:3041345/1‑101 (MQ=255)
                                                      tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                <  2:3010976/101‑1 (MQ=255)
                                                      tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTAGCGCTCAGCCGTGCCGCCTgc                                                >  2:229863/1‑101 (MQ=255)
                                                       ggCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgcg                                               >  2:4064849/1‑101 (MQ=255)
                                                         ccGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGcc                                             <  2:3633053/101‑1 (MQ=255)
                                                          cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                            >  1:1976969/1‑101 (MQ=255)
                                                           gagaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGt                                           <  1:2768329/101‑1 (MQ=255)
                                                            agaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTg                                          <  1:3233292/101‑1 (MQ=255)
                                                              aCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGt                                        <  1:1201677/101‑1 (MQ=255)
                                                                cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                      >  1:2428800/1‑101 (MQ=255)
                                                                cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                      <  1:4650698/101‑1 (MQ=255)
                                                                cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                      <  1:3857900/101‑1 (MQ=255)
                                                                    aGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGACAgcgc                                  >  2:618886/1‑101 (MQ=255)
                                                                     gCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGcc                                 >  1:2081092/1‑101 (MQ=255)
                                                                      cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                >  2:2710580/1‑101 (MQ=255)
                                                                      cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCACCa                                >  1:333341/1‑101 (MQ=255)
                                                                          gCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGAt                            <  1:458781/101‑1 (MQ=255)
                                                                           cGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATg                           >  2:678698/1‑101 (MQ=255)
                                                                            gATATCCGCGGTTTCCAGCGCCGCGATCTGAGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                          >  1:4152820/1‑101 (MQ=255)
                                                                            gATATCCGCGGTGTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTTCGTCGCGCGCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                          <  2:795478/101‑1 (MQ=255)
                                                                            cagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                          >  1:329604/4‑101 (MQ=2)
                                                                            cagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                          >  2:3509874/4‑101 (MQ=2)
                                                                              tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                        <  1:2005976/101‑1 (MQ=255)
                                                                              gatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                        <  1:171829/100‑1 (MQ=2)
                                                                               atCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                       >  1:848824/1‑101 (MQ=2)
                                                                               atCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                       <  1:3501679/101‑1 (MQ=2)
                                                                                tcggcgGTTTCCAGCGCCGCACCCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                      <  1:380566/98‑2 (MQ=2)
                                                                                tcggcgGTTTCCAGCGCCCCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCTTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                      <  1:3312926/98‑2 (MQ=2)
                                                                                tCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAtt                      <  1:2805828/101‑1 (MQ=255)
                                                                                   gcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGt                   >  2:2878734/1‑101 (MQ=255)
                                                                                     ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                 >  1:3843862/1‑101 (MQ=255)
                                                                                      gTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGcc                <  1:229887/101‑1 (MQ=255)
                                                                                      gTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGcc                >  2:2637587/1‑101 (MQ=255)
                                                                                      gTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGGCAGCGCCACGATGTCATTCGTGcc                >  1:3027675/1‑101 (MQ=255)
                                                                                       tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc               <  1:4509094/101‑1 (MQ=255)
                                                                                       tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc               >  2:342602/1‑101 (MQ=255)
                                                                                       tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc               >  1:2586946/1‑101 (MQ=255)
                                                                                        ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa              >  1:2339620/1‑101 (MQ=255)
                                                                                        ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa              >  2:319666/1‑101 (MQ=255)
                                                                                          ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGc            <  1:42634/101‑1 (MQ=255)
                                                                                          ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGc            <  1:625575/101‑1 (MQ=255)
                                                                                           cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc           <  2:1700308/101‑1 (MQ=255)
                                                                                            aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc          >  2:3137656/1‑101 (MQ=255)
                                                                                             gcgcCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcccc         >  2:634970/1‑101 (MQ=255)
                                                                                                 cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa     >  2:3900663/1‑101 (MQ=255)
                                                                                                   cgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAgc   >  2:200611/1‑101 (MQ=255)
                                                                                                    gATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAgcg  >  2:2581959/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGCGGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCG  >  NC_014034/678671‑678871

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