Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,488,059 T→G 5.5% intergenic (+24/‑173) ydcF → / → aldA DUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,488,0590TG 5.5% 325.8 / 6.6 108intergenic (+24/‑173)ydcF/aldADUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base T (50/51);  new base G (2/4);  total (52/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 108 total aligned reads displayed.

ACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAAC  >  NC_000913/1487911‑1488205
                                                                                                                                                    |                                                                                                                                                  
aCGGTCCGCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCTTGGCAAACGCTTAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAATATTGCCCGTTTGTGAACCACTTg                                                                                                                                                   <  2:122640/150‑1 (MQ=255)
aCGGTCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGt                                                                                                                                                  <  1:811138/150‑1 (MQ=255)
         gcgGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAAc                                                                                                                                          <  2:1784109/150‑1 (MQ=255)
         gcgGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAAc                                                                                                                                          <  2:80057/150‑1 (MQ=255)
           gggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGCTGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACgg                                                                                                                                        <  2:284280/150‑1 (MQ=255)
            ggCGAGATTGTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACggg                                                                                                                                       <  2:870310/150‑1 (MQ=255)
             gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                      <  2:686787/150‑1 (MQ=255)
             gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                      <  1:1261128/150‑1 (MQ=255)
              cGAGATTTTATCGTTCACGTTGATTATCCGGCAGAAGTCATCCATGCATGGCAAACGCTGACACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                     <  2:223201/150‑1 (MQ=255)
               gagaTTTTATCGTGCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGATACGTTAAAAATTGCCCGTATGCAAACAAGTGGTTCACAAACGGGCaa                                                                                                                                    <  1:286036/150‑2 (MQ=255)
                         cGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAc                                                                                                                          >  1:558732/1‑150 (MQ=255)
                          gTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACt                                                                                                                         <  2:41192/150‑1 (MQ=255)
                          gTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACt                                                                                                                         <  2:502637/150‑1 (MQ=255)
                            tCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACttt                                                                                                                       <  2:836889/150‑1 (MQ=255)
                            tCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACAGGATGCGGGGCGCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACttt                                                                                                                       <  1:862840/150‑1 (MQ=255)
                               cGTTGATTTTCCGGCAGAAGTCATCCCTGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                                    >  1:299415/1‑150 (MQ=255)
                                gTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAtt                                                                                                                   >  2:629571/1‑150 (MQ=255)
                                  tGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTc                                                                                                                 <  2:1487650/150‑1 (MQ=255)
                                    aTTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAGTTTTATTTCTg                                                                                                               >  1:1744145/1‑150 (MQ=255)
                                    aTTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTg                                                                                                               <  1:1594952/150‑1 (MQ=255)
                                         ccGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCttt                                                                                                          >  1:183673/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAATCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACAACTTGTTTGCAAACGGGCATGACTCATGACTTTTATTTCTGCCtatt                                                                                                         >  1:120817/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                         >  1:917733/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                         >  1:998387/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                         >  2:1773708/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                         >  1:505289/1‑150 (MQ=255)
                                          cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTGCTGCCgttt                                                                                                         >  2:1273952/1‑150 (MQ=255)
                                           ggCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTa                                                                                                        <  1:1060481/150‑1 (MQ=255)
                                            gCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAt                                                                                                       <  2:399157/150‑1 (MQ=255)
                                               gAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTcc                                                                                                    >  1:1490941/1‑150 (MQ=255)
                                                   tCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                                <  1:1750745/150‑1 (MQ=255)
                                                   tCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGGCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                                <  1:348544/150‑1 (MQ=255)
                                                   gcATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                                <  1:1316025/149‑1 (MQ=255)
                                                      tCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGACCCACTTGTTTGCACACGGGCATGCCTGCTGACTTTTATTTCTGCCTTTTATTCCTTTTaca                                                                                             >  2:487310/1‑150 (MQ=255)
                                                        catgcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACt                                                                                           >  1:461240/1‑150 (MQ=255)
                                                           gcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTTGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGt                                                                                        >  2:880068/1‑150 (MQ=255)
                                                             atgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGttt                                                                                      >  2:834923/1‑150 (MQ=255)
                                                                  aaaCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGTCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATg                                                                                 <  2:1021508/150‑1 (MQ=255)
                                                                   aaCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGa                                                                                <  2:1406257/150‑1 (MQ=255)
                                                                     cGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAg                                                                              >  1:43308/1‑150 (MQ=255)
                                                                      gCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGc                                                                             >  1:230131/1‑150 (MQ=255)
                                                                       cTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGcc                                                                            >  1:1740283/1‑150 (MQ=255)
                                                                           aaCATGATTCGGTGCTCATAGAGGCGATGGAAATTCGCTCGTTACGTTAAAAAATGCACGTATGTGAACCACTTGTTAGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTc                                                                        <  1:976086/150‑1 (MQ=255)
                                                                             cATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTcac                                                                      >  1:1350741/1‑150 (MQ=255)
                                                                               tgatgCGGTGCTCCTCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                                    <  2:325709/150‑1 (MQ=255)
                                                                               tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                                    <  2:1490941/150‑1 (MQ=255)
                                                                               tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                                    <  2:786007/150‑1 (MQ=255)
                                                                               tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCGTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                                    >  2:905772/1‑150 (MQ=255)
                                                                                  tgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAAt                                                                 <  1:829907/150‑1 (MQ=255)
                                                                                   gCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAAtt                                                                >  1:185832/1‑150 (MQ=255)
                                                                                    cGGTGCTCATCGAGGCGATGGAAAGTCTCTAGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTg                                                               <  1:345542/150‑1 (MQ=255)
                                                                                      cttcTCATCGAGGCGAGATAAAGTAGTTCGTTACGGTAAAAATTGCGCGTTTGTTAACAACTAGTAAGCAAACGTTCATGACTCCTTTCTTTTATTTCTGCCTATTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTc                                                             <  1:228309/147‑1 (MQ=255)
                                                                                       tGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTcc                                                            >  1:908172/1‑150 (MQ=255)
                                                                                       tGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTcc                                                            >  1:902365/1‑150 (MQ=255)
                                                                                           cATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTc                                                        >  1:386841/1‑150 (MQ=255)
                                                                                           cATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTc                                                        <  2:1249287/150‑1 (MQ=255)
                                                                                               gAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGCCTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGa                                                    >  2:365024/1‑150 (MQ=255)
                                                                                                aGGCGATGGAAAGTCGCTCGTTACGATAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAACACTTCACAGAATTGTCCTTTCACGAt                                                   >  2:908009/1‑150 (MQ=255)
                                                                                                    gATGGAAAGTCGTTCGTTACGTTAAAAATTGCCCGTTTGTGAAACACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACTCTTGATTTTATGAATCCCTTCACACAATTGTCCTTTCACGATTCCg                                               >  2:1413223/1‑150 (MQ=255)
                                                                                                    gATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                               <  1:87481/150‑1 (MQ=255)
                                                                                                    gATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                               <  2:505289/150‑1 (MQ=255)
                                                                                                    gATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                               >  2:1297673/1‑150 (MQ=255)
                                                                                                    gATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                               >  1:753206/1‑150 (MQ=255)
                                                                                                      tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                             <  2:507274/150‑1 (MQ=255)
                                                                                                      tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                             >  2:624736/1‑150 (MQ=255)
                                                                                                      tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                             >  2:443437/1‑150 (MQ=255)
                                                                                                      tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                             >  2:1669968/1‑150 (MQ=255)
                                                                                                       ggAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtct                                            >  1:825960/1‑150 (MQ=255)
                                                                                                        gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAAGCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                           <  1:321627/150‑1 (MQ=255)
                                                                                                         aaaGTCGCTCGTTAAGTTAAAAATTGCCCGTTTGTGCACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAATTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGgctct                                          <  2:9197/150‑1 (MQ=255)
                                                                                                          aaGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctctc                                         <  2:1760615/150‑1 (MQ=255)
                                                                                                              cGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgat                                     >  1:795688/1‑150 (MQ=255)
                                                                                                               gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                                    >  1:564949/1‑150 (MQ=255)
                                                                                                                cTCGTTACGTTAAAATTTGCCCGTTTGTGAAACACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatga                                   <  1:1218507/150‑1 (MQ=255)
                                                                                                                   gttacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattg                                <  1:465604/150‑1 (MQ=255)
                                                                                                                     tacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattgat                              >  2:1368771/1‑150 (MQ=255)
                                                                                                                     tacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattgat                              <  2:230131/150‑1 (MQ=255)
                                                                                                                      acgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGTCATGACTCCTGACTTTTATTTCTGCATTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                             <  1:1154666/150‑1 (MQ=255)
                                                                                                                      acgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                             <  1:17228/150‑1 (MQ=255)
                                                                                                                       cgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTGCTGCCTTTTATTCATTTTAAACTTGTTTTTATGAAGCCCTTCACAGAATAGTCATTTCACGATTCCGTCTATCTTATGATTGATGt                            >  2:340946/1‑150 (MQ=255)
                                                                                                                        gttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                           <  2:1326798/150‑1 (MQ=255)
                                                                                                                        gttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCATTTCACGATTCCGTCTCTCTGATGATTGATGtt                           >  2:1225455/1‑150 (MQ=255)
                                                                                                                         ttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtta                          <  2:1392694/150‑1 (MQ=255)
                                                                                                                           aaaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaat                        <  1:834923/150‑1 (MQ=255)
                                                                                                                           aaaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTAACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaat                        <  2:677/150‑1 (MQ=255)
                                                                                                                           aaaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTATATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCATTTCACTATTCAGTCTCTCTGATGATTGATGttaat                        >  1:1550842/1‑150 (MQ=255)
                                                                                                                              tcatGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaattaa                     >  1:743034/4‑150 (MQ=38)
                                                                                                                              aaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaattaa                     >  1:1238996/1‑150 (MQ=255)
                                                                                                                               aTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc                    <  1:1773708/150‑1 (MQ=255)
                                                                                                                                ttGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACa                   >  1:1762/1‑150 (MQ=255)
                                                                                                                                ttGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACa                   >  1:1439061/1‑150 (MQ=255)
                                                                                                                                atGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACa                   >  1:1318895/2‑150 (MQ=39)
                                                                                                                                   cccGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATg                <  1:1359212/150‑1 (MQ=255)
                                                                                                                                     cGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTa              <  2:1022168/150‑1 (MQ=255)
                                                                                                                                      gTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTAt             <  2:1305639/150‑1 (MQ=255)
                                                                                                                                         tgtgAACCACTTGTTTGCAAACGGGCATTACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCTTCTCTCTGATGATTGATGTTAATTAACAATGTATTCa          <  1:487310/150‑1 (MQ=255)
                                                                                                                                           tgAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCAcc        <  2:1087385/150‑1 (MQ=255)
                                                                                                                                               ccACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACATTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGaaa    >  1:6489/1‑150 (MQ=255)
                                                                                                                                               ccACTTGTTTGCAAACGGGAATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGaaa    <  1:163963/150‑1 (MQ=255)
                                                                                                                                                 aCTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAaaac  <  2:917733/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAAC  >  NC_000913/1487911‑1488205

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