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.

TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGC  >  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)
     ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                            <  1:1612470/100‑1 (MQ=255)
      gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                           >  2:1789888/1‑100 (MQ=255)
       ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCa                                                                                          >  1:2593353/1‑100 (MQ=255)
        tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGCGGAATGAGCGTTCAATGCATATATTACAAAAGAGGTTGTCTTGCAGTACCAGCT‑GCa                                                                                          >  1:1567600/1‑99 (MQ=255)
        tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                         >  1:1035666/1‑100 (MQ=255)
        tCAGAGGAAGATTACCACAGGAAAGAGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                         >  2:420374/1‑100 (MQ=255)
          agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                       <  2:3025142/100‑1 (MQ=255)
          agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCCGAGGAATGAGCGCTCACTGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCT‑GCActg                                                                                       <  2:2733397/100‑3 (MQ=255)
           gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                      >  1:914912/1‑100 (MQ=255)
           gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATAATCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGc                                                                                      <  2:2088874/100‑1 (MQ=255)
            agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGTAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                     <  1:907107/100‑1 (MQ=255)
            agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                     <  1:3659331/100‑1 (MQ=255)
            agGAAGATTACCACAGAAAAGCGGGATGAAATGGCGGAGGAATAAGCGTACAATGCATATATTGCGAGCGCGGTTGTCTTGCCGTGCCAGCT‑GCaaagca                                                                                     >  2:2546260/1‑95 (MQ=255)
             ggAAGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCa                                                                                     <  1:1702887/99‑1 (MQ=255)
             ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                    <  2:3362360/100‑1 (MQ=255)
              gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                   <  2:1063221/100‑1 (MQ=39)
              gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                   <  2:430576/100‑1 (MQ=39)
               aaGATTACCACAGGAAAGCGGTATGAAATGGCGGAGGAATTAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                    <  1:157510/98‑1 (MQ=38)
                aGATTACCACAGGAAAGCGGGATGAAATGGCTGAGGAATGAGTGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgac                                                                                 <  1:425757/100‑2 (MQ=38)
                aGATTACCACAGGAAAGCGGGATGAAATGGCGGATGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                 <  2:3584958/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGTAATGTGTGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                             <  1:2025990/100‑1 (MQ=255)
                    tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                             <  1:2862420/100‑1 (MQ=255)
                     aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgagatcgg                                                                            >  2:1625400/1‑94 (MQ=255)
                      ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                           >  1:1328141/1‑100 (MQ=255)
                      ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                           >  1:2979088/1‑100 (MQ=255)
                       cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCGTATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGTT‑GCACGGCAAGatgatgn                                                                            <  2:1472361/98‑2 (MQ=255)
                       cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAc                                                                          >  2:2092686/1‑100 (MQ=255)
                        acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                         >  1:2805414/1‑100 (MQ=255)
                           ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaa                                                                      >  2:1960488/1‑100 (MQ=255)
                            gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaa                                                                     >  1:1784522/1‑100 (MQ=255)
                                  gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                 <  2:827157/98‑1 (MQ=255)
                                  gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GAACGGCAAGATGATGAACAGAAAAATCAg                                                               >  1:1896557/1‑100 (MQ=255)
                                     aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGAT‑GCACGGCAAGATGATGAACAGAAAAATCAGTcc                                                            >  2:2052087/1‑100 (MQ=255)
                                      tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                           <  1:3928108/100‑1 (MQ=255)
                                      tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                           >  2:3724762/1‑100 (MQ=255)
                                        aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                         <  2:2157249/100‑1 (MQ=255)
                                          aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                       <  1:1680002/100‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCTAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                      <  1:2256755/100‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                       <  1:643120/99‑1 (MQ=255)
                                           tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                      >  2:3716152/1‑100 (MQ=255)
                                             gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                    >  1:3345634/1‑100 (MQ=255)
                                              cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GAACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                   >  1:2115685/1‑100 (MQ=255)
                                               ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTacac                                                  <  2:2992326/100‑1 (MQ=255)
                                                  ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                               <  1:1585025/100‑1 (MQ=255)
                                                     aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                            <  1:2618391/100‑1 (MQ=255)
                                                     aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                            >  1:2523336/1‑100 (MQ=255)
                                                      tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                           >  2:594369/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                          >  2:3426859/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                          >  2:1324838/1‑100 (MQ=255)
                                                       gAGCGTTCAATGCATATATTCCGACAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACACAggcgg                                          >  2:1750816/1‑100 (MQ=255)
                                                        aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                         >  1:3141658/1‑100 (MQ=255)
                                                         gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGgcggcg                                        <  2:3806133/100‑1 (MQ=255)
                                                         gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCATGCCAGCT‑GCACGGCAAAATGATGAACAGAAAAATAAAAACAGATATACACACAggcag                                          >  1:1527994/1‑96 (MQ=255)
                                                           gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                      <  1:2159204/100‑1 (MQ=255)
                                                           gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                      <  1:1464725/100‑1 (MQ=255)
                                                             tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑AAACGGCAAGATGATGAACAGAAAAAACAGTCCAGCTGTACACCCAGGCggcgggcg                                    >  2:671421/1‑100 (MQ=255)
                                                              cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                   >  2:323034/1‑100 (MQ=255)
                                                               aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                  >  1:17663/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                >  2:2176447/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                >  1:1088092/1‑100 (MQ=255)
                                                                 tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                <  1:1642536/100‑1 (MQ=255)
                                                                  gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                               >  2:2337856/1‑100 (MQ=255)
                                                                   cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                              <  2:132002/100‑1 (MQ=255)
                                                                   cATATATTCCGAGAGCGGTTGTCTGGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                              <  2:1035655/100‑1 (MQ=255)
                                                                    atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                             <  2:1539173/100‑1 (MQ=255)
                                                                    atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTATACACCCAGGCGGCGGGCGACGGCtt                             >  2:2899864/1‑100 (MQ=255)
                                                                     tatatTCCGAGTGCGGTTGTCTTGCCGTTGCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACTCAGTCGGCGGGCGACGGCtt                             <  2:3109334/99‑1 (MQ=255)
                                                                     tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                            <  1:368478/100‑1 (MQ=255)
                                                                       tatTCCTAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                          <  2:3789894/100‑1 (MQ=255)
                                                                       tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                          <  1:2956207/100‑1 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         >  2:2174530/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         >  1:2037387/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         >  2:1879520/1‑100 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         <  1:2358608/100‑1 (MQ=255)
                                                                        atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGGCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                         <  1:1168530/100‑1 (MQ=255)
                                                                           ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTAATAGGc                      >  2:1351369/1‑100 (MQ=255)
                                                                           ccGAGAGCGGTTGTCCTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                      >  1:3629211/1‑100 (MQ=255)
                                                                            cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                     <  1:597871/100‑1 (MQ=255)
                                                                              agagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                   >  1:1210002/1‑100 (MQ=255)
                                                                                agCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATACGCTTCGa                 <  1:2052105/100‑1 (MQ=255)
                                                                                 gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                >  1:1741444/1‑100 (MQ=255)
                                                                                  tcatcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc               >  1:3247983/7‑100 (MQ=37)
                                                                                  cGGTTGTCTTGCCGTTCCAGCT‑ACAAGGCAAAATGTTGAACAGAAAAATCAGTCCAGCTGTAAAACCAGGCGGCGCGAGACCGCTTCATAGGCTTCg                  >  2:2970719/1‑97 (MQ=255)
                                                                                   ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa              <  1:2305809/100‑1 (MQ=255)
                                                                                   ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa              <  1:1324852/100‑1 (MQ=255)
                                                                                     ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg            <  1:2110543/100‑1 (MQ=255)
                                                                                      tGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           <  1:1422833/100‑1 (MQ=255)
                                                                                      catCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           >  1:1371416/3‑100 (MQ=38)
                                                                                        tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa         >  2:1840532/1‑100 (MQ=39)
                                                                                         ctggCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTCTACACCCAGGCGGCGGGCGACGGCTTCATAGGCCTCGATCAGGCCa         >  2:3156613/4‑99 (MQ=255)
                                                                                         cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc           <  2:2409970/97‑1 (MQ=255)
                                                                                          ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc       >  2:858063/1‑100 (MQ=255)
                                                                                           tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc       >  2:3409649/1‑99 (MQ=255)
                                                                                           tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:3412079/1‑100 (MQ=255)
                                                                                           tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:3768944/1‑100 (MQ=255)
                                                                                           tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg      >  1:767132/1‑100 (MQ=255)
                                                                                           tGCCGTG‑CAGCTGGCACGACAAGATAATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCATGCCAcc       >  2:951252/1‑99 (MQ=255)
                                                                                                  |                                                                                                 
TAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGC  >  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: ‑