Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,110,719 G→C 17.1% intergenic (+34/+21) cdh → / ← tpiA CDP‑diacylglycerol phosphotidylhydrolase/triosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,110,7190GC17.1% 454.2 / 85.5 176intergenic (+34/+21)cdh/tpiACDP‑diacylglycerol phosphotidylhydrolase/triosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (80/65);  new base C (14/16);  total (94/82)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 178 total aligned reads displayed.

CTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAG  >  NC_000913/4110620‑4110817
                                                                                                   |                                                                                                  
cTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATg                                                                                                    <  1:339499/100‑1 (MQ=255)
 tACTGACGCTTAACCGTGCTTCAGACGTGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGa                                                                                                   >  2:658029/1‑100 (MQ=255)
    tGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGATAAGGAAAGTAAGTGCCGGATATGAAAt                                                                                                >  2:2795484/1‑100 (MQ=255)
     gACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATc                                                                                               >  2:2731437/1‑100 (MQ=255)
     gACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATc                                                                                               >  1:1012130/1‑100 (MQ=255)
       cGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCg                                                                                             >  1:430863/1‑100 (MQ=255)
       cGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCg                                                                                             <  2:1870276/100‑1 (MQ=255)
         cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGc                                                                                           <  1:3245356/100‑1 (MQ=255)
          ttAACCGTGCTTCAGCCGAGGAAATGCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCa                                                                                          <  2:3375218/100‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAc                                                                                         <  2:473087/100‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAc                                                                                         <  2:1879139/100‑1 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCATGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAcc                                                                                        <  2:2888770/100‑1 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAcc                                                                                        <  1:1451133/100‑1 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAATTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAcc                                                                                        <  2:1948408/100‑1 (MQ=255)
              ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                      >  1:3001094/1‑100 (MQ=255)
              ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                      <  1:204930/100‑1 (MQ=255)
              ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                      >  1:3508123/1‑100 (MQ=255)
               cGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtac                                                                                     <  1:579983/100‑5 (MQ=38)
                   cTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtactttc                                                                                 <  1:457935/100‑9 (MQ=37)
                    ttCAGCCGAGGAAGTTCAGGATCATCAGTGTGAGATTTAGCGTTAAGGCGAAGAGTTAAGGAAAGTAAATGCAGGATATGAAATCCGGAACCTGTCAg                                                                                  >  2:2245932/1‑98 (MQ=255)
                    ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActcactttcc                                                                                <  1:1911390/100‑10 (MQ=37)
                    ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGAc                                                                                <  1:3480428/100‑1 (MQ=255)
                    ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGAc                                                                                <  1:1564320/100‑1 (MQ=255)
                     tCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACt                                                                               >  2:1403000/1‑100 (MQ=255)
                      cAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACtt                                                                              >  2:3350352/1‑100 (MQ=255)
                       aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTa                                                                             >  1:810388/1‑100 (MQ=255)
                        gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTaa                                                                            >  1:1869978/1‑100 (MQ=255)
                        gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTaa                                                                            <  2:252999/100‑1 (MQ=255)
                         ccGGGGAAATTCAGGATCATCAGTGTGAGAGTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTACGTGCCGGCTTTCATATGCGGCActgactttcc                                                                                <  2:2059604/95‑10 (MQ=255)
                          cGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGc                                                                          <  2:431879/100‑1 (MQ=255)
                           gAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGcc                                                                         <  1:1465047/100‑1 (MQ=255)
                                  ttCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTa                                                                   >  1:764847/1‑99 (MQ=255)
                                  ttCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAg                                                                  <  2:3491242/100‑1 (MQ=255)
                                   tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGc                                                                 >  2:2740495/1‑100 (MQ=255)
                                    cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                <  2:693933/100‑1 (MQ=255)
                                    cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                <  2:5025/100‑1 (MQ=255)
                                    cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                <  2:24146/100‑1 (MQ=255)
                                     aGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                <  1:1950575/99‑1 (MQ=255)
                                      ggATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGc                                                              >  1:2900924/1‑100 (MQ=255)
                                         tcatcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTc                                                           >  2:568987/1‑100 (MQ=255)
                                           atcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTg                                                         <  1:2404728/100‑1 (MQ=255)
                                           atcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTg                                                         >  2:2531983/1‑100 (MQ=255)
                                               gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCa                                                       <  1:1646115/98‑1 (MQ=255)
                                                 gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                   <  1:2423374/100‑1 (MQ=255)
                                                 gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                   >  1:728179/1‑100 (MQ=255)
                                                   gagaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGACATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCTGCTTCTGCAGCTTTa                                                 >  2:365748/1‑100 (MQ=255)
                                                   gagaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTa                                                 >  1:1747627/1‑100 (MQ=255)
                                                     gaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAAc                                               <  2:1747619/100‑1 (MQ=255)
                                                      aTTTTGCGTTTAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACg                                              <  1:3020182/100‑1 (MQ=255)
                                                      aTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACg                                              <  2:3019552/100‑1 (MQ=255)
                                                      aTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACg                                              >  1:2149519/1‑100 (MQ=255)
                                                       ttttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGa                                             <  2:3242666/100‑1 (MQ=255)
                                                       ntttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACg                                              >  2:1679468/2‑99 (MQ=255)
                                                         ttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGAtt                                           <  2:819487/100‑1 (MQ=255)
                                                           gCGTTAAGGCGAAGAGCTAAGGAAAGTAAGTGCCGGATATGAAATGCGGCACCTGTCAGAGTTAAGCCGGTTTCTACGCTGCTGCAGCTTTAACGAtt                                           >  2:209392/1‑98 (MQ=255)
                                                            cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACt                                        >  2:526599/1‑100 (MQ=255)
                                                             gTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTg                                       <  2:3390225/100‑1 (MQ=255)
                                                               tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCg                                     <  2:3508098/100‑1 (MQ=255)
                                                               tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCg                                     <  1:931803/100‑1 (MQ=255)
                                                               tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGAAATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCg                                     >  2:2189475/1‑100 (MQ=255)
                                                                  ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                  >  1:2181761/1‑100 (MQ=255)
                                                                  ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                  >  2:626847/1‑100 (MQ=255)
                                                                  ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                  >  1:2748829/1‑100 (MQ=255)
                                                                   gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAgg                                 <  1:3242842/100‑1 (MQ=255)
                                                                   gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAgg                                 <  2:2752903/100‑1 (MQ=255)
                                                                    cGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGc                                >  1:393044/1‑100 (MQ=255)
                                                                         agTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCg                               >  2:264794/1‑96 (MQ=255)
                                                                          gTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                          >  1:2861083/1‑100 (MQ=255)
                                                                             aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCttt                       <  1:208332/100‑1 (MQ=255)
                                                                             aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCttt                       >  1:3048875/1‑100 (MQ=255)
                                                                               ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                     >  1:1619154/1‑100 (MQ=255)
                                                                                 aaaGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCaga                   >  1:111181/1‑100 (MQ=255)
                                                                                  aagtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCagan                  <  1:1679476/100‑2 (MQ=255)
                                                                                   agtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAgaga                 >  1:2503440/1‑100 (MQ=255)
                                                                                    gtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCagag                  >  1:788214/1‑98 (MQ=255)
                                                                                    gtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGaa                >  1:1411757/1‑100 (MQ=255)
                                                                                     taagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAg               >  1:1649126/1‑100 (MQ=255)
                                                                                     caggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAg               >  1:2230610/4‑100 (MQ=38)
                                                                                      aggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGc              >  1:3540111/3‑100 (MQ=255)
                                                                                      aagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGc              <  2:592622/100‑1 (MQ=255)
                                                                                       ggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca             >  2:76439/2‑100 (MQ=255)
                                                                                       ggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca             >  1:649538/2‑100 (MQ=255)
                                                                                       ggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGATTTAACGATTACTGCGAAGGCGTCAGCTTTATgagactc              >  1:504167/2‑96 (MQ=255)
                                                                                       agtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca             <  2:1292343/100‑1 (MQ=255)
                                                                                       agtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca             <  2:355114/100‑1 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCCGCTTTCAGAGAAGc              >  2:413725/1‑98 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac            >  1:1039575/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac            >  1:1351834/1‑100 (MQ=255)
                                                                                        gtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac            >  1:2413493/1‑100 (MQ=255)
                                                                                         tGCCGGATTTCATATTCGGCACCTGTCAGACTTACCCCTGTTTAGCCGCTACTGCAGCTTGAACGATCACTACGAAGGCGTCAGCTTTTAGAGAAGca             >  2:1500063/1‑98 (MQ=255)
                                                                                         tGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc           >  1:2687197/1‑100 (MQ=255)
                                                                                          gCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGATTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca          >  2:1221711/1‑100 (MQ=255)
                                                                                          gCCGGATATGAAATCCTGCACCAGACAGACTTAAGCCTGTTTAGCAGCGTCTGCAGCTGTAACGATTACTGCGCAGGCGCCAGCTATCAGAGAAGcac            >  2:2919884/1‑98 (MQ=255)
                                                                                          gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacta          <  1:2472719/100‑3 (MQ=255)
                                                                                          gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca          >  1:2394724/1‑100 (MQ=255)
                                                                                                tatGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacn    <  2:2153170/100‑2 (MQ=255)
                                                                                                tatGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCATCAGCTTTCAGAGAAGCACCaccaacc    >  1:3200044/1‑100 (MQ=255)
                                                                                                  tGAAATCCGGCACCTGTCAGCCTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTGTCAGAGAAGCACCACCAACCAg  >  2:2572234/1‑100 (MQ=255)
                                                                                                  tGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAg  >  1:2983637/1‑100 (MQ=255)
                                                                                                  tGAAATCCGGCACCCGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAg  >  2:66059/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAG  >  NC_000913/4110620‑4110817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: