Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,110,721 A→T 17.3% intergenic (+36/+19) 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,7210AT17.3% 550.4 / 103.0 225intergenic (+36/+19)cdh/tpiACDP‑diacylglycerol phosphotidylhydrolase/triosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (89/97);  new base T (22/17);  total (111/114)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 228 total aligned reads displayed.

ACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGC  >  NC_000913/4110622‑4110820
                                                                                                   |                                                                                                   
aCTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGaa                                                                                                     <  1:1556068/100‑1 (MQ=255)
    aCGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATcc                                                                                                 <  1:3472239/100‑1 (MQ=255)
     cGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCg                                                                                                <  1:521780/100‑1 (MQ=255)
     cGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCg                                                                                                <  1:2604644/100‑1 (MQ=255)
      ggttAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCgg                                                                                               >  1:2654597/3‑100 (MQ=255)
      gCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCgg                                                                                               >  1:3262255/1‑100 (MQ=255)
       cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGc                                                                                              <  2:2717949/100‑1 (MQ=255)
       cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGc                                                                                              <  1:2439660/100‑1 (MQ=255)
       cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGc                                                                                              <  1:613271/100‑1 (MQ=255)
         tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAc                                                                                            <  2:987351/100‑1 (MQ=255)
         tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAc                                                                                            <  2:2462078/100‑1 (MQ=255)
         tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAc                                                                                            <  1:2550235/100‑1 (MQ=255)
         tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAc                                                                                            >  1:649780/1‑100 (MQ=255)
         tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAc                                                                                            <  1:947496/100‑1 (MQ=255)
          aaCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAct                                                                                           <  1:2335011/100‑2 (MQ=255)
          aaCCGTGCTTCAGCCGAGGAAATTCAAGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTTCCGGATTTCATATCCGGCAct                                                                                           <  2:2931681/100‑2 (MQ=255)
           aCCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActt                                                                                          <  1:1259589/100‑3 (MQ=255)
            ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAACTCCGGc                                                                                              >  2:3715617/1‑95 (MQ=255)
            ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                         >  2:591549/1‑100 (MQ=255)
             cGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtac                                                                                        <  2:1609641/100‑5 (MQ=38)
              gTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTc                                                                                       <  2:3172591/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)
                    cAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACtt                                                                                 >  1:2414759/1‑100 (MQ=255)
                      gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTaa                                                                               >  1:278169/1‑100 (MQ=255)
                       ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTaa                                                                               >  2:1555427/1‑99 (MQ=255)
                       ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAg                                                                              <  2:1887134/100‑1 (MQ=255)
                         gAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGcc                                                                            <  1:3715645/100‑1 (MQ=255)
                          aGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCt                                                                           >  1:1679355/1‑100 (MQ=255)
                           ggAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCt                                                                           <  2:1217065/99‑1 (MQ=255)
                           ggAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTg                                                                          >  1:2197033/1‑100 (MQ=255)
                           ggAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATACGGCACCTGTCAGACTTAAGCCTg                                                                          >  2:2696728/1‑100 (MQ=255)
                             aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGCAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                        >  1:1720535/1‑100 (MQ=255)
                             aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                        >  1:1720552/1‑100 (MQ=255)
                              aaatCAAGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                       <  2:1166217/97‑1 (MQ=255)
                              aaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                       >  2:2493307/1‑100 (MQ=255)
                              aaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                       <  2:1066413/100‑1 (MQ=255)
                               aTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTa                                                                      <  2:2896723/100‑1 (MQ=255)
                                ttCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAg                                                                     <  1:3390683/100‑1 (MQ=255)
                                 tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTCAGTGCCGGATATGAAATCCGGTACCTGTCAGACTTAAGCCTGTTTAGc                                                                    <  2:91130/100‑1 (MQ=255)
                                 tCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGc                                                                    <  2:3751503/100‑1 (MQ=255)
                                  cAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGcc                                                                   <  1:2528496/100‑1 (MQ=255)
                                   aGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCg                                                                  <  1:3672970/100‑1 (MQ=255)
                                     gATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGc                                                                 >  2:3109130/1‑99 (MQ=255)
                                       tcatcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTc                                                              <  2:2010942/100‑1 (MQ=255)
                                         atcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTg                                                            >  1:822712/1‑100 (MQ=255)
                                          tcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGc                                                           <  2:1951663/100‑1 (MQ=255)
                                           caGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCa                                                          <  2:3031951/100‑1 (MQ=255)
                                             gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAga                                                        <  2:1720526/100‑2 (MQ=255)
                                             gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGc                                                        <  1:940978/100‑1 (MQ=255)
                                             gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGc                                                        <  2:1720543/100‑1 (MQ=255)
                                               gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                      >  2:3892624/1‑100 (MQ=255)
                                               gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCtt                                                      >  1:1318614/1‑100 (MQ=255)
                                                 gagaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTa                                                    >  2:3243669/1‑100 (MQ=255)
                                                 gagaTTTTGCGTTAAGGCGAAGAGTTAAGCAAATTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTa                                                    <  1:2211491/100‑1 (MQ=255)
                                                   gaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAAc                                                  <  1:3761224/100‑1 (MQ=255)
                                                          cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACt                                           <  1:1751002/100‑1 (MQ=255)
                                                             tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCGTCTGCAGCTTTAACGATTACTGCg                                        >  1:1334699/1‑100 (MQ=255)
                                                                ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                     <  2:1754027/100‑1 (MQ=255)
                                                                 gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGAACCTGTCAGACTTAAGCCTCTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAg                                     >  2:2153624/1‑99 (MQ=255)
                                                                   gAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCg                                  <  1:1391647/100‑1 (MQ=255)
                                                                   gAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCg                                  >  1:417140/1‑100 (MQ=255)
                                                                      gagTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCa                               >  2:1104987/1‑100 (MQ=255)
                                                                        gTTAAGGAAAGTAAGTGCCGGATATGAAATGCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                             >  2:2093968/1‑100 (MQ=255)
                                                                        gTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCCGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                             >  1:3845103/1‑100 (MQ=255)
                                                                        gTTAAGGAAAGTAAGTGCCGGATATGAAAGCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                             <  2:412079/100‑1 (MQ=255)
                                                                         ttAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCt                            <  2:3495734/100‑1 (MQ=255)
                                                                           aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCttt                          <  2:1409025/100‑1 (MQ=255)
                                                                            aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTc                         <  1:1804976/100‑1 (MQ=255)
                                                                            aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTc                         <  2:2339314/100‑1 (MQ=255)
                                                                             ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                        >  1:877612/1‑100 (MQ=255)
                                                                             ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                        >  1:3641964/1‑100 (MQ=255)
                                                                             ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                        <  1:2986156/100‑1 (MQ=255)
                                                                             ggAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCa                        <  1:3174201/100‑1 (MQ=255)
                                                                                aagtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCagag                     <  1:2000975/100‑1 (MQ=255)
                                                                                 agtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAgaga                    <  1:107944/100‑1 (MQ=255)
                                                                                   caggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAg                  >  2:3587580/4‑100 (MQ=38)
                                                                                    aagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGc                 <  2:2877883/100‑1 (MQ=255)
                                                                                     ggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca                >  1:1914416/2‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAATGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  2:1251967/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  1:674594/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  1:3892110/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  1:1090104/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  2:427271/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac               >  1:1578046/1‑100 (MQ=255)
                                                                                      gtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCCGCAGCTTTAACGATTACTGCCAAGGCGTCAGCTTTCAGAGAAGcac               >  1:3089175/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)
                                                                                       tGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc              >  2:3379351/1‑100 (MQ=255)
                                                                                       tGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacc              <  1:2148766/100‑1 (MQ=255)
                                                                                        gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcacca             >  2:2040312/1‑100 (MQ=255)
                                                                                          cGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCaccacc           <  1:1691894/100‑1 (MQ=255)
                                                                                            gATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaa         >  1:1058657/1‑100 (MQ=255)
                                                                                              tatGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacc       <  2:1292354/100‑1 (MQ=255)
                                                                                              tatGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacc       >  1:3124827/1‑100 (MQ=255)
                                                                                              tatGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacc       <  1:3457067/100‑1 (MQ=255)
                                                                                                tGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAg     >  1:2761684/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGC  >  NC_000913/4110622‑4110820

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