Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,110,721 A→T 7.0% intergenic (+36/+19) cdh → / ← tpiA CDP‑diacylglycerol diphosphatase/triose‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,110,7210AT 7.0% 307.7 / 8.2 98intergenic (+36/+19)cdh/tpiACDP‑diacylglycerol diphosphatase/triose‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (44/47);  new base T (4/3);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTT  >  NC_000913/4110578‑4110868
                                                                                                                                               |                                                                                                                                                   
cAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCgg                                                                                                                                               >  1:53231/1‑150 (MQ=255)
cAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCgg                                                                                                                                               <  1:1007257/150‑1 (MQ=255)
 agagTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGc                                                                                                                                              <  1:1232524/150‑1 (MQ=255)
 agagTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGc                                                                                                                                              <  2:824575/150‑1 (MQ=255)
 agagGGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGc                                                                                                                                              <  2:1591008/150‑1 (MQ=255)
  gagTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCa                                                                                                                                             >  1:687709/1‑150 (MQ=255)
    gTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCAct                                                                                                                                           <  1:1632728/150‑2 (MQ=255)
    gTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAcc                                                                                                                                           <  2:572114/150‑1 (MQ=255)
    gTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCAcc                                                                                                                                           <  1:576732/150‑1 (MQ=255)
      gATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTGCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTg                                                                                                                                         <  2:718513/150‑1 (MQ=255)
        tAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTc                                                                                                                                       <  2:1506176/150‑1 (MQ=255)
          acagaTTTGTTCTGCAGGCGACCCAACGCAATCTACTTCCGGTTAACCGGGCATCACCCGACGCAAATCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAg                                                                                                                                     <  2:1040468/146‑1 (MQ=255)
          aCTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATTTCATATCCGGCActtacttt                                                                                                                                     <  2:156285/150‑8 (MQ=255)
          aCTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAg                                                                                                                                     >  2:1567397/1‑150 (MQ=255)
          aCTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAg                                                                                                                                     <  1:1588854/150‑1 (MQ=255)
             cATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACt                                                                                                                                  <  2:1571547/150‑1 (MQ=255)
             cATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACt                                                                                                                                  >  2:538525/1‑150 (MQ=255)
                  gTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGc                                                                                                                             >  2:1082314/1‑150 (MQ=255)
                   tATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGcc                                                                                                                            <  2:53231/150‑1 (MQ=255)
                      tGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGt                                                                                                                         <  1:153390/150‑1 (MQ=255)
                       gCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                                                                        <  2:1058815/150‑1 (MQ=255)
                       gCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGtt                                                                                                                        >  1:144607/1‑150 (MQ=255)
                        cTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                                                                       <  1:1573520/150‑1 (MQ=255)
                        cTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                                                                       <  1:471399/150‑1 (MQ=255)
                        cTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGttt                                                                                                                       <  2:687709/150‑1 (MQ=255)
                          ggCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAg                                                                                                                     >  1:36578/1‑150 (MQ=255)
                             gACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCg                                                                                                                  >  1:1295235/1‑150 (MQ=255)
                               cacaACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCt                                                                                                                >  2:1600294/1‑150 (MQ=255)
                                acaACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGATATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCtt                                                                                                               >  1:267825/1‑150 (MQ=255)
                                   aCGAAATCAACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTACGTGCCGGATATGAAATCCGGCCCCTGTCAGCCTTAAGCCTGTTTAGCCGCTTCTg                                                                                                            >  2:1232470/1‑150 (MQ=255)
                                           tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTa                                                                                                    >  1:306671/1‑150 (MQ=255)
                                             cTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAAc                                                                                                  >  1:291075/1‑150 (MQ=255)
                                              tGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACg                                                                                                 <  1:1368814/150‑1 (MQ=255)
                                               gACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGa                                                                                                >  2:1600239/1‑150 (MQ=255)
                                                   cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTAc                                                                                            <  2:1362959/150‑1 (MQ=255)
                                                    ttAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACt                                                                                           >  1:1143743/1‑150 (MQ=255)
                                                     tAACCGTGCTTCCGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTg                                                                                          >  2:1061935/1‑150 (MQ=255)
                                                     tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTg                                                                                          >  2:1067489/1‑150 (MQ=255)
                                                     tAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTg                                                                                          <  1:1226767/150‑1 (MQ=255)
                                                        ccGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGa                                                                                       >  1:503217/1‑150 (MQ=255)
                                                           tGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAgg                                                                                    >  1:308955/1‑150 (MQ=255)
                                                             cTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCg                                                                                  >  1:1210412/1‑150 (MQ=255)
                                                             cTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCg                                                                                  >  1:129803/1‑150 (MQ=255)
                                                              ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGt                                                                                 >  2:314737/1‑150 (MQ=255)
                                                                 aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAg                                                                              <  2:1295235/150‑1 (MQ=255)
                                                                  gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGc                                                                             >  1:694286/1‑150 (MQ=255)
                                                                   ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAAGctgacgcc                                                                            <  1:1091343/150‑8 (MQ=255)
                                                                        gAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCag                                                                       <  2:528011/150‑1 (MQ=255)
                                                                        gAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCag                                                                       <  2:890859/150‑1 (MQ=255)
                                                                        gAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCag                                                                       <  1:513391/150‑1 (MQ=255)
                                                                         aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCaga                                                                      >  2:1107513/1‑150 (MQ=255)
                                                                         aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCaga                                                                      >  1:323938/1‑150 (MQ=255)
                                                                         aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCaga                                                                      >  1:330001/1‑150 (MQ=255)
                                                                               aGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca                                                                <  1:1082314/150‑1 (MQ=255)
                                                                               aGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca                                                                <  2:291075/150‑1 (MQ=255)
                                                                               aGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGca                                                                <  1:649186/150‑1 (MQ=255)
                                                                                ggATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcac                                                               >  2:298182/1‑150 (MQ=255)
                                                                                   tcatcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcaccac                                                            <  2:350782/150‑1 (MQ=255)
                                                                                      tcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaa                                                         >  2:285387/1‑150 (MQ=255)
                                                                                       caGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaac                                                        >  1:485277/1‑150 (MQ=255)
                                                                                        aGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacc                                                       <  1:1567397/150‑1 (MQ=255)
                                                                                        aGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCaccaacc                                                       <  2:503217/150‑1 (MQ=255)
                                                                                           gtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCCAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAgc                                                    >  1:1118785/1‑150 (MQ=255)
                                                                                            tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAgcg                                                   >  1:1580074/1‑150 (MQ=255)
                                                                                                   ttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACGTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTc                                            <  2:1317733/150‑1 (MQ=255)
                                                                                                          aaGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATcc                                     >  2:47683/1‑150 (MQ=255)
                                                                                                          aaGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATcc                                     <  1:1107513/150‑1 (MQ=255)
                                                                                                             gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGc                                  >  1:891423/1‑150 (MQ=255)
                                                                                                             gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGc                                  <  2:36578/150‑1 (MQ=255)
                                                                                                               gAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTg                                <  2:1210412/150‑1 (MQ=255)
                                                                                                               gAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTg                                <  1:399227/150‑1 (MQ=255)
                                                                                                                  cagTTAAGGAAAGTCAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACGAGCGCGCCGTCGATATCCGGCTGAGc                             <  2:1086497/149‑1 (MQ=255)
                                                                                                                      tAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAAc                         <  1:285387/150‑1 (MQ=255)
                                                                                                                      tAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAAc                         <  2:669446/150‑1 (MQ=255)
                                                                                                                      tAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAAc                         <  1:1410287/150‑1 (MQ=255)
                                                                                                                       aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACa                        >  2:1006014/1‑150 (MQ=255)
                                                                                                                        aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAg                       >  1:8149/1‑150 (MQ=255)
                                                                                                                           aaaGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTc                    >  1:241653/1‑150 (MQ=255)
                                                                                                                             agtaagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTg                  <  2:144607/150‑1 (MQ=255)
                                                                                                                                aggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAg               >  1:651956/3‑150 (MQ=255)
                                                                                                                                aagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAg               >  2:11464/1‑150 (MQ=255)
                                                                                                                                 ggtGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGc              >  1:661734/2‑150 (MQ=255)
                                                                                                                                 agtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGc              <  2:485277/150‑1 (MQ=255)
                                                                                                                                  gtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCGAACAGTTCTGCAGCg             <  1:298182/150‑1 (MQ=255)
                                                                                                                                   tGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGt            >  2:1104561/1‑150 (MQ=255)
                                                                                                                                    gCCGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGtt           >  1:1122397/1‑150 (MQ=255)
                                                                                                                                    gCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGtt           >  2:421669/1‑150 (MQ=255)
                                                                                                                                      cGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAg         <  1:1289230/150‑1 (MQ=255)
                                                                                                                                      cGGATTTCATATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAg         >  2:1426196/1‑150 (MQ=255)
                                                                                                                                       ggATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGa        >  2:561334/1‑150 (MQ=255)
                                                                                                                                           atGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGAcgcg    <  2:594299/150‑1 (MQ=255)
                                                                                                                                           atGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGAcgcg    <  1:314737/150‑1 (MQ=255)
                                                                                                                                           atGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGAcgcg    <  1:1600239/150‑1 (MQ=255)
                                                                                                                                           atGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGAcgcg    >  1:81525/1‑150 (MQ=255)
                                                                                                                                            tGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGt   >  2:982542/1‑150 (MQ=255)
                                                                                                                                             gAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATCTCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGtt  <  1:1600294/150‑1 (MQ=255)
                                                                                                                                             gAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGtt  >  2:1261753/1‑150 (MQ=255)
                                                                                                                                             gAAATCCGGCACCTGTCAGACTTAAGCCTGTATAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGtt  >  1:982243/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTT  >  NC_000913/4110578‑4110868

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