Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,079 G→T 6.3% C215* (TGC→TGA 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,0790GT 6.3% 254.8 / 6.2 80C215* (TGC→TGApde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (35/40);  new base T (3/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGC  >  NC_002947/5589934‑5590225
                                                                                                                                                   |                                                                                                                                                  
ggCTTTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGc                                                                                                                                                   <  4:109103/147‑1 (MQ=255)
ggCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAg                                                                                                                                                 >  5:234869/1‑149 (MQ=255)
       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                                                                                                    >  3:205784/1‑149 (MQ=255)
                                             ccccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                                    <  8:343285/149‑1 (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:124757/149‑1 (MQ=255)
                                                                                          tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                        <  4:205784/149‑1 (MQ=255)
                                                                                          tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                        <  4:208436/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)
                                                                                                                                                   |                                                                                                                                                  
GGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGC  >  NC_002947/5589934‑5590225

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: