Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,095,998 G→T 5.3% intergenic (+69/+29) tpm → / ← htpX thiopurine S‑methyltransferase/protease HtpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,095,9980GT 5.3% 691.9 / 18.1 200intergenic (+69/+29)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base G (93/96);  new base T (6/5);  total (99/101)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

GATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095849‑2096146
                                                                                                                                                     |                                                                                                                                                    
gATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCg                                                                                                                                                      >  2:3016194/1‑150 (MQ=255)
     gATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTc                                                                                                                                                 >  2:2415028/1‑150 (MQ=255)
      aTTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCg                                                                                                                                                >  1:755971/1‑150 (MQ=255)
      aTTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCg                                                                                                                                                >  2:2406590/1‑150 (MQ=255)
          aGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAACAAAGGGCGACTTCGGTCGCCCt                                                                                                                                            <  1:1963788/150‑1 (MQ=255)
          aGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCt                                                                                                                                            <  1:1940823/150‑1 (MQ=255)
            gACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACTTCGGTCGCCCttt                                                                                                                                          <  1:2799618/150‑1 (MQ=255)
             aCAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTg                                                                                                                                         <  2:2128916/150‑1 (MQ=255)
               aGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtt                                                                                                                                       <  1:1785592/150‑1 (MQ=255)
                gCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGttt                                                                                                                                      <  1:4303624/150‑1 (MQ=255)
                 cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtttt                                                                                                                                     >  1:1268277/1‑150 (MQ=255)
                 cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtttt                                                                                                                                     >  1:2498922/1‑150 (MQ=255)
                 cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtttt                                                                                                                                     >  2:65767/1‑150 (MQ=255)
                  ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                    >  1:1477096/1‑150 (MQ=255)
                   cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                   >  2:2408564/1‑150 (MQ=255)
                    gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                  >  1:3577826/1‑150 (MQ=255)
                      aGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCa                                                                                                                                <  1:1526663/150‑1 (MQ=255)
                        tttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAgg                                                                                                                              <  2:4088969/150‑1 (MQ=255)
                         ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                             <  1:981743/150‑1 (MQ=255)
                           gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                           <  2:2761048/150‑1 (MQ=255)
                             gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGAt                                                                                                                         <  1:4521584/150‑1 (MQ=255)
                             gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTCGTTTTTATCAGGTGGAt                                                                                                                         >  2:3519756/1‑150 (MQ=255)
                              cAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGCGCGACCGAAGTCGCCCTTTGTTATTATCAGGTGGATc                                                                                                                        >  1:1634142/1‑150 (MQ=255)
                                 gCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGc                                                                                                                     >  1:31192/1‑150 (MQ=255)
                                    ggggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcgc                                                                                                                  <  1:3036022/150‑1 (MQ=255)
                                      ggCGTCGAGTGTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCgcgc                                                                                                                <  2:390800/150‑1 (MQ=255)
                                          tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:1985191/1‑150 (MQ=255)
                                          tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:3181766/1‑150 (MQ=255)
                                            gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCAGCGCTGACGc                                                                                                          >  1:4002011/1‑150 (MQ=255)
                                              gTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAg                                                                                                        <  1:262182/150‑1 (MQ=255)
                                               tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAgg                                                                                                       >  2:3111326/1‑150 (MQ=255)
                                                   ttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCt                                                                                                   <  1:719392/150‑1 (MQ=255)
                                                    tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  >  1:3414963/1‑150 (MQ=255)
                                                    tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  >  2:300089/1‑150 (MQ=255)
                                                    tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  <  1:564714/150‑1 (MQ=255)
                                                     ggAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCg                                                                                                 >  1:4068886/1‑150 (MQ=255)
                                                       aGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGAt                                                                                               >  1:2526891/1‑150 (MQ=255)
                                                        gCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATg                                                                                              <  2:4414407/150‑1 (MQ=255)
                                                        gCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATg                                                                                              <  1:509314/150‑1 (MQ=255)
                                                          gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                            <  2:2593148/150‑1 (MQ=255)
                                                          gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                            <  1:3787223/150‑1 (MQ=255)
                                                                aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                      <  2:3794037/150‑1 (MQ=255)
                                                                 tCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTcc                                                                                     >  2:4389687/1‑150 (MQ=255)
                                                                   aGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAg                                                                                   >  1:633732/1‑150 (MQ=255)
                                                                    ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                  <  1:3703890/150‑1 (MQ=255)
                                                                       cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               >  2:715332/1‑150 (MQ=255)
                                                                       cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:198947/150‑1 (MQ=255)
                                                                       cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:3949695/150‑1 (MQ=255)
                                                                       cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:4490049/150‑1 (MQ=255)
                                                                       cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCGTAGTTTTTATCAGGTGGATCAGCCGCGCAGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               >  2:1110153/1‑150 (MQ=255)
                                                                           cGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggt                                                                           >  1:1372402/1‑150 (MQ=255)
                                                                           cGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggt                                                                           >  2:3918956/1‑150 (MQ=255)
                                                                            gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggtg                                                                          >  2:635526/1‑150 (MQ=255)
                                                                             aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                         <  2:2148495/150‑1 (MQ=255)
                                                                               aaGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACTTCGGTCGCCCTttttttattc                                                                                                                                  <  1:1711672/91‑5 (MQ=25)
                                                                                aGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtc                                                                      >  2:1850667/1‑150 (MQ=255)
                                                                                aGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGCGGCtc                                                                      >  1:1952195/1‑150 (MQ=255)
                                                                                 gCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCt                                                                                                                                            >  2:1658858/1‑79 (MQ=255)
                                                                                  cAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcat                                                                    <  1:2408564/150‑1 (MQ=255)
                                                                                     gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAg                                                                 >  2:2229214/1‑150 (MQ=255)
                                                                                     gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAg                                                                 >  2:4052023/1‑150 (MQ=255)
                                                                                           cAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                           >  2:1216757/1‑150 (MQ=255)
                                                                                            aGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCg                                                          >  2:776971/1‑150 (MQ=255)
                                                                                              gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGc                                                        <  2:3432278/150‑1 (MQ=255)
                                                                                              gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGc                                                        <  2:3432678/150‑1 (MQ=255)
                                                                                                    ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                  <  2:606330/150‑1 (MQ=255)
                                                                                                    ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                  >  1:2154467/1‑150 (MQ=255)
                                                                                                     cGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCg                                                 >  2:12998/1‑150 (MQ=255)
                                                                                                      gCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGt                                                <  1:1231685/150‑1 (MQ=255)
                                                                                                          ccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                            >  1:1293219/1‑150 (MQ=255)
                                                                                                            acaGAAGGCTGCGTAAGGCCGGGAATAAAAATAGGGCGACCGAAGTCGCCCTTTGTTTTTATCCGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCCGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                          <  2:4522395/150‑1 (MQ=255)
                                                                                                            acaGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                          >  2:3191160/1‑150 (MQ=255)
                                                                                                                  ggCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                    >  2:808520/1‑150 (MQ=255)
                                                                                                                  ggCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                    <  1:2211164/150‑1 (MQ=255)
                                                                                                                   gCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgcc                                   >  2:155672/1‑150 (MQ=255)
                                                                                                                    cTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGCGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCa                                  <  2:1952195/150‑1 (MQ=255)
                                                                                                                        gTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGa                              >  2:2198745/1‑150 (MQ=255)
                                                                                                                        gTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGa                              >  2:2841165/1‑150 (MQ=255)
                                                                                                                         tAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                             <  2:3901405/150‑1 (MQ=255)
                                                                                                                               cGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGaa                       <  2:1724941/150‑1 (MQ=255)
                                                                                                                                gggAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                      >  1:4313591/1‑150 (MQ=255)
                                                                                                                                   aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAGTGATGCCGAAGGcc                   >  2:3912725/1‑150 (MQ=255)
                                                                                                                                     tAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCtt                 <  2:2492074/150‑1 (MQ=255)
                                                                                                                                      aaacaaaGGGCGACTTCGGTCGCCCTTTTTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTc                >  2:3084055/1‑150 (MQ=255)
                                                                                                                                       aacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa               >  2:689090/4‑150 (MQ=255)
                                                                                                                                       aaaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa               >  1:825095/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  1:419077/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              <  1:2406590/150‑1 (MQ=255)
                                                                                                                                         aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  1:4308497/150‑1 (MQ=255)
                                                                                                                                         aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  2:2047345/150‑1 (MQ=255)
                                                                                                                                          aaggggCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATgg            >  2:4461289/4‑150 (MQ=255)
                                                                                                                                           aaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt           >  1:3171016/1‑150 (MQ=255)
                                                                                                                                           aaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGACCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt           >  2:2260645/1‑150 (MQ=255)
                                                                                                                                               gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc       >  2:1124524/1‑150 (MQ=255)
                                                                                                                                               gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc       >  2:3978951/1‑150 (MQ=255)
                                                                                                                                                 gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg     >  1:3316155/1‑150 (MQ=255)
                                                                                                                                                  aCCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc    >  2:1377599/1‑150 (MQ=255)
                                                                                                                                                  aCCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc    >  1:611893/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCa   >  2:4134679/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCa   >  2:879429/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095849‑2096146

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