Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,083 G→T 7.1% T214N (ACC→AAC)  pde ← 3',5'‑cyclic‑nucleotide phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,590,0830GT 7.1% 246.9 / 6.2 79T214N (ACC→AAC) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base T (3/3);  total (38/41)
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.70e-01

GAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATCC  >  NC_002947/5589941‑5590230
                                                                                                                                                |                                                                                                                                                   
gAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGGCGACg                                                                                                                                               <  2:151990/149‑1 (MQ=255)
    tCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGt                                                                                                                                           <  6:313713/149‑1 (MQ=255)
    gcGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGTCGCAACCTGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcggcgc                                                                                                                                           <  2:29607/148‑2 (MQ=255)
     cGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTg                                                                                                                                          >  1:242068/1‑149 (MQ=255)
      gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTgg                                                                                                                                         >  5:410504/1‑149 (MQ=255)
      gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTg                                                                                                                                          <  7:245254/148‑1 (MQ=255)
        gCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGTGCCGTCGACCTGCATTCAGTTCgcggcgcgcag                                                                                                                                       <  6:129832/149‑6 (MQ=255)
        gCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGc                                                                                                                                        >  4:51591/1‑148 (MQ=255)
           gCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCaa                                                                                                                                     >  6:309458/1‑148 (MQ=255)
            cACCTTGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAAc                                                                                                                                    <  5:218606/148‑1 (MQ=255)
            cACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCaa                                                                                                                                     <  8:20433/147‑1 (MQ=255)
              ccTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCTCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAg                                                                                                                                  >  1:59122/1‑149 (MQ=255)
               cTCGAAATCCAGTGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                                >  1:198601/1‑149 (MQ=255)
                      tCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGccc                                                                                                                         >  2:316052/1‑149 (MQ=255)
                        cAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCct                                                                                                                       >  8:133326/1‑149 (MQ=255)
                           cgcCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCtctc                                                                                                                    <  5:248281/149‑1 (MQ=255)
                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                                <  6:396312/134‑1 (MQ=255)
                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                                >  5:396312/1‑134 (MQ=255)
                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                                 >  5:321665/1‑149 (MQ=255)
                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGAGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                                 >  6:158854/1‑149 (MQ=255)
                               cGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcgg                                                                                                                                               >  8:21670/1‑116 (MQ=255)
                               cGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCTAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCCGCTCTCGAAc                                                                                                                <  2:242068/149‑1 (MQ=255)
                               cGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAAc                                                                                                                <  5:184572/149‑1 (MQ=255)
                                gCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAAcc                                                                                                               >  7:304210/1‑149 (MQ=255)
                                    cACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCa                                                                                                            <  7:74981/148‑1 (MQ=255)
                                      ccccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGCACCtcctcc                                                                                                         >  7:315233/1‑149 (MQ=255)
                                      ccccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                                         <  8:343285/149‑1 (MQ=255)
                                      ccccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                                         >  3:205784/1‑149 (MQ=255)
                                      ccccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCACCCCTCTCGAACCTCATcc                                                                                                         <  3:51591/149‑1 (MQ=255)
                                       cccGGTCTCCAGCCGCCCATCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                                         >  3:124757/1‑148 (MQ=255)
                                        ccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCa                                                                                                       >  8:86504/1‑149 (MQ=255)
                                           gTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTc                                                                                                    <  8:304210/149‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCTCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                               <  2:382061/149‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCt                                                                                                  <  5:309458/146‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                                <  6:234869/148‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                                <  6:149290/148‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                               <  6:321665/149‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                               <  6:264092/149‑1 (MQ=255)
                                                 cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGAGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                                <  5:158854/148‑1 (MQ=255)
                                                       cccGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCTAACTGAATGCAGGTCGACGGCGTGGCTAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                                         <  1:316052/149‑1 (MQ=255)
                                                       cccGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                                         >  3:208436/1‑149 (MQ=255)
                                                             gtcATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGccca                                                                                    <  8:220704/146‑2 (MQ=255)
                                                               cATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCAt                                                                                 <  3:161445/149‑1 (MQ=255)
                                                                            ccAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTa                                                                    >  5:222180/1‑149 (MQ=255)
                                                                                cGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCtg                                                                <  2:59122/149‑1 (MQ=255)
                                                                                   tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                             <  4:205784/149‑1 (MQ=255)
                                                                                   tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                             <  4:208436/149‑1 (MQ=255)
                                                                                   tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                             <  4:124757/149‑1 (MQ=255)
                                                                                     ccccGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGa                                                            >  8:20264/1‑148 (MQ=255)
                                                                                       ccGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGatat                                                         <  8:118201/149‑1 (MQ=255)
                                                                                             gCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTc                                                    >  1:235246/1‑148 (MQ=255)
                                                                                              ctctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCaa                                                  <  2:198601/149‑1 (MQ=255)
                                                                                                ctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAAc                                                 >  1:344931/1‑148 (MQ=255)
                                                                                                      tCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATc                                           >  2:129716/1‑148 (MQ=255)
                                                                                                      tCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATCCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATc                                           >  8:46349/1‑148 (MQ=255)
                                                                                                      tCACCTGGAAGTCTTCGCTGCCCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGCCGCCGCCCCTCTCGAACCACATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGataac                                           >  4:171527/1‑146 (MQ=255)
                                                                                                        aCCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGaa                                        >  3:309326/1‑149 (MQ=255)
                                                                                                          cTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACa                                      <  6:222180/149‑1 (MQ=255)
                                                                                                           tGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAg                                     >  4:159552/1‑149 (MQ=255)
                                                                                                            ggAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCAACTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCCTAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAgc                                    >  7:200640/1‑149 (MQ=255)
                                                                                                                     cGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTgg                              <  7:20264/146‑1 (MQ=255)
                                                                                                                       gtGCGCGGCGCGTACTGAATGCAGGTCGACGGCGTGGGCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCTAAGCCCGTACCTGTGGATATCGTTCAACGATATCGAACAGCGACTGGGCAtt                         <  3:171527/148‑1 (MQ=255)
                                                                                                                       cTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt                         <  7:86504/149‑1 (MQ=255)
                                                                                                                       cTGCGCGCCGCGAACTGAATGCAGGTCGAAGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt                         <  7:46349/149‑1 (MQ=255)
                                                                                                                       atGCGCGCCGCTAACTGAATGCAGTGCGACTGCGTGGCCACCAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt                         <  8:315233/148‑1 (MQ=255)
                                                                                                                        tGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTg                        >  4:194480/1‑149 (MQ=255)
                                                                                                                         gcgcgcCGCGAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCTCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgc                       <  4:101477/149‑1 (MQ=255)
                                                                                                                          cgcgcCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg                      >  2:241172/1‑149 (MQ=255)
                                                                                                                           acgccgtcgACCTGCATTCAGTTCGCGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg                      >  4:131139/8‑148 (MQ=255)
                                                                                                                           acgccgtcgACCTGCATTCAGTTCGCGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg                      >  3:290730/8‑148 (MQ=255)
                                                                                                                               cgcgAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGc                 <  1:129716/149‑1 (MQ=255)
                                                                                                                                    aaCTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAt             <  6:89947/149‑1 (MQ=255)
                                                                                                                                    aaCTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAt             >  5:89947/1‑149 (MQ=255)
                                                                                                                                        gAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGgc         >  8:119942/1‑149 (MQ=255)
                                                                                                                                          aTGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCAACCCATTCCTGATGCACAAGCCCCCATAGCAAAGCCCGTACCTGTGGAAATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGgcgc       >  5:272776/1‑149 (MQ=255)
                                                                                                                                           tGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                                         >  7:373634/1‑66 (MQ=255)
                                                                                                                                           tGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                                         <  8:373634/66‑1 (MQ=255)
                                                                                                                                           tGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGgcgc       >  6:40320/1‑148 (MQ=255)
                                                                                                                                              aGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATc   >  3:109343/1‑149 (MQ=255)
                                                                                                                                               ggTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATcc  <  2:247701/149‑1 (MQ=255)
                                                                                                                                                gTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCaa                   <  3:386283/131‑1 (MQ=255)
                                                                                                                                                gTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCaa                   >  4:386283/1‑131 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATCC  >  NC_002947/5589941‑5590230

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: