Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 15.0% intergenic (‑142/+26) 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.G15.0% 409.0 / 76.2 120intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (44/58);  new base G (8/10);  total (52/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 = 8.77e-01

Only 100 of 128 total aligned reads displayed.

GGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596780‑2596977
                                                                                                   |                                                                                                   
ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                    >  1:1369555/1‑100 (MQ=255)
 gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                                   <  2:671450/100‑1 (MQ=255)
  ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCa                                                                                                  >  1:4012353/1‑100 (MQ=255)
   tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                                 <  2:4029296/100‑1 (MQ=255)
    cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                <  1:3372287/100‑1 (MQ=255)
     agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATACCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                      >  2:4008988/1‑94 (MQ=255)
     agagGAAGAGTCCCACAGGAAAGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                <  2:312324/99‑1 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              >  1:1598223/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              >  2:4111737/1‑100 (MQ=255)
      gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              >  1:2450010/1‑100 (MQ=255)
       agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCTAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                             <  1:2910523/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                            <  1:3957111/100‑1 (MQ=255)
        ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                            <  1:3954041/100‑1 (MQ=255)
         gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                           <  2:2818642/100‑1 (MQ=39)
          aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                          <  2:3629197/100‑1 (MQ=39)
          aaGACTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgn                                                                                          <  2:3435717/100‑2 (MQ=38)
             aTTACCACAGTAAAGCGGGATGAAATGGCGGAGTAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgacaa                                                                                       <  2:3079042/100‑4 (MQ=37)
             aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                       <  2:1422764/100‑1 (MQ=255)
               tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                                     <  1:1442667/100‑1 (MQ=255)
                aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                    <  2:900973/100‑1 (MQ=255)
                 ccACAGGAAAGCGGGATGAAATGTCGGGGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGCCTTTCCATGCCAGCT‑GCACGGCAAGCTGATGaa                                                                                   >  2:839972/1‑100 (MQ=255)
                 ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                   >  2:2980374/1‑100 (MQ=255)
                 ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                   <  1:1541919/100‑1 (MQ=255)
                 ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATGCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                   <  2:3286752/100‑1 (MQ=255)
                   acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                                 <  2:2319215/100‑1 (MQ=255)
                       gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                                 <  2:3442211/96‑1 (MQ=255)
                       gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaa                                                                             <  2:982567/100‑1 (MQ=255)
                         aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                           <  2:2321515/100‑1 (MQ=255)
                          aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAt                                                                          <  1:1670894/100‑1 (MQ=255)
                             gtgttGAAACGGCGCAAGAAAGTGCGCCCAATACATATATTCCAACAGCAGAAAACTTTCATAGCCACCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:3842819/94‑1 (MQ=255)
                             gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:808416/98‑1 (MQ=255)
                              ggATGAAATGGCGGAGGAATGAGCGTTCATTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACCGCAAGATGTTGAACAGAAAAATCAg                                                                       >  2:2967746/1‑99 (MQ=255)
                              ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                                      <  2:1782950/100‑1 (MQ=255)
                              ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                                      >  2:609992/1‑100 (MQ=255)
                               gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                     >  2:2123658/1‑100 (MQ=255)
                               gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                     >  2:3845826/1‑100 (MQ=255)
                                    aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                >  1:2105899/1‑100 (MQ=255)
                                    aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                <  1:1781411/100‑1 (MQ=255)
                                     aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  2:1369546/100‑1 (MQ=255)
                                      tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                              >  2:2424265/1‑100 (MQ=255)
                                       ggCGGGGGAATGAGCGTTCAATGCATATATGCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:3546406/100‑1 (MQ=255)
                                        gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            >  1:3991109/1‑100 (MQ=255)
                                             ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       >  2:1131190/1‑100 (MQ=255)
                                               aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                                     <  1:2767857/100‑1 (MQ=255)
                                                  gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                                  >  1:4104667/1‑100 (MQ=255)
                                                   aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAACAATCAGTCCAGCTGTACACCCAggcggc                                                 >  2:405365/1‑100 (MQ=255)
                                                   aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAggcggc                                                 >  1:3645608/1‑100 (MQ=255)
                                                     cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                               <  2:3497952/100‑1 (MQ=255)
                                                       ttCAATGCATATATTCCGAGAGCGGGTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  2:2119529/99‑1 (MQ=255)
                                                          aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          <  1:4009014/100‑1 (MQ=255)
                                                            tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  2:1491075/1‑98 (MQ=255)
                                                            tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  1:4101242/1‑100 (MQ=255)
                                                             gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACATAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:2259981/100‑1 (MQ=255)
                                                             gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:905289/100‑1 (MQ=255)
                                                              cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                      <  1:3852959/100‑1 (MQ=255)
                                                              cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                      >  1:593488/1‑100 (MQ=255)
                                                               atatatTCCGAGAGCGGTTGTCTTGCCGTGCCCGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCAGCGGGCGACGGCtt                                     >  1:1191729/1‑100 (MQ=255)
                                                               atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     <  1:3086132/100‑1 (MQ=255)
                                                               atatacTCCGAGAGCGGTTGTCTTGCAGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     <  2:1934799/100‑1 (MQ=255)
                                                                tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACATCCAGGCGGCGGGCGACGGCTTc                                    <  2:2397303/100‑1 (MQ=255)
                                                                  tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                  <  1:4045802/100‑1 (MQ=255)
                                                                   atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  1:1480350/1‑100 (MQ=255)
                                                                     tCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               <  2:2298218/100‑1 (MQ=255)
                                                                      ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                              >  1:185462/1‑100 (MQ=255)
                                                                      ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                              >  1:353403/1‑100 (MQ=255)
                                                                      ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                              <  1:45794/100‑1 (MQ=255)
                                                                       cGAGAGCGGTTGTCTTGCCGTGCCCGAT‑GCACGGCAAGATGATGAACAGAAAACTCAGTCCCGCTGTACACCCACGCGGCGGGCGACCGCTTCATa                                 >  2:1937782/1‑96 (MQ=255)
                                                                       cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             <  1:3125359/100‑1 (MQ=255)
                                                                        gagaGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                            >  2:2983361/1‑100 (MQ=255)
                                                                         agagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           <  1:3126502/100‑1 (MQ=255)
                                                                          tgttcatcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                          >  1:1625306/10‑100 (MQ=25)
                                                                              ggTTGTCTTGCCGTGCCAGCG‑GCACGGCAAGATGATGAACAGAATAATCAGTCCAGCTGTACACCCCGGCGGCGGGCGACGGCTGCATAGGCTTCGagc                       <  2:1586762/99‑3 (MQ=255)
                                                                               atcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGTACAGAAAAATCAGTCCCGCTGTACACCCAGGCGGCGGGCGACGGCTTCATCGGCTTCGATCa                      >  2:1896401/5‑99 (MQ=25)
                                                                                ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCatg                    <  1:2790746/100‑3 (MQ=255)
                                                                                ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:328248/100‑1 (MQ=255)
                                                                                ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  1:1481684/100‑1 (MQ=255)
                                                                                 tGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   >  1:2258994/1‑100 (MQ=255)
                                                                                 tGTCTTGCCGTGCCAGCT‑CCACGGCAAGAGGATGAACAGAAAGATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  2:679220/95‑1 (MQ=255)
                                                                                  gTCTTGCAGTGCCAGCT‑GCACGTCTAGATGATGAACAGAAAAAGCACTCCAGCTGTTCACCCAGACTGCGGGCGACGGCTTTTTAGGCTTCGATCAggtc                  <  2:4057725/100‑3 (MQ=255)
                                                                                  atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGCCGCCCCCATCGGCGTCGATCACGcc                  >  1:3736486/2‑100 (MQ=25)
                                                                                  atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  >  1:4065558/2‑100 (MQ=39)
                                                                                      tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAGCg              >  1:1929419/1‑100 (MQ=255)
                                                                                      tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:1307454/1‑100 (MQ=255)
                                                                                      tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:2969488/1‑100 (MQ=255)
                                                                                      tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCGGTACACCCAGGCGGCGGGCGACTGTTTCATAGACTTCGATCAGGCCACCg              <  2:1331228/100‑1 (MQ=255)
                                                                                      tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCAAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:63719/1‑100 (MQ=255)
                                                                                        ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTTATAGGCTTCGATCAGGCCACCGAg            <  2:565697/100‑1 (MQ=255)
                                                                                         cGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg           <  1:3674851/100‑1 (MQ=255)
                                                                                          gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  1:3419854/100‑1 (MQ=255)
                                                                                             ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct       >  1:1260836/1‑100 (MQ=255)
                                                                                               aGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCATCCTAGGCTTCGATCAGGCCACCGAGGCGCTgg     >  2:3689139/1‑100 (MQ=255)
                                                                                                 cT‑GGACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:3137885/100‑1 (MQ=255)
                                                                                                 cT‑GCACGGCAAGATGATGAACAGAACAATCAGTCCAGCTGTACACCCAGGCAGCGGGCGACGGCTACATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   >  2:2294045/1‑100 (MQ=255)
                                                                                                 cT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:3811643/100‑1 (MQ=255)
                                                                                                 cT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   >  1:3412622/1‑100 (MQ=255)
                                                                                                  t‑GCGCGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCGGGCGGCGGGCGACGGCTTCAGAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:1731904/99‑1 (MQ=255)
                                                                                                  t‑GCACGGCAAGATGCTGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   >  2:2020268/1‑99 (MQ=255)
                                                                                                  t‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  >  1:1790143/1‑100 (MQ=255)
                                                                                                  t‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  >  2:2408008/1‑100 (MQ=255)
                                                                                                  t‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  >  1:508399/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596780‑2596977

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