Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,874 Δ1 bp 16.1% 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.1% 405.4 / 77.7 137intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (53/62);  new base . (11/11);  total (64/73)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 141 total aligned reads displayed.

GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCT  >  NC_000913/2596774‑2596972
                                                                                                   |                                                                                                    
  aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                    >  2:1705850/1‑100 (MQ=255)
  aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                    <  1:633614/100‑1 (MQ=255)
   ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                   >  1:1685050/1‑100 (MQ=255)
      ggtTCAGAGGAAGATTACCACAGGAAAGCGCGACGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                               >  1:1917050/1‑100 (MQ=255)
       gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                               >  1:2589122/1‑99 (MQ=255)
       gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                              <  2:1637175/100‑1 (MQ=255)
       gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                              <  2:2304844/100‑1 (MQ=255)
        ttCAGAGGAAGATTACCACACGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAGTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCa                                                                                             >  1:2041033/1‑100 (MQ=255)
         tCAGAGGCAGGTTACCACAGGAAAGCGGGATGAAATCGCGGAGGAATGTGTGTGCAATGCATATATTCCGGGAACGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                            <  1:441863/100‑1 (MQ=255)
          cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                           >  2:2334680/1‑100 (MQ=255)
            gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                         >  1:3908/1‑100 (MQ=255)
            gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                         >  1:949943/1‑100 (MQ=255)
             agGAAGATTTCCACAGGAACGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                        <  1:2415009/100‑1 (MQ=255)
             agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCa                                                                                        <  2:929689/100‑1 (MQ=255)
             agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                        <  2:2030778/100‑1 (MQ=255)
             agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                        <  1:2327472/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                       <  1:804755/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                       <  2:1088142/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                       <  2:2357695/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                       <  2:2541183/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                       <  2:717183/100‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGAGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCGCGGCaa                                                                                       >  1:558423/1‑100 (MQ=255)
               gAAGATTCCCACAGGAAAGCGGGATGAAATGGAGAAAGAATGAGCGTTCCATCCATATATTCCGAGAGAGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaaa                                                                                      <  1:969521/100‑2 (MQ=255)
               gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAg                                                                                      <  1:2222405/100‑1 (MQ=255)
               gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAg                                                                                      <  1:2332404/100‑1 (MQ=255)
                aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGCGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                            >  1:2722339/1‑93 (MQ=255)
                 aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                    <  2:1710253/100‑1 (MQ=255)
                 aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgac                                                                                    <  2:767666/100‑2 (MQ=39)
                   aTTACCACAGGAAAGCGGGATTCAAGGGCGGAGGAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                    <  1:1338764/98‑1 (MQ=255)
                   aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                  <  1:1614946/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                               >  2:516565/1‑100 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                               <  1:182786/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                               >  2:1723403/1‑100 (MQ=255)
                       ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                              >  1:2673889/1‑100 (MQ=255)
                       ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                              >  2:1875401/1‑100 (MQ=255)
                        cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAc                                                                             >  1:238373/1‑100 (MQ=255)
                         acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                            <  2:450234/100‑1 (MQ=255)
                                aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAt                                                                     >  2:79057/1‑100 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATTCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGTACAGAAAAATc                                                                    <  1:1529476/100‑1 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATTAACAGAAAAATc                                                                    <  2:1917059/100‑1 (MQ=255)
                                   gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                  <  2:3908/100‑1 (MQ=255)
                                   cgaaTGAAATGGCGGATGAATGAGCGTTCAACCCATCTCTTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAACTc                                                                    <  1:1357579/95‑1 (MQ=255)
                                     gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                <  2:168323/100‑1 (MQ=255)
                                      aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTcc                                                               <  2:110001/100‑1 (MQ=255)
                                        gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                             <  2:631777/100‑1 (MQ=255)
                                        gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                             >  1:2503617/1‑100 (MQ=255)
                                          aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                           <  2:1471952/100‑1 (MQ=255)
                                           aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                          <  2:1456164/100‑1 (MQ=255)
                                           aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                          <  1:538742/100‑1 (MQ=255)
                                             ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                        >  2:1050385/1‑100 (MQ=255)
                                             ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                        <  2:1310088/100‑1 (MQ=255)
                                              gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                       <  2:610398/100‑1 (MQ=255)
                                              gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                       <  2:613965/100‑1 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                      >  2:1044559/1‑100 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                      <  1:2268582/100‑1 (MQ=255)
                                                ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTacac                                                     >  1:2601302/1‑100 (MQ=255)
                                                   ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                  >  2:256465/1‑100 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAATAGAAAAATCAGGCCAGCTGTACACCCAgg                                                <  1:1195244/100‑1 (MQ=255)
                                                      aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                               >  2:635074/1‑100 (MQ=255)
                                                         aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                               >  2:875820/1‑97 (MQ=255)
                                                               cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                      <  2:2035442/100‑1 (MQ=255)
                                                               cAATGCATATATCCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                      >  1:903644/1‑100 (MQ=255)
                                                               cAATACATATATTCCGAGAGCGGTTGACTTGCCGTGCCAGCT‑GCACGGCAAGAAAAAAAAAAAAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                       <  1:2868482/99‑1 (MQ=255)
                                                                  tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                   >  2:122278/1‑100 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGTGGCGGGCGACGGCt                                 >  1:1986729/1‑100 (MQ=255)
                                                                     atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                >  1:2059699/1‑100 (MQ=255)
                                                                     atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                >  2:2687622/1‑100 (MQ=255)
                                                                      tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                               >  2:2369802/1‑100 (MQ=255)
                                                                      tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                               >  2:1417897/1‑100 (MQ=255)
                                                                      tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                               <  1:927508/100‑1 (MQ=255)
                                                                         atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGGACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                              >  1:398251/1‑98 (MQ=255)
                                                                         atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                            <  1:1901394/100‑1 (MQ=255)
                                                                          ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                           >  1:874699/1‑100 (MQ=255)
                                                                          ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                           <  1:343902/100‑1 (MQ=255)
                                                                              gagaGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                       >  1:1666613/1‑100 (MQ=255)
                                                                               agagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                      <  1:2269440/100‑1 (MQ=255)
                                                                                  gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                   >  2:1655084/1‑100 (MQ=255)
                                                                                   cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGAGGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                  <  1:1052972/100‑1 (MQ=255)
                                                                                     gTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                <  1:1713357/100‑1 (MQ=255)
                                                                                     gTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                >  2:156604/1‑100 (MQ=255)
                                                                                      ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGATGGCTTCATAGGCTTCGATCAgg               <  1:1691728/100‑1 (MQ=255)
                                                                                      ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg               <  2:1685057/100‑1 (MQ=255)
                                                                                      ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg               <  2:2794309/100‑1 (MQ=255)
                                                                                       tGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc              <  2:1299484/100‑1 (MQ=255)
                                                                                       catCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc              >  1:2398215/3‑100 (MQ=38)
                                                                                        gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc             <  2:1632600/100‑1 (MQ=255)
                                                                                        atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAACAATCAGTCCAGCTGTACACGCAGCCGGCGGGCCACGGCTTCCTATGCTTCCATCa                 >  1:1793839/2‑96 (MQ=21)
                                                                                         tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa            >  2:1634274/1‑100 (MQ=39)
                                                                                           ttGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc          >  1:1358486/1‑100 (MQ=255)
                                                                                           ttGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc          <  2:1121068/100‑1 (MQ=255)
                                                                                            tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcca         >  2:2741320/1‑99 (MQ=255)
                                                                                            tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg         >  2:2442514/1‑100 (MQ=255)
                                                                                            tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg         >  1:2094448/1‑100 (MQ=255)
                                                                                             gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa        <  2:829679/100‑1 (MQ=255)
                                                                                             gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa        >  2:975622/1‑100 (MQ=255)
                                                                                                gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc     <  2:2090040/100‑1 (MQ=255)
                                                                                                gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc     <  1:516562/100‑1 (MQ=255)
                                                                                                   ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct  <  2:2519760/100‑1 (MQ=255)
                                                                                                   ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct  <  2:238374/100‑1 (MQ=255)
                                                                                                   ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct  >  2:1166850/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCT  >  NC_000913/2596774‑2596972

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