Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,110,719 G→C 17.2% 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.2% 601.2 / 113.9 227intergenic (+34/+21)cdh/tpiACDP‑diacylglycerol phosphotidylhydrolase/triosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (88/100);  new base C (22/17);  total (110/117)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 231 total aligned reads displayed.

CTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGC  >  NC_000913/4110620‑4110818
                                                                                                   |                                                                                                   
 tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGa                                                                                                    <  2:1192970/100‑1 (MQ=255)
 tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGa                                                                                                    <  2:2532393/100‑1 (MQ=255)
  aCTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGaa                                                                                                   <  1:1556068/100‑1 (MQ=255)
       cGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCg                                                                                              <  1:521780/100‑1 (MQ=255)
         cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGc                                                                                            <  1:613271/100‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGGGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCa                                                                                           <  2:2086108/99‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAc                                                                                          <  2:987351/100‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAc                                                                                          <  1:947496/100‑1 (MQ=255)
           tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAc                                                                                          >  1:649780/1‑100 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAct                                                                                         <  1:2335011/100‑2 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAct                                                                                         <  2:1798734/100‑2 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAct                                                                                         <  1:1410940/100‑2 (MQ=255)
            aaCCGTGCTTCAGCCGAGGAAATTCAAGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTTCCGGATTTCATATCCGGCAct                                                                                         <  2:2931681/100‑2 (MQ=255)
             aCCGTGGTTCAGACGATGAAAGCCAGGTTCATCAGTGTGGGATGTTGCGTTGAGACGGAGAGTTAAGGAAAGTAAGTGACGGTTGTCATATCCGGCAc                                                                                          <  2:2481860/98‑1 (MQ=255)
             aCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActt                                                                                        <  1:1259589/100‑3 (MQ=255)
             aCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCagat                                                                                        >  2:1286010/1‑97 (MQ=255)
              ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                       >  2:591549/1‑100 (MQ=255)
                gTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtact                                                                                     <  2:2894949/100‑6 (MQ=38)
                gTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTc                                                                                     <  2:3172591/100‑1 (MQ=255)
                gTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTc                                                                                     <  1:3797722/100‑1 (MQ=255)
                gTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTc                                                                                     <  1:1938385/100‑1 (MQ=255)
                 cgatctaGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCa                                                                                    <  2:301276/94‑1 (MQ=255)
                  gCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAg                                                                                   <  1:27492/100‑1 (MQ=255)
                  gCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTCAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAg                                                                                   >  1:770036/1‑100 (MQ=255)
                   cTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtactttc                                                                                  <  1:1489651/100‑9 (MQ=37)
                    ttCAGCCGAGGAAGTTCAGGATCATCAGTGTGAGATTTGGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGAc                                                                                 <  1:3710033/100‑1 (MQ=255)
                    ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGAc                                                                                 <  1:3469516/100‑1 (MQ=255)
                    ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGAc                                                                                 <  2:2382798/100‑1 (MQ=255)
                       aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGCAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTa                                                                              <  2:327093/100‑1 (MQ=255)
                       aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCaagatcg                                                                              >  1:301276/1‑94 (MQ=255)
                         ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTaa                                                                             >  2:1555427/1‑99 (MQ=255)
                           gAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGcc                                                                          >  1:267743/1‑100 (MQ=255)
                            aGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCt                                                                         <  1:3848314/100‑1 (MQ=255)
                            aGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCt                                                                         <  1:172172/100‑1 (MQ=255)
                            aGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCt                                                                         >  1:1679355/1‑100 (MQ=255)
                             ggAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATACGGCACCTGTCAGACTTAAGCCTg                                                                        >  2:2696728/1‑100 (MQ=255)
                               aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                      >  1:1720536/1‑100 (MQ=255)
                               aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                      <  2:135371/100‑1 (MQ=255)
                                aaatCAAGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                     <  2:1166217/97‑1 (MQ=255)
                                aaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                     >  1:3439800/1‑100 (MQ=255)
                                  ttCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAg                                                                   >  2:1751166/1‑100 (MQ=255)
                                   tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTCAGTGCCGGATATGAAATCCGGTACCTGTCAGACTTAAGCCTGTTTAGc                                                                  <  2:91130/100‑1 (MQ=255)
                                   tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGc                                                                  <  1:2924222/100‑1 (MQ=255)
                                   tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGc                                                                  <  2:3751503/100‑1 (MQ=255)
                                    cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                 <  1:2528496/100‑1 (MQ=255)
                                    cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGc                                                                  <  2:1614374/99‑1 (MQ=255)
                                       gATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGc                                                               >  2:3109130/1‑99 (MQ=255)
                                       gATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCt                                                              <  2:3216617/100‑1 (MQ=255)
                                           atcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTg                                                          <  2:2004721/100‑1 (MQ=255)
                                           atcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTg                                                          >  1:822712/1‑100 (MQ=255)
                                            tcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTGAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGc                                                         >  1:3745307/1‑100 (MQ=255)
                                               gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGc                                                      <  2:1720527/100‑1 (MQ=255)
                                               gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGc                                                      <  1:940978/100‑1 (MQ=255)
                                                 gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                    >  2:3892624/1‑100 (MQ=255)
                                                 gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                    >  1:1318614/1‑100 (MQ=255)
                                                 gtgAGATTTTGCGTTAAGGCGAAGAGCTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTCAGCCTGTTTAGCCGCTTCTGCAGc                                                      >  2:1086007/1‑98 (MQ=255)
                                                   gagaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTa                                                  >  1:2627327/1‑100 (MQ=255)
                                                         ttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGAtt                                            <  1:91473/100‑1 (MQ=255)
                                                          tGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTa                                           >  2:1946299/1‑100 (MQ=255)
                                                            cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACt                                         <  1:1751002/100‑1 (MQ=255)
                                                              ttAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGc                                       <  2:160267/100‑1 (MQ=255)
                                                              ttAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGACTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGc                                       <  1:2781417/100‑1 (MQ=255)
                                                               tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCg                                      >  2:369230/1‑100 (MQ=255)
                                                               tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCGTCTGCAGCTTTAACGATTACTGCg                                      >  1:1334699/1‑100 (MQ=255)
                                                                  ggCGAAGAGTTAAGGAGAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                   <  1:3109150/100‑1 (MQ=255)
                                                                  ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                   >  1:1614119/1‑100 (MQ=255)
                                                                   gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGAACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAgg                                  >  2:2152780/1‑100 (MQ=255)
                                                                    cGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGc                                 <  1:2557557/100‑1 (MQ=255)
                                                                        gagTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCa                             >  2:1104988/1‑100 (MQ=255)
                                                                          gTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                           >  2:3031844/1‑100 (MQ=255)
                                                                           ttAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCt                          <  2:3495734/100‑1 (MQ=255)
                                                                             aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCttt                        >  1:2427161/1‑100 (MQ=255)
                                                                              aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTc                       <  1:1804976/100‑1 (MQ=255)
                                                                               ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                      >  1:3216884/1‑100 (MQ=255)
                                                                               ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                      <  2:286908/100‑1 (MQ=255)
                                                                               ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                      <  1:3623934/100‑1 (MQ=255)
                                                                                  aagtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCagag                   <  1:2000975/100‑1 (MQ=255)
                                                                                   agtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAgaga                  <  1:107944/100‑1 (MQ=255)
                                                                                   agtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAgaga                  <  1:3301326/100‑1 (MQ=255)
                                                                                     caggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAg                >  2:3587580/4‑100 (MQ=38)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAATGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  2:1251967/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGTCGTCAGCTGTCAGAGAAGcac             >  2:2014441/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  1:414568/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  1:1090104/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  2:427271/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  1:674594/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  1:3664824/1‑100 (MQ=255)
                                                                                        gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  2:2186451/1‑100 (MQ=255)
                                                                                        gtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac             >  1:116461/1‑100 (MQ=255)
                                                                                        gtGCAGGATATGAAATCCGGCCCCAGTTAGACTTAACCTCGTTTAGCCCCTTCGGCAGCTTGAACCAATGCGGCGAAGCCGTCAGCTTTc                       >  2:503981/1‑90 (MQ=255)
                                                                                         tGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc            >  2:624995/1‑100 (MQ=255)
                                                                                         tGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc            <  1:2148766/100‑1 (MQ=255)
                                                                                         tGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc            <  2:400841/100‑1 (MQ=255)
                                                                                          gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca           >  1:3135849/1‑100 (MQ=255)
                                                                                          gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca           >  1:2443938/1‑100 (MQ=255)
                                                                                          gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca           >  2:2798562/1‑100 (MQ=255)
                                                                                           ccGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcaccac          <  1:2207481/100‑1 (MQ=255)
                                                                                            cGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCaccacc         >  1:2225570/1‑100 (MQ=255)
                                                                                            cGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCaccacc         <  1:1691894/100‑1 (MQ=255)
                                                                                              gATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaa       >  1:1058657/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGC  >  NC_000913/4110620‑4110818

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: