Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,878:1 +G 8.2% 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.2% 343.4 / 29.0 85intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (38/40);  new base G (5/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  NC_000913/2596730‑2597023
                                                                                                                                                     |                                                                                                                                                  
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                     >  1:221637/1‑149 (MQ=255)
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                     <  2:66897/149‑1 (MQ=255)
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTCGCCGTGCCAGCt                                                                                                                                                     >  3:185384/1‑149 (MQ=255)
  aCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑c                                                                                                                                                 >  1:916163/1‑149 (MQ=255)
   cGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑Ca                                                                                                                                                <  3:137094/149‑1 (MQ=255)
      aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACgg                                                                                                                                             >  4:705615/1‑149 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCaa                                                                                                                                          <  3:38102/149‑1 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCaa                                                                                                                                          <  1:227562/149‑1 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTATCTTGCCGTG‑CAGCTGG‑CACGGCaa                                                                                                                                          <  2:106188/149‑1 (MQ=255)
         gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCaa                                                                                                                                          >  3:619124/1‑149 (MQ=255)
         gATACCTCAACAATATATTCATTATACAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCaa                                                                                                                                          <  2:284679/149‑1 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatga                                                                                                                                     >  1:490324/1‑149 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatga                                                                                                                                     >  3:538202/1‑149 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatgat                                                                                                                                    <  2:833417/149‑1 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatgat                                                                                                                                    >  2:857963/1‑149 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatgat                                                                                                                                    >  1:721963/1‑149 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatgat                                                                                                                                    >  3:426490/1‑149 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGAGGTATGAAATGGCGGAGGAATGAGCGTTGAATGCAGATATTCCGAGAGCGCTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgatgat                                                                                                                                    >  4:686092/1‑149 (MQ=255)
                caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGCAATGGCGGAGGAATGAGCGTTCAATGCATATACTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGatgatg                                                                                                                                   <  3:100968/149‑1 (MQ=255)
                         atTCAATATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCACTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGaaaa                                                                                                                          <  4:388969/149‑1 (MQ=255)
                          ttCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGaaaaa                                                                                                                         >  1:901635/1‑149 (MQ=255)
                           tCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAAt                                                                                                                        >  2:21329/1‑149 (MQ=255)
                            cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATc                                                                                                                       >  2:100136/1‑149 (MQ=255)
                             aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      <  2:803790/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  2:176731/149‑1 (MQ=255)
                                ataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   <  1:731760/149‑1 (MQ=255)
                                    aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                               >  1:140776/1‑149 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  1:397968/1‑149 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCaa                                                                                                                                          >  4:967579/1‑121 (MQ=255)
                                     aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCT‑GACACGGCaa                                                                                                                                          <  3:967579/121‑1 (MQ=255)
                                                gtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                                                                                   >  3:127330/1‑149 (MQ=255)
                                                    ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                                                               >  3:873574/1‑149 (MQ=255)
                                                        gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTATCTTGCCGTG‑CAGCTGG‑CACGGCaa                                                                                                                                          >  1:599840/1‑102 (MQ=39)
                                                        gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTATCTTGCCGTG‑CAGCTGG‑CACGGCaa                                                                                                                                          <  2:599840/102‑1 (MQ=39)
                                                          ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑A‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                                                                         >  1:825496/1‑149 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAAc                                                                                                                                >  1:978718/1‑109 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAAc                                                                                                                                <  2:978718/109‑1 (MQ=255)
                                                                 tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                                                                  <  2:715989/149‑1 (MQ=255)
                                                                    cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                               <  3:705615/149‑1 (MQ=255)
                                                                       aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                            <  2:901635/149‑1 (MQ=255)
                                                                         gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTCGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                          <  4:185384/149‑1 (MQ=255)
                                                                          aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATGCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                                                                         <  4:628270/149‑1 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                      <  4:101370/149‑1 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                      <  4:538202/149‑1 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                      <  2:819724/149‑1 (MQ=255)
                                                                                ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                                                                   <  2:443740/149‑1 (MQ=255)
                                                                                 gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTC‑TGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                 >  4:178799/1‑149 (MQ=255)
                                                                                 gATGAAATGGCGGAGGAATGAGCGTGCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                                                                  <  3:688877/149‑1 (MQ=255)
                                                                                       aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTC‑TGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa                                                           <  3:178799/149‑1 (MQ=255)
                                                                                            ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGct                                                       >  4:342968/1‑149 (MQ=255)
                                                                                              aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct                                                     >  4:732671/1‑149 (MQ=255)
                                                                                               ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTg                                                    <  1:42612/149‑1 (MQ=255)
                                                                                               ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTg                                                    <  1:808512/149‑1 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                 >  3:17068/1‑149 (MQ=255)
                                                                                                    gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGa                                               >  2:96570/1‑149 (MQ=255)
                                                                                                      gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                             >  2:909598/1‑149 (MQ=255)
                                                                                                              tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCTACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                     <  4:873574/149‑1 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                    >  4:969841/1‑149 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                    <  2:887722/149‑1 (MQ=255)
                                                                                                                 atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTc                                  >  1:163274/1‑149 (MQ=255)
                                                                                                                  tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcc                                 >  1:643997/1‑149 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCAGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                                >  3:695649/1‑149 (MQ=255)
                                                                                                                     atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAtt                              >  3:127770/1‑149 (MQ=255)
                                                                                                                      ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAttt                             <  1:1054380/149‑1 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                           <  2:397968/149‑1 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                           <  1:535144/149‑1 (MQ=255)
                                                                                                                          gagaGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                         >  3:7082/1‑149 (MQ=255)
                                                                                                                           agagCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTcc                        <  2:140776/149‑1 (MQ=255)
                                                                                                                            gagCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCa                       <  4:7082/149‑1 (MQ=255)
                                                                                                                              gCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                     >  1:116606/1‑149 (MQ=255)
                                                                                                                                  ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCGGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                 >  2:798202/1‑149 (MQ=255)
                                                                                                                                  ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                 >  2:716977/1‑149 (MQ=255)
                                                                                                                                  ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                 <  1:83691/149‑1 (MQ=255)
                                                                                                                                  ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTCCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                 >  2:588370/1‑149 (MQ=255)
                                                                                                                                  tcatCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                 >  4:44454/4‑149 (MQ=255)
                                                                                                                                   tGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTc                <  2:16593/149‑1 (MQ=255)
                                                                                                                                    gTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt               <  2:490324/149‑1 (MQ=255)
                                                                                                                                     tCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCtt              >  1:1001556/1‑149 (MQ=255)
                                                                                                                                       ttGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg            >  3:651723/1‑149 (MQ=255)
                                                                                                                                       ttGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg            >  4:164558/1‑149 (MQ=255)
                                                                                                                                        tGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt           <  1:19046/149‑1 (MQ=255)
                                                                                                                                            gTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCca       >  3:797360/1‑149 (MQ=255)
                                                                                                                                             tGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCcac      <  4:619124/149‑1 (MQ=255)
                                                                                                                                             tGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCcac      <  1:100136/149‑1 (MQ=255)
                                                                                                                                               ccAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAg    >  1:563411/1‑149 (MQ=255)
                                                                                                                                                 aGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAggc  <  2:643997/149‑1 (MQ=255)
                                                                                                                                                  gCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAg                                                          >  3:992325/1‑92 (MQ=255)
                                                                                                                                                  gCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCAGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAg                                                          <  4:992325/92‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  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: ‑