Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,095,999 A→C 5.6% intergenic (+70/+28) 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,9990AC 5.6% 690.6 / 20.9 198intergenic (+70/+28)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base A (92/95);  new base C (6/5);  total (98/100)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 200 total aligned reads displayed.

TGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095851‑2096146
                                                                                                                                                    |                                                                                                                                                   
tgtgATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGaa                                                                                                                                                    <  2:701756/150‑1 (MQ=255)
   gATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTc                                                                                                                                                 >  2:2415028/1‑150 (MQ=255)
    aTTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCg                                                                                                                                                >  2:2406590/1‑150 (MQ=255)
    aTTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCg                                                                                                                                                >  1:755971/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)
             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)
                ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                    >  1:1477096/1‑150 (MQ=255)
                 cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGACTAAAAAAAGGTCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                   >  2:820617/1‑150 (MQ=255)
                 cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                   >  2:2408564/1‑150 (MQ=255)
                  gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                  >  1:3577826/1‑150 (MQ=255)
                   aaGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATc                                                                                                                                 >  2:4444864/1‑150 (MQ=255)
                    aGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCa                                                                                                                                <  1:1526663/150‑1 (MQ=255)
                       ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                             <  1:981743/150‑1 (MQ=255)
                        tgtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTg                                                                                                                            >  2:4521584/1‑150 (MQ=255)
                         gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                           <  2:2761048/150‑1 (MQ=255)
                         gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                           <  2:3110080/150‑1 (MQ=255)
                          tgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGa                                                                                                                          >  2:4308497/1‑150 (MQ=255)
                           gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGAt                                                                                                                         >  2:2115764/1‑150 (MQ=255)
                           gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGAt                                                                                                                         <  1:4521584/150‑1 (MQ=255)
                            cAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGCGCGACCGAAGTCGCCCTTTGTTATTATCAGGTGGATc                                                                                                                        >  1:1634142/1‑150 (MQ=255)
                               gCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGc                                                                                                                     >  1:31192/1‑150 (MQ=255)
                                  ggggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcgc                                                                                                                  <  1:3036022/150‑1 (MQ=255)
                                   gggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcgcg                                                                                                                 >  2:4054921/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:1985191/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:3181604/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:3181802/1‑150 (MQ=255)
                                          gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCAGCGCTGACGc                                                                                                          >  1:4002011/1‑150 (MQ=255)
                                            gTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAg                                                                                                        <  1:262182/150‑1 (MQ=255)
                                             tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAgg                                                                                                       <  2:3065494/150‑1 (MQ=255)
                                                 ttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCt                                                                                                   <  1:719392/150‑1 (MQ=255)
                                                  tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  >  2:300089/1‑150 (MQ=255)
                                                  tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  <  2:3456816/150‑1 (MQ=255)
                                                  tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  <  1:564714/150‑1 (MQ=255)
                                                  tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  >  1:3414963/1‑150 (MQ=255)
                                                   ggAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCg                                                                                                 >  1:4068886/1‑150 (MQ=255)
                                                     aGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGAt                                                                                               >  1:2526891/1‑150 (MQ=255)
                                                      gCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATg                                                                                              <  1:509314/150‑1 (MQ=255)
                                                        gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                            <  2:2593148/150‑1 (MQ=255)
                                                        gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                            <  1:3802263/150‑1 (MQ=255)
                                                        gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                            <  1:3787223/150‑1 (MQ=255)
                                                          gtgtATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGt                                                                                          >  2:4410234/1‑150 (MQ=255)
                                                                 aGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAg                                                                                   >  1:633732/1‑150 (MQ=255)
                                                                  ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                  <  1:3703890/150‑1 (MQ=255)
                                                                     cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:198947/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)
                                                                          gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggtg                                                                          >  2:791075/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)
                                                                                  ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                  <  2:3607074/150‑1 (MQ=255)
                                                                                  ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                  <  2:2043079/150‑1 (MQ=255)
                                                                                   gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAg                                                                 >  2:2229214/1‑150 (MQ=255)
                                                                                         cAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                           >  2:1216757/1‑150 (MQ=255)
                                                                                            gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGc                                                        <  2:3432531/150‑1 (MQ=255)
                                                                                                  ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                  >  1:2154467/1‑150 (MQ=255)
                                                                                                   cGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCg                                                 >  2:12998/1‑150 (MQ=255)
                                                                                                   cGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCAAAGccc                                                 <  2:688127/150‑2 (MQ=255)
                                                                                                    gCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGt                                                <  1:1231685/150‑1 (MQ=255)
                                                                                                        ccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                            >  1:1293219/1‑150 (MQ=255)
                                                                                                        ccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                            >  2:3339231/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:2841165/1‑150 (MQ=255)
                                                                                                                      gTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGa                              >  2:2198745/1‑150 (MQ=255)
                                                                                                                       tAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                             <  2:3901618/150‑1 (MQ=255)
                                                                                                                         aGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGc                           <  2:651797/150‑1 (MQ=255)
                                                                                                                             cGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGaa                       <  2:1724941/150‑1 (MQ=255)
                                                                                                                              gggAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                      >  1:4313591/1‑150 (MQ=255)
                                                                                                                                 aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGcc                   <  2:3836051/150‑1 (MQ=255)
                                                                                                                                 aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGcc                   >  2:3912690/1‑150 (MQ=255)
                                                                                                                                 aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGcc                   >  2:3912749/1‑150 (MQ=255)
                                                                                                                                   tAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCtt                 <  2:2492074/150‑1 (MQ=255)
                                                                                                                                    aaacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTc                >  2:4025939/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa               >  1:825095/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              <  2:896742/150‑1 (MQ=255)
                                                                                                                                      aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  2:581703/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  2:4122766/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  1:419077/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              <  1:2406590/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  2:3941168/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  1:4308497/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  2:727032/150‑1 (MQ=255)
                                                                                                                                         aaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt           >  1:3171016/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGACCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt           >  2:2260645/1‑150 (MQ=255)
                                                                                                                                             gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc       >  2:1124524/1‑150 (MQ=255)
                                                                                                                                               gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg     >  1:3316155/1‑150 (MQ=255)
                                                                                                                                                aCCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc    >  2:1377599/1‑150 (MQ=255)
                                                                                                                                                aCCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc    <  2:61193/150‑1 (MQ=255)
                                                                                                                                                aCCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc    >  1:611893/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095851‑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: