Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 16.9% 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.G16.9% 506.4 / 120.3 160intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (67/66);  new base G (12/15);  total (79/81)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

Only 100 of 165 total aligned reads displayed.

TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977
                                                                                                    |                                                                                                   
tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                      >  2:3524673/1‑100 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGTCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                                   >  1:96709/1‑100 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                    >  1:25204/1‑99 (MQ=255)
   ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCa                                                                                                  <  1:1060659/100‑1 (MQ=255)
   ttCAGAGGAAGATTACCACAGGAAAGCGGGATAAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                                   >  1:3741583/1‑99 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                                 >  1:800857/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                                 >  1:1494604/1‑100 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              <  1:3369146/100‑1 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              >  1:386676/1‑100 (MQ=255)
        gggAAGATTGCCACAGGAAAGCGGGATGAAAGGGCGGAGGAATGAGCGTTCAAGGCATAGATTCGGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                             <  2:3767080/99‑1 (MQ=255)
        agGAAGATTACCACAGGAGAGCGTGATGAAATGGCGGGGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCG‑GGCAGCTGGCACGGc                                                                                              <  2:2491587/99‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                             <  2:2147656/100‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                             <  1:3903524/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                            >  1:2242826/1‑100 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                            <  1:364917/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                            >  1:2952494/1‑100 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                            <  2:434316/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                            <  2:3621087/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                            <  1:3298281/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGTGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAg                                                                                           <  2:3017045/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                           <  2:1140621/100‑1 (MQ=39)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                           <  2:31282/100‑1 (MQ=39)
            aGATTACCACCGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAga                                                                                          >  2:2809434/1‑99 (MQ=255)
            aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                         <  1:2042129/100‑1 (MQ=255)
               ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatgat                                                                                      <  1:1875026/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAGAGGGCGGAGGAGTGAGCGGTCAAGGCAGATATGCCGAGAGCGGGTGTCTTGCCGTGCCAGCT‑GGACGGCCAGAggctga                                                                                    >  2:3608430/1‑100 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                    <  2:259812/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                    >  2:875541/1‑100 (MQ=255)
                  ctacaGGAAAGCGGGATGAAATGGCGGAGGAGTGAGCGTTCAGTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                   >  1:2666857/3‑100 (MQ=255)
                  ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                   >  1:1627294/1‑100 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAc                                                                                  >  2:563849/1‑100 (MQ=255)
                   cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAc                                                                                  <  1:2048302/100‑1 (MQ=255)
                    acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                                 >  1:1765095/1‑100 (MQ=255)
                          aaGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGCGCCTGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                           >  2:3459956/1‑100 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                         >  1:2885925/1‑100 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                       <  2:2952476/100‑1 (MQ=255)
                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                       <  2:2956487/100‑1 (MQ=255)
                               ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                                      >  2:2214694/1‑100 (MQ=255)
                               ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                                      >  1:944833/1‑100 (MQ=255)
                               ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                                      >  2:3843670/1‑100 (MQ=255)
                                  tGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                   <  1:1653387/100‑1 (MQ=255)
                                   gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGGTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                   >  2:3651220/1‑99 (MQ=255)
                                    aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                 >  1:1922855/1‑100 (MQ=255)
                                     aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                <  1:3030722/100‑1 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAACTCAGTCCAGCTGt                                                              >  2:594859/1‑100 (MQ=255)
                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCCGTCCAGCTGt                                                              >  2:595017/1‑100 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  1:2214706/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            <  1:3493945/100‑1 (MQ=255)
                                          cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  1:2791981/1‑100 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  1:3807371/100‑1 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         >  1:3569204/1‑100 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  1:121218/100‑1 (MQ=255)
                                             aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                        <  2:474644/100‑1 (MQ=255)
                                             aggaATCAGCGTTCAATGCATATATTCCGAGAGCCGTTCTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                        <  2:2567726/100‑1 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       >  1:757888/1‑100 (MQ=255)
                                               gaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAg                                                      >  1:1948457/1‑100 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 >  1:698462/1‑100 (MQ=255)
                                                      cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                               >  2:3310561/1‑100 (MQ=255)
                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                              <  1:3637574/100‑1 (MQ=255)
                                                         tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             <  1:102655/99‑1 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGGACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          <  2:129400/100‑1 (MQ=255)
                                                            aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                         <  1:518989/100‑1 (MQ=255)
                                                             tGCATCTATTCCTAGAGCGGTTTTCTTGCCGTTTCAGCT‑CCACGACAAGATGATGAACAGAAAAATCAATCCTGCTGTACATCCAGGCGGCGGGAGACgg                                        <  2:3522052/100‑1 (MQ=255)
                                                               cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                      >  2:1875633/1‑100 (MQ=255)
                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     >  2:2640200/1‑100 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTTTTGCCGTGCGAGCT‑GCACGGGCAGATGGTGAGCAGGAAAATCAGTCCAGCTGGACACCCAGGCGGCGGGGGa                                           >  2:2706536/1‑93 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTAAACCCAGGCGGCGGGCGACGGCTTc                                    >  2:1031670/1‑97 (MQ=255)
                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAACATCAGTCCAGCTGTACACCCAGTCGGCGGGCGACGGCTTCATAg                                <  2:1370387/100‑1 (MQ=255)
                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                <  2:2722672/100‑1 (MQ=255)
                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                >  2:1421151/1‑100 (MQ=255)
                                                                       ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCGGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               <  1:2646581/99‑1 (MQ=255)
                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             <  1:436769/100‑1 (MQ=255)
                                                                          agagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           <  2:2257296/100‑1 (MQ=255)
                                                                            agCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                         <  2:3623679/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCACGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  2:1886328/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  1:3083358/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  2:3225549/100‑1 (MQ=255)
                                                                              cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGacc                       >  1:3822377/1‑98 (MQ=255)
                                                                              cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       <  2:3893073/100‑1 (MQ=255)
                                                                              cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       <  1:1715934/100‑1 (MQ=255)
                                                                              cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATAAGTGCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                         >  1:2497795/1‑98 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      >  1:2055396/1‑100 (MQ=255)
                                                                               ggTTGTCTGGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACACACGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       <  2:2248215/99‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  1:1135492/100‑1 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCT‑GCACGGCAACATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGATTCGATCAgg                    <  2:3294573/100‑1 (MQ=255)
                                                                                   gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  >  1:2579561/1‑100 (MQ=255)
                                                                                   atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGTACAGAAAAATCCGTCCAGCTGTACACGCAGGCGGCGGGAGACGGCTTCATCGGCTTCGATCAg                     >  2:2043707/2‑97 (MQ=25)
                                                                                   atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  >  1:2100365/2‑100 (MQ=39)
                                                                                   atCTTGCCGTG‑CAGCTGGCACCGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCAAGCGACGGCTTCATAGGCTTCGATCAgg                    >  2:1565932/2‑98 (MQ=38)
                                                                                    tCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  1:306725/1‑100 (MQ=255)
                                                                                    tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  1:1354027/1‑100 (MQ=39)
                                                                                     cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc                >  1:99691/1‑100 (MQ=255)
                                                                                      ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  2:1068744/1‑100 (MQ=255)
                                                                                          cGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg           >  1:561489/1‑100 (MQ=255)
                                                                                           gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  2:2002743/100‑1 (MQ=255)
                                                                                           gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          >  1:153128/1‑100 (MQ=255)
                                                                                                 gCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc    >  1:1631278/1‑100 (MQ=255)
                                                                                                 gCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc    <  1:1566027/100‑1 (MQ=255)
                                                                                                  cT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:1261278/100‑1 (MQ=255)
                                                                                                   t‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCCCTGGCgg  >  1:711207/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑