Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,878:1 +G 8.0% intergenic (‑141/+27) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8781.G 8.0% 512.9 / 40.3 125intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (53/62);  new base G (4/5);  total (57/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 130 total aligned reads displayed.

CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  NC_000913/2596730‑2597023
                                                                                                                                                     |                                                                                                                                                 
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                    <  1:843132/149‑1 (MQ=255)
 aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                                                                  >  1:434730/1‑149 (MQ=255)
 aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                                                                  >  4:149496/1‑149 (MQ=255)
    gCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                                                               >  2:871589/1‑149 (MQ=255)
    gCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGACAGCTGCAc                                                                                                                                               >  4:850320/1‑149 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATGTTCCGAGAGCGGTTGTCTTCCCGTGCAGCTAGCACGGCaa                                                                                                                                          <  2:790017/149‑1 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                                                                          <  2:103414/149‑1 (MQ=255)
           tACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATACATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAga                                                                                                                                        <  2:413742/149‑1 (MQ=255)
            aCCTCAACAATATGTTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                                                                       <  4:685691/149‑1 (MQ=255)
             ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatg                                                                                                                                      >  2:903901/1‑149 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                     >  3:199798/1‑149 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                     >  3:21098/1‑149 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGCCTTGCCGTGCAGCTGGCACGGCaa                                                                                                                                          <  3:1009771/144‑1 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGGGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                                                                          >  4:1009771/1‑144 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatgat                                                                                                                                    <  4:543325/149‑1 (MQ=255)
                caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                                                                                   <  1:1064410/149‑1 (MQ=255)
                caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTACACGGCAAGatgatg                                                                                                                                   <  4:281582/149‑1 (MQ=255)
                   caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAAc                                                                                                                                >  4:201609/1‑149 (MQ=255)
                     atatatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAg                                                                                                                              >  2:959983/1‑149 (MQ=255)
                         atTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGATCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCGGCTGCACGGCAAGATTATGAACAGaaaa                                                                                                                          >  2:843396/1‑149 (MQ=255)
                          ttCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaaa                                                                                                                         >  1:64443/1‑149 (MQ=255)
                           tCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAt                                                                                                                        <  1:959983/149‑1 (MQ=255)
                           tCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAt                                                                                                                        <  2:809877/149‑1 (MQ=255)
                             aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      <  3:872030/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     >  3:488241/1‑149 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  1:2713/149‑1 (MQ=255)
                                 taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTcc                                                                                                                  >  2:730170/1‑149 (MQ=255)
                                   gAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                                                                >  4:890032/1‑149 (MQ=255)
                                   gAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                                                                        <  4:660496/125‑1 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  1:7683/149‑1 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  2:310530/1‑149 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  1:80882/149‑1 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  4:453552/1‑149 (MQ=255)
                                      aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                                                                             >  3:381886/1‑149 (MQ=255)
                                       tGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGTGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     >  3:951492/1‑140 (MQ=255)
                                       tGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGTGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  4:951492/140‑1 (MQ=255)
                                        gATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                                                                       >  1:297808/1‑121 (MQ=255)
                                        gATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                                                                       <  2:297808/121‑1 (MQ=255)
                                          tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                                                                         >  1:31708/1‑149 (MQ=255)
                                            gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                                                                       >  1:207666/1‑149 (MQ=255)
                                               ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAg                                                                                                    <  4:313880/149‑1 (MQ=255)
                                                        gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                                                                           <  3:890032/149‑1 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                        <  4:576673/149‑1 (MQ=255)
                                                            aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                                                                       <  1:1015269/149‑1 (MQ=255)
                                                             aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                                                                      <  2:443264/149‑1 (MQ=255)
                                                              gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGACAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                                                                     >  2:849790/1‑149 (MQ=255)
                                                              gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTCGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                                                                     <  4:906424/149‑1 (MQ=255)
                                                               aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                                                                    >  4:520952/1‑149 (MQ=255)
                                                                ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                                   >  4:422681/1‑149 (MQ=255)
                                                                  aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatg                                                                                                                                      <  2:548377/96‑1 (MQ=255)
                                                                  aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatg                                                                                                                                      >  1:548377/1‑96 (MQ=255)
                                                                   ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                                                                <  4:403163/149‑1 (MQ=255)
                                                                   ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                                                                <  4:488241/149‑1 (MQ=255)
                                                                        ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   <  4:675430/109‑1 (MQ=255)
                                                                        ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                                                                           >  3:560288/1‑149 (MQ=255)
                                                                        ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCATGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   >  3:675430/1‑109 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     >  3:434740/1‑102 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  4:434740/102‑1 (MQ=255)
                                                                                      aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc                                                             >  4:441308/1‑149 (MQ=255)
                                                                                      aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc                                                             >  3:161388/1‑149 (MQ=255)
                                                                                         ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  2:229302/1‑97 (MQ=255)
                                                                                         ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  1:229302/97‑1 (MQ=255)
                                                                                           cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc                                                        >  1:305004/1‑149 (MQ=255)
                                                                                            ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  3:987999/1‑94 (MQ=255)
                                                                                                gaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAGAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg                                                   >  4:965621/1‑149 (MQ=255)
                                                                                                   tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg                                                <  3:496691/149‑1 (MQ=255)
                                                                                                      gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                             >  3:438749/1‑149 (MQ=255)
                                                                                                         ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                          >  1:459788/1‑149 (MQ=255)
                                                                                                         ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                          <  3:850320/149‑1 (MQ=255)
                                                                                                         ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                          <  3:125992/149‑1 (MQ=255)
                                                                                                           cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                        >  1:364125/1‑149 (MQ=255)
                                                                                                            aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                            <  3:494161/112‑1 (MQ=255)
                                                                                                            aaTGCATAGATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATATGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCt                                       <  4:860457/149‑1 (MQ=255)
                                                                                                              tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACTGAGGCTCTGGCGGAAACGGACttt                                     <  1:937750/149‑1 (MQ=255)
                                                                                                              tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                     <  3:441308/149‑1 (MQ=255)
                                                                                                              tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                     <  2:459788/149‑1 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                    >  1:938848/1‑149 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                    <  1:374004/149‑1 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                    <  1:903901/149‑1 (MQ=255)
                                                                                                                  tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                 <  4:742513/117‑1 (MQ=255)
                                                                                                                  tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                 >  3:742513/1‑117 (MQ=255)
                                                                                                                  tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcc                                 <  2:1072833/149‑1 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                                >  4:348359/1‑149 (MQ=255)
                                                                                                                          gagaGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                         >  2:413803/1‑149 (MQ=255)
                                                                                                                              gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                     <  1:730170/149‑1 (MQ=255)
                                                                                                                                    gTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt               <  3:552223/149‑1 (MQ=255)
                                                                                                                                     tCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCtt              <  3:629259/149‑1 (MQ=255)
                                                                                                                                      cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCttt             <  3:520952/149‑1 (MQ=255)
                                                                                                                                       ttGCCGTGCAGCTGGCACGGCAAGATGACGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                        <  4:358401/121‑1 (MQ=255)
                                                                                                                                       ttGCCGTGCAGCTGGCACGGCAAGATGACGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                        >  3:358401/1‑121 (MQ=255)
                                                                                                                                        tGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt           <  2:64443/149‑1 (MQ=255)
                                                                                                                                        tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt           >  4:377809/1‑149 (MQ=255)
                                                                                                                                         gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc          <  1:104781/149‑1 (MQ=255)
                                                                                                                                         gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc          >  4:21219/1‑149 (MQ=255)
                                                                                                                                          ccGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc         <  3:786513/149‑1 (MQ=255)
                                                                                                                                          ccGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGCCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc         <  3:178930/149‑1 (MQ=255)
                                                                                                                                            gTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCca       >  1:1076618/1‑149 (MQ=255)
                                                                                                                                                cAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAgg   >  3:231350/1‑149 (MQ=255)
                                                                                                                                                 aGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAggc  >  2:215979/1‑149 (MQ=255)
                                                                                                                                                 aGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAggc  >  2:616203/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  NC_000913/2596730‑2597023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: