Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,144 +T 1.9% coding (26/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,1441.→T 1.9% 8.1 ‑3.2 ‑0.1213coding (26/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (116/93):  total (116/97)
Fisher's exact test strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.50e-01

Only 100 of 229 total aligned reads displayed.

GCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCC  >  NC_014034/1396044‑1396243
                                                                                                                                                                                                        
gcgGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTgg                                                                                                     <  1:3829247/101‑1 (MQ=255)
 cgGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGc                                                                                                    <  1:1898545/101‑1 (MQ=255)
   ggggCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAg                                                                                                  >  2:1642514/1‑101 (MQ=255)
   ggggCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAg                                                                                                  <  1:668138/101‑1 (MQ=255)
    gggCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAgg                                                                                                 <  1:92064/101‑1 (MQ=255)
     ggCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGc                                                                                                >  2:4217669/1‑101 (MQ=255)
       cGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCaccctc                                                                                              <  2:3583439/101‑6 (MQ=255)
        gTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTGGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTcc                                                                                             >  1:3583034/1‑100 (MQ=255)
        gTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTcg                                                                                             >  2:1368982/1‑101 (MQ=255)
         tGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTcgc                                                                                            >  2:1539222/1‑101 (MQ=255)
          ggATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTcgcg                                                                                           >  1:1993935/1‑101 (MQ=255)
          ggATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTcgcg                                                                                           <  2:1653784/101‑1 (MQ=255)
             tCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGcgccg                                                                                        >  1:1841090/1‑101 (MQ=255)
             tCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGcgccg                                                                                        >  2:2630928/1‑101 (MQ=255)
               gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCgccgcc                                                                                      <  2:362169/101‑1 (MQ=255)
               gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCcccccc                                                                                      >  1:3189825/1‑101 (MQ=255)
                aaaCAGGCGCATCTCATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCgccgccc                                                                                     >  1:3355976/1‑100 (MQ=255)
                 aaCAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCAc                                                                                    >  1:685672/1‑101 (MQ=255)
                   cAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGc                                                                                  >  1:2881622/1‑101 (MQ=255)
                       cgcATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCg                                                                              >  2:516757/1‑101 (MQ=255)
                            cTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCCCGCTGCGTcccc                                                                         >  2:128131/1‑98 (MQ=255)
                             tGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTc                                                                        >  2:2904870/1‑101 (MQ=255)
                              gATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTcc                                                                       >  2:2139666/1‑101 (MQ=255)
                                 gTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTccccc                                                                    <  1:425058/101‑1 (MQ=255)
                                 gTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTccccc                                                                    <  2:448397/101‑1 (MQ=255)
                                  tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCCCGCTGCGTCGCTcccccc                                                                   >  2:4085451/1‑101 (MQ=255)
                                  tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                   <  1:115749/101‑1 (MQ=255)
                                  tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                   <  1:1675210/101‑1 (MQ=255)
                                  tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                   >  1:2255804/1‑101 (MQ=255)
                                   atatGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACCCTGCGTCGCTccccccc                                                                  >  2:2305155/1‑100 (MQ=255)
                                     atGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCcgcg                                                                >  2:3533939/1‑101 (MQ=255)
                                     atGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCCCGCCCGCACGCTGCGTCGCCTCCCCcgcg                                                                >  2:210750/1‑101 (MQ=37)
                                       gCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCggg                                                              >  2:4164026/1‑101 (MQ=255)
                                             tGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGcaca                                                        <  2:1619838/101‑4 (MQ=255)
                                             tGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTg                                                        >  2:2316166/1‑101 (MQ=255)
                                             tGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTg                                                        <  2:3289075/101‑1 (MQ=255)
                                             tGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTg                                                        <  1:1086573/101‑1 (MQ=255)
                                              gagaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTgg                                                       >  2:1086479/1‑101 (MQ=255)
                                              gagaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTgg                                                       <  1:2714707/101‑1 (MQ=255)
                                                gaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                     <  2:2075428/101‑1 (MQ=255)
                                                gaCGGCTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                     <  2:819677/101‑1 (MQ=255)
                                                 aCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCCGCGTCGCCCCCCCCCCGGGGGCGGGGGTc                                                    >  1:3007248/1‑101 (MQ=37)
                                                   ggCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTTGGGCGCTCCCAACGCGGGGGCGGGGGTCGa                                                  >  2:709024/1‑101 (MQ=25)
                                                   ggCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCCCGGGGGCGGGGGTCgg                                                  >  1:3259775/1‑100 (MQ=255)
                                                   ggCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCCCGCCCCCACGCCGCGTCGCTCCCCCCGCGGGGGCCTGGATCGa                                                  >  1:2105006/1‑101 (MQ=37)
                                                    gCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                 >  1:339317/1‑101 (MQ=255)
                                                    gCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                 <  1:1396583/101‑1 (MQ=255)
                                                    gCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                 <  1:2936005/101‑1 (MQ=255)
                                                     cTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                >  2:4053027/1‑101 (MQ=255)
                                                       gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                              <  1:974321/101‑1 (MQ=255)
                                                           cTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGGGCCGCCACGCTGCGGCGCTCCCCCCCCGGGGGCGTGGATCGAAACCCCGc                                          >  1:4160433/1‑101 (MQ=38)
                                                           cTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCCCGCTGCGTCGCCCCCCCCCCGGGGGCGCGGATCGAAACCCGGc                                          >  1:2504507/1‑101 (MQ=38)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGGTCCAAACCCCGCg                                         >  2:2098809/1‑101 (MQ=255)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCg                                         >  1:3061192/1‑101 (MQ=255)
                                                             ccccccGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGa                                        <  1:1290191/101‑1 (MQ=255)
                                                             ccccccGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGa                                        >  2:1176414/1‑101 (MQ=255)
                                                               ccccGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGaaa                                      <  1:1768591/101‑1 (MQ=255)
                                                               ccccGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCCCGGGGGCGTGGATCGAAACCCGGCGaaa                                      >  2:1406281/1‑101 (MQ=255)
                                                                  cgcgGGGGCGTGGATTGAAACTTTTCGGCATCTGGCAGGCCCGCGCCGCCACGCCGGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAAc                                   <  2:2835161/101‑1 (MQ=37)
                                                                  cgcgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAAc                                   <  2:2727225/101‑1 (MQ=255)
                                                                   gcgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACg                                  >  2:1926622/1‑101 (MQ=255)
                                                                   gcgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACg                                  >  1:405004/1‑101 (MQ=255)
                                                                    cgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACgg                                 >  2:2806614/1‑101 (MQ=255)
                                                                    cgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACgg                                 >  2:3471234/1‑101 (MQ=255)
                                                                    cgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACgg                                 <  1:1020445/101‑1 (MQ=255)
                                                                    cgGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACgg                                 <  1:1642629/101‑1 (MQ=255)
                                                                     gggggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAACACGGc                                >  2:3790591/1‑101 (MQ=255)
                                                                     gggggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGc                                <  1:2320054/101‑1 (MQ=255)
                                                                        ggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACCCGGCGGAAAAAGGCCAt                             >  1:2078549/1‑101 (MQ=255)
                                                                        ggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCAt                             <  1:2130753/101‑1 (MQ=255)
                                                                        ggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCAt                             <  1:779782/101‑1 (MQ=255)
                                                                        ggCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCAt                             >  2:2710934/1‑101 (MQ=255)
                                                                           gTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCCTGcc                          >  1:1587076/1‑101 (MQ=255)
                                                                             ggATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTc                        >  2:3229300/1‑101 (MQ=255)
                                                                              gATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTcc                       >  2:3710461/1‑101 (MQ=255)
                                                                              gATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTcc                       >  1:3960247/1‑101 (MQ=255)
                                                                              gATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTcc                       <  2:3160746/101‑1 (MQ=255)
                                                                              gATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTcc                       <  2:3331581/101‑1 (MQ=255)
                                                                               aTCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCg                      <  1:2639456/101‑1 (MQ=255)
                                                                                tCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCgg                     >  1:3683468/1‑101 (MQ=255)
                                                                                  gAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCgggg                   >  1:4164411/1‑100 (MQ=255)
                                                                                  gAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGgt                   <  2:2255645/101‑1 (MQ=255)
                                                                                      cTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCt               <  1:3969412/101‑1 (MQ=255)
                                                                                      cTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCt               <  2:3439199/101‑1 (MQ=255)
                                                                                      cTTTTCGTCATCTGGCAGGCTCGCGCCGCCACCCTGCGTCGCTCCCCCCGCGGGGGCCTGGATCGGAACCCGGCGAAAAAAGGCCCTGCCTCCGGGCGTcg               >  2:2385449/1‑100 (MQ=37)
                                                                                      cTTTTCATCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCt               <  1:622052/101‑1 (MQ=255)
                                                                                       ttttCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTg              >  1:4116742/1‑101 (MQ=255)
                                                                                       ttttCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGGGTCTg              >  1:179640/1‑101 (MQ=255)
                                                                                        tttCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTgg             <  1:1218836/101‑1 (MQ=255)
                                                                                        tttCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTgg             <  2:1569824/101‑1 (MQ=255)
                                                                                           cGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCAAAAAACGGCCATGCCTCCGGGGGTCTGGCTg          >  2:178666/1‑101 (MQ=255)
                                                                                             tCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTc        >  2:866276/1‑101 (MQ=255)
                                                                                             tCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCAGGGTGTCTGGCTGTc        <  1:3534192/101‑1 (MQ=255)
                                                                                             gcATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTc        <  1:1021087/100‑1 (MQ=255)
                                                                                               aTCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGGGTCTGGCTGTCGc      >  1:3314720/1‑101 (MQ=255)
                                                                                                tCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCt     <  1:380690/101‑1 (MQ=255)
                                                                                                 cTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTc    <  2:1758311/101‑1 (MQ=255)
                                                                                                  tGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTcc   >  2:3929891/1‑101 (MQ=255)
                                                                                                  tGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTcc   <  1:3922258/101‑1 (MQ=255)
                                                                                                   ggCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTccc  <  2:1993809/101‑1 (MQ=255)
                                                                                                                                                                                                        
GCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTCATCTGGCAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCC  >  NC_014034/1396044‑1396243

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG