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.

TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977
                                                                                                    |                                                                                                    
tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                       <  1:2562783/100‑1 (MQ=255)
 ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                     <  1:1612470/100‑1 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                                   >  2:1789888/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGCGGAATGAGCGTTCAATGCATATATTACAAAAGAGGTTGTCTTGCAGTACCAGCTGCa                                                                                                  >  1:1567600/1‑99 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                 >  1:1035666/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGAGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                 >  2:420374/1‑100 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATTAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                               <  1:1763535/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                               <  2:3025142/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCCGAGGAATGAGCGCTCACTGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCTGCActg                                                                                               <  2:2733397/100‑3 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                              >  1:914912/1‑100 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATAATCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGc                                                                                              <  2:2088874/100‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                             <  1:3659331/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCa                                                                                             <  1:1702887/99‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCaa                                                                                            <  2:3362360/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                           <  2:1063221/100‑1 (MQ=39)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                           <  1:1789344/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATTAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                            <  1:157510/98‑1 (MQ=38)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                          <  2:112122/100‑1 (MQ=39)
            aGATTACCACAGTAAAGCGGGATGAAATGGCGGAGAAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTTTTGTCGTGCAGCTGTCACGGCAAga                                                                                          <  2:1340360/99‑1 (MQ=255)
            aGATTACCACAGGAAAGCGGGATGAAATGGCTGAGGAATGAGTGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgac                                                                                         <  1:425757/100‑2 (MQ=38)
            aGATTACCACAGGAAAGCGGGATGAAATGGCGGATGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                         <  2:3584958/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGTAATGTGTGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                                     <  1:2025990/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGTTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                                     <  2:2922871/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                                     <  1:2862420/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgagatcgg                                                                                    >  2:1625400/1‑94 (MQ=255)
                  ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                   >  1:2979088/1‑100 (MQ=255)
                  ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                   >  1:1328141/1‑100 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCGTATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGTTGCACGGCAAGatgatgn                                                                                    <  2:1472361/98‑2 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAAc                                                                                  >  2:2092686/1‑100 (MQ=255)
                    acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACa                                                                                 >  1:2805414/1‑100 (MQ=255)
                       ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaa                                                                              >  2:1960488/1‑100 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:827157/98‑1 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGAACGGCAAGATGATGAACAGAAAAATCAg                                                                       >  1:1896557/1‑100 (MQ=255)
                                  tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                   <  1:3928108/100‑1 (MQ=255)
                                      aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:1680002/100‑1 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:643120/99‑1 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                              >  2:3716152/1‑100 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:819551/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            <  1:2089497/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            >  1:3345634/1‑100 (MQ=255)
                                          cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGAACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  1:2115685/1‑100 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  2:1002369/100‑1 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       >  2:3878792/1‑100 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       <  1:1585025/100‑1 (MQ=255)
                                                aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  2:2181783/1‑94 (MQ=255)
                                                 aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    >  1:2523336/1‑100 (MQ=255)
                                                 aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                    <  1:2618391/100‑1 (MQ=255)
                                                  tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                   >  2:594369/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                                  >  2:1324838/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                                  >  2:3426859/1‑100 (MQ=255)
                                                   gAGCGTTCAATGCATATATTCCGACAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACACAggcgg                                                  >  2:1750816/1‑100 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  1:3141658/1‑100 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  2:3733112/1‑100 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGgcggcg                                                <  1:2347987/100‑1 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCATGCCAGCTGCACGGCAAAATGATGAACAGAAAAATAAAAACAGATATACACACAggcag                                                  >  1:1527994/1‑96 (MQ=255)
                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  1:1464725/100‑1 (MQ=255)
                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  1:2159204/100‑1 (MQ=255)
                                                         tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             >  2:2050963/1‑99 (MQ=255)
                                                          cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                           >  2:323034/1‑100 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  1:17663/1‑100 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        <  1:1642536/100‑1 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  1:1088092/1‑100 (MQ=255)
                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        >  1:2254469/1‑100 (MQ=255)
                                                              gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:549585/100‑1 (MQ=255)
                                                                 tatatTCCGAGTGCGGTTGTCTTGCCGTTGCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACTCAGTCGGCGGGCGACGGCtt                                     <  2:3109334/99‑1 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  1:368478/100‑1 (MQ=255)
                                                                   tatTCCTAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                  <  2:3789894/100‑1 (MQ=255)
                                                                   tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                  <  1:2956207/100‑1 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  1:2037387/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:1879520/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 <  1:618008/100‑1 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:2174530/1‑100 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGGCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 <  1:1168530/100‑1 (MQ=255)
                                                                      tCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               <  2:1481503/100‑1 (MQ=255)
                                                                       ccGAGAGCGGTTGTCCTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                              >  1:3629211/1‑100 (MQ=255)
                                                                          agagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           >  1:1210002/1‑100 (MQ=255)
                                                                           gagCGGTTGTCTTGCCGTGCCAGCTGCACAGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACTCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                          <  2:2346134/100‑1 (MQ=255)
                                                                            agCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATACGCTTCGa                         <  1:2052105/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        >  1:802164/1‑100 (MQ=255)
                                                                              tcatcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  1:3247983/7‑100 (MQ=37)
                                                                              cGGTTGTCTTGCCGTTCCAGCTACAAGGCAAAATGTTGAACAGAAAAATCAGTCCAGCTGTAAAACCAGGCGGCGCGAGACCGCTTCATAGGCTTCg                          >  2:2970719/1‑97 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  1:1324852/100‑1 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  1:2305809/100‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:3018103/100‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  1:2110543/100‑1 (MQ=255)
                                                                                 tcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATAAGTCCAGCTGCACACCCAGGCGGCGGCCAACGGCTACATAGGCTTCGAACa                      >  2:2164891/4‑98 (MQ=21)
                                                                                  tGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   <  1:1422833/100‑1 (MQ=255)
                                                                                  catCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   >  1:1371416/3‑100 (MQ=38)
                                                                                    tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:1840532/1‑100 (MQ=39)
                                                                                     cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                   <  2:2409970/97‑1 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  2:3409649/1‑99 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:3768944/1‑100 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  1:3412079/1‑100 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGACAAGATGATGAACACAAAACTCAGTCCAGCTGTACACCCAGGCGGCGCGCGACGGCGTCATAGGCTTCCATCATGcccc                >  2:313550/1‑96 (MQ=255)
                                                                                       tGCCGTGCAGCTGGCACGACAAGATAATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCATGCCAcc               >  2:951252/1‑99 (MQ=255)
                                                                                                aGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg     >  2:2271079/1‑100 (MQ=255)
                                                                                                  cTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:1608732/100‑1 (MQ=255)
                                                                                                  cTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:21518/100‑1 (MQ=255)
                                                                                                   tGCACGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  2:2700030/100‑1 (MQ=255)
                                                                                                   tGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  2:3716019/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  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: