Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,095,997 C→T 5.5% intergenic (+68/+30) 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,9970CT 5.5% 642.0 / 17.1 198intergenic (+68/+30)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base C (92/95);  new base T (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.98e-01

Only 100 of 202 total aligned reads displayed.

GGATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095848‑2096146
                                                                                                                                                     |                                                                                                                                                     
ggATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGTGAATAAAAAAAGGGCGAcc                                                                                                                                                       <  2:3331081/150‑1 (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                                                                                                                                     >  2:65767/1‑150 (MQ=255)
                  cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtttt                                                                                                                                     >  1:2498922/1‑150 (MQ=255)
                   ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                    >  1:1477096/1‑150 (MQ=255)
                    cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                   >  2:2408557/1‑150 (MQ=255)
                     gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                  >  1:3577826/1‑150 (MQ=255)
                     gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                  <  2:260028/150‑1 (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:2202188/150‑1 (MQ=255)
                            gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                           <  2:3110080/150‑1 (MQ=255)
                              gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGAt                                                                                                                         <  1:4521584/150‑1 (MQ=255)
                              gCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGAt                                                                                                                         >  2:2115764/1‑150 (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:3181802/1‑150 (MQ=255)
                                           tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                            >  2:3181604/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)
                                                   gttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGcc                                                                                                    <  2:2911957/150‑1 (MQ=255)
                                                    ttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCt                                                                                                   <  1:719392/150‑1 (MQ=255)
                                                     tgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                  >  1:3414963/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                                                                                            <  1:3787223/150‑1 (MQ=255)
                                                                 aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                      <  2:3794037/150‑1 (MQ=255)
                                                                 aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                      <  2:1919208/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:1517414/150‑1 (MQ=255)
                                                                        cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               >  2:715332/1‑150 (MQ=255)
                                                                        cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:4490049/150‑1 (MQ=255)
                                                                        cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:3949695/150‑1 (MQ=255)
                                                                        cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                               <  1:198947/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:300214/150‑1 (MQ=255)
                                                                             gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggtg                                                                          >  2:635526/1‑150 (MQ=255)
                                                                              aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                         <  2:2148547/150‑1 (MQ=255)
                                                                                aaGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACTTCGGTCGCCCTttttttattc                                                                                                                                  >  2:1711672/1‑87 (MQ=25)
                                                                                aaGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACTTCGGTCGCCCTttttttattc                                                                                                                                  <  1:1711672/91‑5 (MQ=25)
                                                                                 aGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGCGGCtc                                                                      >  1:1952195/1‑150 (MQ=255)
                                                                                   cAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcat                                                                    <  1:2408564/150‑1 (MQ=255)
                                                                                    aGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatc                                                                   <  2:2265561/150‑1 (MQ=255)
                                                                                     tgTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCATGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                  <  2:1740447/149‑1 (MQ=255)
                                                                                     ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                  <  2:2043079/150‑1 (MQ=255)
                                                                                      gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAg                                                                 >  2:4052023/1‑150 (MQ=255)
                                                                                            cAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                           >  2:1216757/1‑150 (MQ=255)
                                                                                             aGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCg                                                          <  2:2836512/150‑1 (MQ=255)
                                                                                             aGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCg                                                          >  2:776971/1‑150 (MQ=255)
                                                                                               gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGc                                                        <  2:3432678/150‑1 (MQ=255)
                                                                                               gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGc                                                        <  2:3432278/150‑1 (MQ=255)
                                                                                                     ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                  <  2:606330/150‑1 (MQ=255)
                                                                                                     ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                  >  1:2154467/1‑150 (MQ=255)
                                                                                                      cGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCg                                                 >  2:2449355/1‑150 (MQ=255)
                                                                                                      cGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCg                                                 >  2:12998/1‑150 (MQ=255)
                                                                                                       gCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGt                                                <  1:1231685/150‑1 (MQ=255)
                                                                                                          cccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCt                                             <  2:2410204/150‑1 (MQ=255)
                                                                                                           ccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                            >  1:1293219/1‑150 (MQ=255)
                                                                                                             acaGAAGGCTGCGTAAGGCCGGGAATAAAAATAGGGCGACCGAAGTCGCCCTTTGTTTTTATCCGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCCGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                          <  2:4522395/150‑1 (MQ=255)
                                                                                                             acaGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                          <  2:158522/150‑1 (MQ=255)
                                                                                                                   ggCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                    <  2:2526891/150‑1 (MQ=255)
                                                                                                                   ggCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                    <  1:2211164/150‑1 (MQ=255)
                                                                                                                   ggCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                    >  2:808520/1‑150 (MQ=255)
                                                                                                                     cTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGCGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCa                                  <  2:1952379/150‑1 (MQ=255)
                                                                                                                          tAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                             <  2:3901405/150‑1 (MQ=255)
                                                                                                                                 gggAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                      >  1:4313591/1‑150 (MQ=255)
                                                                                                                                    aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAGTGATGCCGAAGGcc                   >  2:3912725/1‑150 (MQ=255)
                                                                                                                                        aacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa               >  2:689090/4‑150 (MQ=255)
                                                                                                                                        aaaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa               >  1:825095/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  2:2144562/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              >  1:419077/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt              <  1:2406590/150‑1 (MQ=255)
                                                                                                                                          aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  2:2254931/150‑1 (MQ=255)
                                                                                                                                          aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg             <  1:4308497/150‑1 (MQ=255)
                                                                                                                                           aaggggCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATgg            >  2:4461289/4‑150 (MQ=255)
                                                                                                                                            aaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt           >  1:3171016/1‑150 (MQ=255)
                                                                                                                                                gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc       >  2:3978951/1‑150 (MQ=255)
                                                                                                                                                gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc       >  2:1124524/1‑150 (MQ=255)
                                                                                                                                                  gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg     >  1:3316155/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)
                                                                                                                                                     |                                                                                                                                                     
GGATGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAT  >  NC_002947/2095848‑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: