Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,874 Δ1 bp 16.7% intergenic (‑137/+31) 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,8730C.16.7% 416.3 / 86.4 150intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (61/63);  new base . (13/12);  total (74/76)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-02

Only 100 of 152 total aligned reads displayed.

TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGC  >  NC_000913/2596775‑2596969
                                                                                                  |                                                                                                 
 aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                 >  1:223889/1‑100 (MQ=255)
  ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGTCGTGCCAg                                                                                                <  1:1902058/100‑1 (MQ=255)
  ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                >  2:1983575/1‑100 (MQ=255)
     ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                            <  1:1612470/100‑1 (MQ=255)
      gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGc                                                                                           >  2:1789888/1‑100 (MQ=255)
       ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCa                                                                                          >  1:2593353/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)
          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)
            agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGTAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                     <  1:907107/100‑1 (MQ=255)
            agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                     <  1:3659331/100‑1 (MQ=255)
            agGAAGATTACCACAGAAAAGCGGGATGAAATGGCGGAGGAATAAGCGTACAATGCATATATTGCGAGCGCGGTTGTCTTGCCGTGCCAGCTGCaaagca                                                                                     >  2:2546260/1‑95 (MQ=255)
             ggAAGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCa                                                                                     <  1:1702887/99‑1 (MQ=255)
             ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCaa                                                                                    <  2:3362360/100‑1 (MQ=255)
              gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                   <  2:1063221/100‑1 (MQ=39)
              gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAg                                                                                   <  2:430576/100‑1 (MQ=39)
               aaGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATTAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                    <  1:157510/98‑1 (MQ=38)
                aGATTACCACAGGAAAGCGGGATGAAATGGCTGAGGAATGAGTGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAgac                                                                                 <  1:425757/100‑2 (MQ=38)
                aGATTACCACAGGAAAGCGGGATGAAATGGCGGATGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgat                                                                                 <  2:3584958/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGTAATGTGTGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                             <  1:2025990/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGatgatg                                                                             <  1:2862420/100‑1 (MQ=255)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgagatcgg                                                                            >  2:1625400/1‑94 (MQ=255)
                      ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                           >  1:1328141/1‑100 (MQ=255)
                      ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                           >  1:2979088/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)
                            gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                     >  1:1784522/1‑100 (MQ=255)
                                  gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                 <  2:827157/98‑1 (MQ=255)
                                  gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGAACGGCAAGATGATGAACAGAAAAATCAg                                                               >  1:1896557/1‑100 (MQ=255)
                                     aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGATGCACGGCAAGATGATGAACAGAAAAATCAGTcc                                                            >  2:2052087/1‑100 (MQ=255)
                                      tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                           <  1:3928108/100‑1 (MQ=255)
                                      tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                           >  2:3724762/1‑100 (MQ=255)
                                        aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                         <  2:2157249/100‑1 (MQ=255)
                                          aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                       <  1:1680002/100‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCTAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                      <  1:2256755/100‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                       <  1:643120/99‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                      >  2:3716152/1‑100 (MQ=255)
                                             gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                    >  1:3345634/1‑100 (MQ=255)
                                              cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGAACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                   >  1:2115685/1‑100 (MQ=255)
                                               ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTacac                                                  <  2:2992326/100‑1 (MQ=255)
                                                  ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                               <  1:1585025/100‑1 (MQ=255)
                                                     aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                            <  1:2618391/100‑1 (MQ=255)
                                                     aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                            >  1:2523336/1‑100 (MQ=255)
                                                      tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                           >  2:594369/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                          >  2:3426859/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                          >  2:1324838/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGACAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACACAggcgg                                          >  2:1750816/1‑100 (MQ=255)
                                                        aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                         >  1:3141658/1‑100 (MQ=255)
                                                         gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGgcggcg                                        <  2:3806133/100‑1 (MQ=255)
                                                         gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCATGCCAGCTGCACGGCAAAATGATGAACAGAAAAATAAAAACAGATATACACACAggcag                                          >  1:1527994/1‑96 (MQ=255)
                                                           gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                      <  1:2159204/100‑1 (MQ=255)
                                                           gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                      <  1:1464725/100‑1 (MQ=255)
                                                             tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTAAACGGCAAGATGATGAACAGAAAAAACAGTCCAGCTGTACACCCAGGCggcgggcg                                    >  2:671421/1‑100 (MQ=255)
                                                              cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                   >  2:323034/1‑100 (MQ=255)
                                                               aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                  >  1:17663/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                >  2:2176447/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                >  1:1088092/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                <  1:1642536/100‑1 (MQ=255)
                                                                  gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                               >  2:2337856/1‑100 (MQ=255)
                                                                   cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                              <  2:132002/100‑1 (MQ=255)
                                                                   cATATATTCCGAGAGCGGTTGTCTGGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                              <  2:1035655/100‑1 (MQ=255)
                                                                    atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                             <  2:1539173/100‑1 (MQ=255)
                                                                    atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTATACACCCAGGCGGCGGGCGACGGCtt                             >  2:2899864/1‑100 (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                         >  2:2174530/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         >  1:2037387/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         >  2:1879520/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         <  1:2358608/100‑1 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGGCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         <  1:1168530/100‑1 (MQ=255)
                                                                           ccGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTAATAGGc                      >  2:1351369/1‑100 (MQ=255)
                                                                           ccGAGAGCGGTTGTCCTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                      >  1:3629211/1‑100 (MQ=255)
                                                                            cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                     <  1:597871/100‑1 (MQ=255)
                                                                              agagCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                   >  1:1210002/1‑100 (MQ=255)
                                                                                agCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATACGCTTCGa                 <  1:2052105/100‑1 (MQ=255)
                                                                                 gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                >  1:1741444/1‑100 (MQ=255)
                                                                                  tcatcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc               >  1:3247983/7‑100 (MQ=37)
                                                                                  cGGTTGTCTTGCCGTTCCAGCTACAAGGCAAAATGTTGAACAGAAAAATCAGTCCAGCTGTAAAACCAGGCGGCGCGAGACCGCTTCATAGGCTTCg                  >  2:2970719/1‑97 (MQ=255)
                                                                                   ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa              <  1:2305809/100‑1 (MQ=255)
                                                                                   ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa              <  1:1324852/100‑1 (MQ=255)
                                                                                     ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg            <  1:2110543/100‑1 (MQ=255)
                                                                                      tGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           <  1:1422833/100‑1 (MQ=255)
                                                                                      catCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           >  1:1371416/3‑100 (MQ=38)
                                                                                        tCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa         >  2:1840532/1‑100 (MQ=39)
                                                                                         ctggCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTCTACACCCAGGCGGCGGGCGACGGCTTCATAGGCCTCGATCAGGCCa         >  2:3156613/4‑99 (MQ=255)
                                                                                         cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           <  2:2409970/97‑1 (MQ=255)
                                                                                          ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc       >  2:858063/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc       >  2:3409649/1‑99 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:3412079/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:3768944/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:767132/1‑100 (MQ=255)
                                                                                           tGCCGTGCAGCTGGCACGACAAGATAATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCATGCCAcc       >  2:951252/1‑99 (MQ=255)
                                                                                                  |                                                                                                 
TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGC  >  NC_000913/2596775‑2596969

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: