Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,068:1 +T 6.8% coding (656/801 nt) 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,0681.T 6.8% 288.4 / 16.1 73coding (656/801 nt)pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base . (42/26);  new base T (2/3);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CTCAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAG  >  NC_002947/5589922‑5590214
                                                                                                                                                   |                                                                                                                                                  
cTCAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGcgccgc                                                                                                                                                   >  8:19891/1‑148 (MQ=255)
cTCAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCcgcg                                                                                                                                                  <  6:329345/149‑1 (MQ=255)
   aaTACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAAc                                                                                                                                               >  6:74045/1‑149 (MQ=255)
   aaTACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAAc                                                                                                                                               >  6:335098/1‑149 (MQ=255)
        ccTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAAt                                                                                                                                          >  5:385145/1‑149 (MQ=255)
            ggCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGc                                                                                                                                        <  7:4221/147‑1 (MQ=255)
            ggCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAg                                                                                                                                      >  4:332527/1‑149 (MQ=255)
             gCTGTCGAAGTCGAGCTGCACCTAAAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAgg                                                                                                                                     >  7:246812/1‑149 (MQ=255)
               tGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTc                                                                                                                                   >  1:142439/1‑149 (MQ=255)
               tGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTc                                                                                                                                   >  1:380428/1‑149 (MQ=255)
                 tCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTc                                                                                                                                   <  8:315854/147‑1 (MQ=255)
                 tCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGa                                                                                                                                 >  8:171501/1‑149 (MQ=255)
                   gAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGAc                                                                                                                                <  4:18437/148‑1 (MQ=255)
                         gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTgg                                                                                                                         <  2:150511/149‑1 (MQ=255)
                         gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTgg                                                                                                                         >  6:321860/1‑149 (MQ=255)
                         gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGGCGACGGCGTgg                                                                                                                         <  4:78873/149‑1 (MQ=255)
                          aGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcggcgcgca                                                                                                                        <  4:35193/149‑5 (MQ=255)
                               cACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAAc                                                                                                                    <  2:142439/148‑1 (MQ=255)
                               cACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACa                                                                                                                   <  5:74045/149‑1 (MQ=255)
                                aCCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTg                                                                                                                          <  5:361737/141‑1 (MQ=255)
                                aCCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTg                                                                                                                          >  6:361737/1‑141 (MQ=255)
                                  cTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgc                                                                                                                                 >  3:247743/1‑132 (MQ=255)
                                  cTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgc                                                                                                                                 <  4:247743/133‑2 (MQ=255)
                                  cTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcggcgcgcag                                                                                                                       <  8:22395/142‑6 (MQ=38)
                                  cTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcggcgcgcag                                                                                                                       >  7:22395/1‑137 (MQ=38)
                                   tCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgcg                                                                                                               >  7:27110/1‑149 (MQ=255)
                                       aaTCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGcgccgc                                                                                                           >  4:50917/1‑149 (MQ=255)
                                        aTCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCgccgcc                                                                                                          >  2:84546/1‑149 (MQ=255)
                                           cAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCct                                                                                                       >  5:7238/1‑149 (MQ=255)
                                            aGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCctc                                                                                                      >  3:308820/1‑149 (MQ=255)
                                             gcgcCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACACCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCctct                                                                                                     >  1:113241/1‑149 (MQ=255)
                                               gcCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCg                                                                                                   <  3:332527/149‑1 (MQ=255)
                                                 ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                 <  2:380428/149‑1 (MQ=255)
                                                 ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                 >  5:100348/1‑149 (MQ=255)
                                                 ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                 >  7:179926/1‑149 (MQ=255)
                                                  cGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAAc                                                                                                <  7:19891/149‑1 (MQ=255)
                                                   gCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCCACAGGCGCCGCCCCTCTCGAAcc                                                                                               >  3:236092/1‑149 (MQ=255)
                                                    cTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCt                                                                                              <  3:50917/149‑1 (MQ=255)
                                                     tCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCt                                                                                              >  7:338621/1‑148 (MQ=255)
                                                        aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                          <  8:27110/149‑1 (MQ=255)
                                                        aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                          <  1:308062/149‑1 (MQ=255)
                                                        aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                          <  8:179926/149‑1 (MQ=255)
                                                          cccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                         >  5:6501/1‑148 (MQ=255)
                                                           ccGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAAc                                                                                                                    >  6:237198/1‑120 (MQ=255)
                                                            cGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCAt                                                                                      >  2:60487/1‑149 (MQ=255)
                                                               tctcCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTcc                                                                                   <  5:335098/149‑1 (MQ=255)
                                                                   cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                <  6:385145/148‑1 (MQ=255)
                                                                     gccgccCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGc                                                                             <  4:308820/149‑1 (MQ=255)
                                                                        gccCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGCCGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGcac                                                                           >  2:116778/1‑148 (MQ=255)
                                                                              cGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGcccc                                                                    <  8:246812/149‑1 (MQ=255)
                                                                                gCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCa                                                                  >  5:400359/1‑149 (MQ=255)
                                                                                gCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCa                                                                  <  5:321860/149‑1 (MQ=255)
                                                                                gCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCACGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAAGCACATGccccc                                                                   >  7:145209/1‑148 (MQ=255)
                                                                                    gCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAg                                                               >  5:16593/1‑148 (MQ=255)
                                                                                       ggCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCaaa                                                           <  1:60487/149‑1 (MQ=255)
                                                                                       ggCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCaaa                                                           <  3:162431/149‑1 (MQ=255)
                                                                                          gcAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGcc                                                        >  7:320263/1‑149 (MQ=255)
                                                                                            aaCCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCg                                                      >  3:267347/1‑149 (MQ=255)
                                                                                                 gCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCt                                                 <  2:113241/149‑1 (MQ=255)
                                                                                                    gTACCCCGCATGCTCTTCGCTAACCTGCAAGTCTACGCTGCGCGCCGCGAACTGAATGCAGGTAGACGGCGAGGCCAACAGGCGCCGCCCATCTCGAACCTCAACCCAATCCTCATCCACAACCCCCCATAGCAAAGCCCGTACCtctg                                              >  1:393797/1‑149 (MQ=255)
                                                                                                        cccGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGata                                          >  8:268554/1‑149 (MQ=255)
                                                                                                          cGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATc                                        >  8:344336/1‑149 (MQ=255)
                                                                                                                 tctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAAc                                 <  8:320263/149‑1 (MQ=255)
                                                                                                                  ctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAAc                                 >  8:103965/1‑148 (MQ=255)
                                                                                                                        tCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATc                           >  2:329618/1‑148 (MQ=255)
                                                                                                                          aCCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGaa                        >  5:364769/1‑149 (MQ=255)
                                                                                                                          aCCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGaa                        >  5:358805/1‑149 (MQ=255)
                                                                                                                                  gTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCt                >  4:183430/1‑149 (MQ=255)
                                                                                                                                    cTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTgg              >  5:302692/1‑149 (MQ=255)
                                                                                                                                     ttCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTggg             <  6:16593/149‑1 (MQ=255)
                                                                                                                                     ttCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTggg             >  3:9758/1‑149 (MQ=255)
                                                                                                                                      tCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTggg             >  8:261933/1‑148 (MQ=255)
                                                                                                                                         cTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt         <  7:171501/149‑1 (MQ=255)
                                                                                                                                            cgcgcCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg      >  6:361024/1‑149 (MQ=255)
                                                                                                                                             gcgcCGCGAACTGAATGCAGGTCGACGGCGTGGCCCACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCCCATGCCCCCATAGCAAAGCCCGTACCGGGGGATATCCTTCAACGATAGCGCACAGCGCCTGGGCATTgcg      >  5:391389/1‑148 (MQ=255)
                                                                                                                                               gccgcGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCaa   <  6:7238/149‑1 (MQ=255)
                                                                                                                                                ccgcGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAg  >  7:105799/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAG  >  NC_002947/5589922‑5590214

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