Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,072 A→C 6.7% F218V (TTC→GTC)  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,0720AC 6.7% 229.4 / 6.3 75F218V (TTC→GTC) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (43/27);  new base C (2/3);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAG  >  NC_002947/5589925‑5590214
                                                                                                                                                    |                                                                                                                                              
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                                                                                                 >  7:179926/1‑149 (MQ=255)
                                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                 <  2:380428/149‑1 (MQ=255)
                                              ccGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                 >  5:100348/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:179926/149‑1 (MQ=255)
                                                     aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                          <  1:308062/149‑1 (MQ=255)
                                                     aCCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATc                                                                                          <  8:27110/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                                                           <  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)
                                                                                                                                                    |                                                                                                                                              
AATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGCGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAG  >  NC_002947/5589925‑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: