Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,083 G→T 6.9% 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 6.9% 232.8 / 6.2 73T214N (ACC→AAC) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (41/27);  new base T (2/3);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

GCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATCCA  >  NC_002947/5589935‑5590231
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                           >  6:361737/1‑141 (MQ=255)
                   aCCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTg                                                                                                                                           <  5:361737/141‑1 (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                                                                                                                  >  5:100348/1‑149 (MQ=255)
                                    ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                                  >  7:179926/1‑149 (MQ=255)
                                    ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                                  <  2:380428/149‑1 (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                                                                                                           <  8:179926/149‑1 (MQ=255)
                                           aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                                           <  1:308062/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:321860/149‑1 (MQ=255)
                                                                   gCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCa                                                                                   >  5:400359/1‑149 (MQ=255)
                                                                   gCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCACGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAAGCACATGccccc                                                                                    >  7:145209/1‑148 (MQ=255)
                                                                       gCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAg                                                                                >  5:16593/1‑148 (MQ=255)
                                                                          ggCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCaaa                                                                            <  3:162431/149‑1 (MQ=255)
                                                                          ggCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCaaa                                                                            <  1:60487/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)
                                                                                                                                           cTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCg            <  6:302692/149‑1 (MQ=255)
                                                                                                                                              aaTGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGgcg         >  3:126178/1‑149 (MQ=255)
                                                                                                                                              aaTGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCAGTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGgcg         >  2:253229/1‑149 (MQ=255)
                                                                                                                                                   aGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATc    >  4:132416/1‑149 (MQ=255)
                                                                                                                                                     gTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATc    >  1:300913/1‑147 (MQ=255)
                                                                                                                                                     gTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATCCa  <  7:103965/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGATCGGCGCGATCCA  >  NC_002947/5589935‑5590231

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: