Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,093,437 C→T 5.9% L42L (CTG→TTG)  tpm → thiopurine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,093,4370CT 5.9% 201.6 / 2.3 68L42L (CTG→TTG) tpmthiopurine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TCATCTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAGGTGCGG  >  NC_002947/2093291‑2093571
                                                                                                                                                  |                                                                                                                                      
tcatcTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGc                                                                                                                                        <  4:110450/147‑1 (MQ=255)
tcatcTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGc                                                                                                                                        >  3:110450/1‑147 (MQ=255)
    cTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCACGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGc                                                                                                                                  >  3:126960/1‑149 (MQ=255)
    cTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGc                                                                                                                                  >  8:17654/1‑149 (MQ=255)
     tAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGcc                                                                                                                                 >  6:147766/1‑149 (MQ=255)
       aaaaTCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGcc                                                                                                                                 <  7:17654/147‑1 (MQ=255)
             gggAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCTCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCTCGTGCTGGTGCCGTTGTGc                                                                                                                          <  8:364463/148‑1 (MQ=255)
                aGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCa                                                                                                                      <  8:50007/149‑1 (MQ=255)
                    aTCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAg                                                                                                                  >  7:287595/1‑149 (MQ=255)
                      cATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGcc                                                                                                                >  5:274398/1‑149 (MQ=255)
                       aTGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCt                                                                                                               <  5:147766/149‑1 (MQ=255)
                        tGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTg                                                                                                              <  3:383883/149‑1 (MQ=255)
                          gAACCAGCGTTCTGGCAGCAGCGGGGGGCCGACAACCAGAGCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGa                                                                                                            <  5:175734/149‑1 (MQ=255)
                               aGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTgg                                                                                                       >  3:14052/1‑149 (MQ=255)
                                gCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGc                                                                                                      <  7:49539/149‑1 (MQ=255)
                                        ggaggagcGGTGGGCCGACAAGCAGATCGGCTTTCACGAAGCGCAGCTGAACCCCTATCTGCAGACGTATTGGCCAGAGTTGCAGTTGGCGGCCGGGAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGc                                                                                                                 <  8:17544/125‑1 (MQ=255)
                                        gcagcaGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGAGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCCGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGc                                                                                                                 >  7:17544/1‑130 (MQ=255)
                                              gctgTGGGCCGACAACCAGATCGGCTTTCCCCAAGCCCAGGTGAACCCCTCTCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTggcagggc                                                                                        <  4:126960/146‑1 (MQ=255)
                                              gcGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTggcaggg                                                                                         >  7:303622/1‑148 (MQ=255)
                                                       cGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATc                                                                               >  3:375856/1‑149 (MQ=255)
                                                       cGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATc                                                                               >  6:188470/1‑149 (MQ=255)
                                                             ccAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGt                                                                                                                                    >  8:1321/1‑90 (MQ=255)
                                                             ccAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGt                                                                                                                                    <  7:1321/90‑1 (MQ=255)
                                                             ccAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTAt                                                                         >  5:353654/1‑149 (MQ=255)
                                                              cAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTAtt                                                                        >  5:60378/1‑149 (MQ=255)
                                                               aGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTg                                                                       >  3:204658/1‑149 (MQ=255)
                                                               aGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTg                                                                       <  6:353654/149‑1 (MQ=255)
                                                               aGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTg                                                                       >  2:15294/1‑149 (MQ=255)
                                                                gATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTgg                                                                      <  6:217254/149‑1 (MQ=255)
                                                                 aTCGTCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTggg                                                                     >  8:41790/1‑149 (MQ=255)
                                                                  tCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGgt                                                                    >  5:162151/1‑149 (MQ=255)
                                                                  tCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGgt                                                                    <  1:349073/149‑1 (MQ=255)
                                                                    ggCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGgtg                                                                   <  3:72579/148‑1 (MQ=255)
                                                                     gCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGtgtg                                                                 <  2:140770/149‑1 (MQ=255)
                                                                         tCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGc                                                             >  1:158134/1‑149 (MQ=255)
                                                                         tCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGGGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGa                                                               <  5:322378/147‑1 (MQ=255)
                                                                          cACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCGTCGCGTATTGGGTGTGGAGCt                                                            >  6:276218/1‑149 (MQ=255)
                                                                             cAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTc                                                         >  2:233585/1‑149 (MQ=255)
                                                                             cAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTc                                                         >  5:222651/1‑149 (MQ=255)
                                                                                  gcAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCgcggc                                                    >  3:235956/1‑149 (MQ=255)
                                                                                  gcAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCgcggc                                                    >  4:283599/1‑149 (MQ=255)
                                                                                       tGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCggcgggcg                                               <  4:53317/149‑1 (MQ=255)
                                                                                            cccTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGTTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTg                                            >  5:273629/1‑147 (MQ=255)
                                                                                            cccTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTgga                                          <  6:60378/149‑1 (MQ=255)
                                                                                            cccTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTg                                            >  5:131222/1‑147 (MQ=255)
                                                                                            accTAACTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCTGCAGCCGCGTGTTGCTGCCGTTGTGCGGAAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTg                                            >  5:273632/2‑147 (MQ=255)
                                                                                                 tCTGCAGACGTATTGGCCACAGGTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTggagga                                       <  6:274398/147‑1 (MQ=255)
                                                                                                     cAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCACGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGAtttttt                                 <  5:376421/149‑1 (MQ=255)
                                                                                                        aCGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATcgcg                                                                            <  1:75853/103‑1 (MQ=255)
                                                                                                        aCGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATcgcg                                                                            >  2:75853/1‑103 (MQ=255)
                                                                                                        aCGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCg                              <  7:41790/149‑1 (MQ=255)
                                                                                                           tATTGGCCACAGGTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTAGTGGGTGTGGAGCTGTCGCGGCGGGCGGGGGAGGAGGTTTTCCGt                             >  2:207764/1‑147 (MQ=255)
                                                                                                            aTTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTcc                               >  6:20002/1‑144 (MQ=255)
                                                                                                            aTTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTcc                               <  5:20002/144‑1 (MQ=255)
                                                                                                            aTTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAg                          >  2:49838/1‑149 (MQ=255)
                                                                                                             ttGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAg                          >  2:299047/1‑148 (MQ=255)
                                                                                                                gCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACg                      <  4:375856/149‑1 (MQ=255)
                                                                                                                     agttgcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCtt                 <  1:49838/149‑1 (MQ=255)
                                                                                                                     agttgcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCtt                 <  4:235956/149‑1 (MQ=255)
                                                                                                                       ttgcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCtt                 >  4:378703/1‑147 (MQ=255)
                                                                                                                       ttgcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGa               <  3:378703/149‑1 (MQ=255)
                                                                                                                       ttgcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGa               <  4:204658/149‑1 (MQ=255)
                                                                                                                         gcagttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAg              <  3:283599/148‑1 (MQ=255)
                                                                                                                           agttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGcc           >  5:373421/1‑149 (MQ=255)
                                                                                                                             ttgGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGcc           >  4:384374/1‑147 (MQ=255)
                                                                                                                              tgtcgcCCTGCAGCCGCCGGTTTGGGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTTTTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAg        <  6:273632/146‑1 (MQ=255)
                                                                                                                               ggCGCCCGGCAGCCGCGTGTTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAg        <  6:273629/148‑1 (MQ=255)
                                                                                                                                gcgcCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAGGt      <  6:162151/149‑1 (MQ=255)
                                                                                                                                gcgcCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATATTTTCCGTGAGCACGGGCTTGAGGCCGAGGt      >  8:368301/1‑149 (MQ=255)
                                                                                                                                    ccGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAGGTGCgg  <  1:53198/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
TCATCTAAAAATCGGGAGTGATCATGGAACCAGCGTTCTGGCAGCAGCGGTGGGCCGACAACCAGATCGGCTTTCACCAAGCGCAGGTGAACCCCTATCTGCAGACGTATTGGCCACAGTTGCAGTTGGCGCCCGGCAGCCGCGTGCTGGTGCCGTTGTGCGGCAAAAGCCTGGACTTGGCCTGGCTGGCAGGGCAGGGCCATCGCGTATTGGGTGTGGAGCTGTCGCGGCGGGCGGTGGAGGATTTTTTCCGTGAGCACGGGCTTGAGGCCGAGGTGCGG  >  NC_002947/2093291‑2093571

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