Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,488,057 C→G 8.9% intergenic (+22/‑175) 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,0570CG 8.9% 244.2 / 10.7 83intergenic (+22/‑175)ydcF/aldADUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base C (33/42);  new base G (5/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAAC  >  NC_000913/1487909‑1488205
                                                                                                                                                    |                                                                                                                                                    
cTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACt                                                                                                                                                     >  1:644037/1‑150 (MQ=255)
    ggTGCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTGATCAATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                                 <  2:946516/150‑1 (MQ=255)
    ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                                 <  2:852607/150‑1 (MQ=255)
    ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                                 >  1:761271/1‑150 (MQ=255)
      tCCCCGCGGGCGAGATTTTATCGTTCACGTGGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTg                                                                                                                                               <  2:318750/150‑1 (MQ=255)
       ccccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGc                                                                                                                                              <  2:557078/150‑1 (MQ=255)
        cccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCa                                                                                                                                             <  1:538579/150‑1 (MQ=255)
        cccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCa                                                                                                                                             >  1:1060540/1‑150 (MQ=255)
          cgcgGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCaaa                                                                                                                                           <  2:1069543/150‑1 (MQ=255)
          cgcgGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCaaa                                                                                                                                           >  1:1008342/1‑150 (MQ=255)
             gggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACgg                                                                                                                                        <  1:680616/150‑1 (MQ=255)
              ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATTCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACggg                                                                                                                                       <  1:407998/150‑1 (MQ=255)
              ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACggg                                                                                                                                       <  1:276583/150‑1 (MQ=255)
              ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACggg                                                                                                                                       <  2:283987/150‑1 (MQ=255)
              ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCCGTGCTCATCGAGGCGATGGAATGTCGCTCGTGACGTTAAAAATTGCCCGTTTGTGAACCACTTGTGTGCAAACggg                                                                                                                                       >  2:192305/1‑150 (MQ=255)
                cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                     <  2:18716/150‑1 (MQ=255)
                        tATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCt                                                                                                                             >  1:553103/1‑150 (MQ=255)
                        tATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAAAGGGCATGACTCCt                                                                                                                             <  1:1021948/150‑1 (MQ=255)
                             ttCACGTTGATTTTACGGCAGAAGTCATCCATGCATGGCAAACGATGCAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACtt                                                                                                                        >  1:966620/1‑150 (MQ=255)
                                aCGTTGAGTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAACCATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTa                                                                                                                     <  1:76418/150‑1 (MQ=255)
                                 cGTTGATTTTCCGGCAGAAGTCATCCCTGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                                    <  1:362445/150‑1 (MQ=255)
                                 cGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                                    >  2:966267/1‑150 (MQ=255)
                                    ggATGTTCCGGCAGATGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTAAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTc                                                                                                                 <  2:820791/149‑1 (MQ=255)
                                     gATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCt                                                                                                                >  1:772769/1‑150 (MQ=255)
                                       ctttCCGGCCGAAGTCATCCATGCATGGCAAACGCTGAAACATGCTGCGGTGCTCATCGAGGCGATGGACAGTCGCTCGTTACGTTAAAAATTGCCCGCTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                              <  1:187424/149‑1 (MQ=255)
                                        tttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCGCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGcc                                                                                                             <  1:192165/150‑1 (MQ=255)
                                                 gAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTcc                                                                                                    >  2:940251/1‑150 (MQ=255)
                                                 gAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTcc                                                                                                    >  1:303254/1‑150 (MQ=255)
                                                  aaGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCt                                                                                                   <  1:306927/150‑1 (MQ=255)
                                                      cATCCATGCATGGCAAACGCTGAAACATGATGCGGTGATCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTACTGACTTTTATTTCTGCCTTTTATTCCTTTTa                                                                                               >  1:424869/1‑150 (MQ=255)
                                                      cATCCATGATTGGCATACGGTTTAACTTGATCCGGTGCTCATCGAGGGGATGGAATTGCGCTCGTTACGTTAAAAATTGCCGGTTTGTGTACCACTTGTTTGCAAACGGTCATGACTCCTGACTTTTATTTCTGCCTTTTATTGCTTTTa                                                                                               <  2:426016/150‑1 (MQ=255)
                                                        tCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTaca                                                                                             >  2:1250095/1‑150 (MQ=255)
                                                            tgcatgGCAATCGCTGAATCATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTg                                                                                         <  1:25134/150‑1 (MQ=255)
                                                            tgcatgGCAAACGCTGAAACATGAAGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTg                                                                                         <  2:633120/150‑1 (MQ=255)
                                                             gcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGt                                                                                        >  1:292970/1‑150 (MQ=255)
                                                              catgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGtt                                                                                       >  2:1059607/1‑150 (MQ=255)
                                                                tgGCAATCGCTGAAACATGATGCGGTGATCATCGAGGCGATGGAATGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTTCAACCGGGCATGCCTCCTGCAGTGTAGTTCGGCCTTTTATTCCTTTTACACTTGtttt                                                                                     >  2:207463/1‑150 (MQ=255)
                                                                  gCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTa                                                                                   >  2:736340/1‑150 (MQ=255)
                                                                  gCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTa                                                                                   >  1:492133/1‑150 (MQ=255)
                                                                  gCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTa                                                                                   <  2:510204/150‑1 (MQ=255)
                                                                  gCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTa                                                                                   >  2:378376/1‑150 (MQ=255)
                                                                  gCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTa                                                                                   <  2:1142379/150‑1 (MQ=255)
                                                                       cGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCATTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAg                                                                              <  2:1181979/150‑1 (MQ=255)
                                                                       cGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAg                                                                              <  2:364178/150‑1 (MQ=255)
                                                                       cGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCGGCCTTTTATTTCTGCCTTTTTTTCCTTTTACACTTGttttatgtagc                                                                              >  1:267073/1‑143 (MQ=255)
                                                                        gCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGc                                                                             >  1:501012/1‑150 (MQ=255)
                                                                         cTGAAACATGCAGCGCGGCTCATCGAGGCGATGGAAAGGCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGcc                                                                            <  2:858482/150‑1 (MQ=255)
                                                                         cTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCACGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTTTTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGcc                                                                            >  2:429204/1‑150 (MQ=255)
                                                                              aCATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTca                                                                       <  1:941940/150‑1 (MQ=255)
                                                                                atgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTAAGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTcaca                                                                     <  1:38502/150‑1 (MQ=255)
                                                                                 tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTAAACTTGTTTTTATGAAGCCCTTCACAg                                                                    >  1:199404/1‑150 (MQ=255)
                                                                                               tCGAGGCGATGGAAAGGCGCTCGTTACGTTAAAATTTGCCCGTTTGTGAACCAGTTGTTTGCAAACGGGCATGACTCCTAACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGTATTGTCCTTTCAc                                                      <  2:470208/150‑1 (MQ=255)
                                                                                                       aTGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCCTGACTCCTGACTTTTAGTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGACGCCCTTCACAGAATTGTCCTTTCACGATTCAGt                                              >  2:244102/1‑150 (MQ=255)
                                                                                                        tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                             <  1:1246380/150‑1 (MQ=255)
                                                                                                          gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                           >  1:416641/1‑150 (MQ=255)
                                                                                                          gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                           >  1:1208750/1‑150 (MQ=255)
                                                                                                            aaGTCGCTCGTTACGTTAAAAATTTCCCGTTTTTTAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctctc                                         <  1:1247604/150‑1 (MQ=255)
                                                                                                                 gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                                    <  2:417668/150‑1 (MQ=255)
                                                                                                                 gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                                    <  1:46878/150‑1 (MQ=255)
                                                                                                                  gtCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACTGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCC‑‑ACACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgatt                                 <  2:127031/149‑1 (MQ=255)
                                                                                                                  cTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatga                                   >  2:872486/1‑150 (MQ=255)
                                                                                                                    cgttacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgatt                                 <  1:1059607/150‑1 (MQ=255)
                                                                                                                      ttacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattga                               >  1:695167/1‑150 (MQ=255)
                                                                                                                        acgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                             <  2:1008342/150‑1 (MQ=255)
                                                                                                                         cgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCATACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGt                            >  1:541630/1‑150 (MQ=255)
                                                                                                                         cgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGt                            <  1:960074/150‑1 (MQ=255)
                                                                                                                           ttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCATTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCAGGATTCCGTCTATCTGATGATTGATGtta                          >  2:76687/1‑150 (MQ=255)
                                                                                                                           ttaAAAATTGCCCGTTTCTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAATCCTTATACACTTGTTTTAATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtta                          <  2:1016962/150‑1 (MQ=255)
                                                                                                                            taAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaa                         >  1:428695/1‑150 (MQ=255)
                                                                                                                               aaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatta                      <  1:571933/150‑1 (MQ=255)
                                                                                                                                aaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaattaa                     <  2:1060540/150‑1 (MQ=255)
                                                                                                                                aaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaattaa                     >  2:331512/1‑150 (MQ=255)
                                                                                                                                  ttGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACa                   <  1:955048/150‑1 (MQ=255)
                                                                                                                                   tGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTATACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACaa                  <  2:553103/150‑1 (MQ=255)
                                                                                                                                   tGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACaa                  >  1:543628/1‑150 (MQ=255)
                                                                                                                                   tGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGGCTCTCTGATGATTGATGTTAATTAACaa                  >  1:982645/1‑150 (MQ=255)
                                                                                                                                   tGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTAATTCTGCCTTTTAATCCTTTTACACTTGTTTTTATGAAGCCCTTCTCAGAATTGTCCTTTCACGATTCCGTGTCTCTGATGATTGATGTTTATTAACaa                  >  1:827909/1‑150 (MQ=255)
                                                                                                                                    gCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAAt                 >  2:505432/1‑150 (MQ=255)
                                                                                                                                      ccGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGt               >  1:713820/1‑150 (MQ=255)
                                                                                                                                         tttGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTTCCTTTTATTCCTTTTACACTAGTTTTTATTAAGCCCTTCACAGAATTGTGCTTTCACGATTCCCTCTCTCTGATGATTGATGTTAATTAACAATGTAtt            <  1:672707/150‑1 (MQ=255)
                                                                                                                                             tgAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCAcc        <  1:1078915/150‑1 (MQ=255)
                                                                                                                                               aaCCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGa      >  2:46009/1‑150 (MQ=255)
                                                                                                                                                   aCTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAaaac  <  1:900964/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAAC  >  NC_000913/1487909‑1488205

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