Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 15.4% 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.4% 498.6 / 94.4 149intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (62/64);  new base G (12/11);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

Only 100 of 153 total aligned reads displayed.

TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977
                                                                                                    |                                                                                                    
tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                       <  1:2562783/100‑1 (MQ=255)
 ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                     <  1:1612470/100‑1 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑c                                                                                                   >  2:1789888/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGCGGAATGAGCGTTCAATGCATATATTACAAAAGAGGTTGTCTTGCAGTACCAGCT‑G‑Ca                                                                                                  >  1:1567600/1‑99 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CAc                                                                                                 >  1:1035666/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGAGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CAc                                                                                                 >  2:420374/1‑100 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATTAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACgg                                                                                               <  1:1763535/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACgg                                                                                               <  2:3025142/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCCGAGGAATGAGCGCTCACTGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCT‑G‑CActg                                                                                               <  2:2733397/100‑3 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGc                                                                                              >  1:914912/1‑100 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATAATCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGc                                                                                              <  2:2088874/100‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCa                                                                                             <  1:3659331/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCa                                                                                             <  1:1702887/99‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCaa                                                                                            <  2:3362360/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCAAg                                                                                           <  2:1063221/100‑1 (MQ=39)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCAAg                                                                                           <  1:1789344/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATTAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCaa                                                                                            <  1:157510/98‑1 (MQ=38)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCAAga                                                                                          <  2:112122/100‑1 (MQ=39)
            aGATTACCACAGTAAAGCGGGATGAAATGGCGGAGAAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTTTTGTCGTG‑CAGCT‑GTCACGGCAAga                                                                                          <  2:1340360/99‑1 (MQ=255)
            aGATTACCACAGGAAAGCGGGATGAAATGGCTGAGGAATGAGTGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGG‑CACGGCAAgac                                                                                         <  1:425757/100‑2 (MQ=38)
            aGATTACCACAGGAAAGCGGGATGAAATGGCGGATGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgat                                                                                         <  2:3584958/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGTAATGTGTGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGatgatg                                                                                     <  1:2025990/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGTTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGatgatg                                                                                     <  2:2922871/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGatgatg                                                                                     <  1:2862420/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAgagatcgg                                                                                    >  2:1625400/1‑94 (MQ=255)
                  ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGaa                                                                                   >  1:2979088/1‑100 (MQ=255)
                  ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGaa                                                                                   >  1:1328141/1‑100 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCGTATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGTT‑G‑CACGGCAAGatgatgn                                                                                    <  2:1472361/98‑2 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAAc                                                                                  >  2:2092686/1‑100 (MQ=255)
                    acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACa                                                                                 >  1:2805414/1‑100 (MQ=255)
                       ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGaa                                                                              >  2:1960488/1‑100 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:827157/98‑1 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑AACGGCAAGATGATGAACAGAAAAATCAg                                                                       >  1:1896557/1‑100 (MQ=255)
                                  tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                   <  1:3928108/100‑1 (MQ=255)
                                      aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:1680002/100‑1 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:643120/99‑1 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                              >  2:3716152/1‑100 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:819551/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            <  1:2089497/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            >  1:3345634/1‑100 (MQ=255)
                                          cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑AACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  1:2115685/1‑100 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  2:1002369/100‑1 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       >  2:3878792/1‑100 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       <  1:1585025/100‑1 (MQ=255)
                                                aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  2:2181783/1‑94 (MQ=255)
                                                 aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    >  1:2523336/1‑100 (MQ=255)
                                                 aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    <  1:2618391/100‑1 (MQ=255)
                                                  tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                   >  2:594369/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                                  >  2:1324838/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                                  >  2:3426859/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGACAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACACAggcgg                                                  >  2:1750816/1‑100 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  1:3141658/1‑100 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  2:3733112/1‑100 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGgcggcg                                                <  1:2347987/100‑1 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCATGCCAGCT‑G‑CACGGCAAAATGATGAACAGAAAAATAAAAACAGATATACACACAggcag                                                  >  1:1527994/1‑96 (MQ=255)
                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  1:1464725/100‑1 (MQ=255)
                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  1:2159204/100‑1 (MQ=255)
                                                         tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             >  2:2050963/1‑99 (MQ=255)
                                                          cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                           >  2:323034/1‑100 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  1:17663/1‑100 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        <  1:1642536/100‑1 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  1:1088092/1‑100 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  1:2254469/1‑100 (MQ=255)
                                                              gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:549585/100‑1 (MQ=255)
                                                                 tatatTCCGAGTGCGGTTGTCTTGCCGTTGCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACTCAGTCGGCGGGCGACGGCtt                                     <  2:3109334/99‑1 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  1:368478/100‑1 (MQ=255)
                                                                   tatTCCTAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                  <  2:3789894/100‑1 (MQ=255)
                                                                   tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                  <  1:2956207/100‑1 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  1:2037387/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:1879520/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 <  1:618008/100‑1 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:2174530/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGGCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 <  1:1168530/100‑1 (MQ=255)
                                                                      tCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               <  2:1481503/100‑1 (MQ=255)
                                                                       ccGAGAGCGGTTGTCCTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                              >  1:3629211/1‑100 (MQ=255)
                                                                          agagCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           >  1:1210002/1‑100 (MQ=255)
                                                                           gagCGGTTGTCTTGCCGTGCCAGCT‑G‑CACAGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACTCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                          <  2:2346134/100‑1 (MQ=255)
                                                                            agCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATACGCTTCGa                         <  1:2052105/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        >  1:802164/1‑100 (MQ=255)
                                                                              tcatcatCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  1:3247983/7‑100 (MQ=37)
                                                                              cGGTTGTCTTGCCGTTCCAGCT‑A‑CAAGGCAAAATGTTGAACAGAAAAATCAGTCCAGCTGTAAAACCAGGCGGCGCGAGACCGCTTCATAGGCTTCg                          >  2:2970719/1‑97 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  1:1324852/100‑1 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  1:2305809/100‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:3018103/100‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  1:2110543/100‑1 (MQ=255)
                                                                                 tcatCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATAAGTCCAGCTGCACACCCAGGCGGCGGCCAACGGCTACATAGGCTTCGAACa                      >  2:2164891/4‑98 (MQ=21)
                                                                                  tGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   <  1:1422833/100‑1 (MQ=255)
                                                                                  catCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   >  1:1371416/3‑100 (MQ=38)
                                                                                    tCTTGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:1840532/1‑100 (MQ=39)
                                                                                     cTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   <  2:2409970/97‑1 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  2:3409649/1‑99 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:3768944/1‑100 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGG‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:3412079/1‑100 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGG‑CACGACAAGATGATGAACACAAAACTCAGTCCAGCTGTACACCCAGGCGGCGCGCGACGGCGTCATAGGCTTCCATCATGcccc                >  2:313550/1‑96 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGG‑CACGACAAGATAATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCATGCCAcc               >  2:951252/1‑99 (MQ=255)
                                                                                                aGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg     >  2:2271079/1‑100 (MQ=255)
                                                                                                  cT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:1608732/100‑1 (MQ=255)
                                                                                                  cT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:21518/100‑1 (MQ=255)
                                                                                                   t‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  2:2700030/100‑1 (MQ=255)
                                                                                                   t‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  2:3716019/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑G‑CACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977

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