Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,590,072 A→C 6.6% 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.6% 230.3 / 6.8 76F218V (TTC→GTC) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (30/41);  new base C (3/2);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAT  >  NC_002947/5589924‑5590219
                                                                                                                                                      |                                                                                                                                                   
cAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGaa                                                                                                                                                     >  3:138820/1‑149 (MQ=255)
 aaTACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGaa                                                                                                                                                     <  2:53366/148‑1 (MQ=255)
 aaTACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAAc                                                                                                                                                    <  4:347359/149‑1 (MQ=255)
  aTACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACt                                                                                                                                                   <  7:195194/149‑1 (MQ=255)
     cccTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGaa                                                                                                                                                <  3:398017/149‑1 (MQ=255)
          ggCTTTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGc                                                                                                                                             <  4:109103/147‑1 (MQ=255)
          ggCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAg                                                                                                                                           >  5:234869/1‑149 (MQ=255)
                 gAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGGCGACg                                                                                                                                    <  2:151990/149‑1 (MQ=255)
                     tCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGt                                                                                                                                <  6:313713/149‑1 (MQ=255)
                     gcGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGTCGCAACCTGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCg‑cggcgc                                                                                                                                <  2:29607/148‑2 (MQ=255)
                      cGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTg                                                                                                                               >  1:242068/1‑149 (MQ=255)
                       gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTgg                                                                                                                              >  5:410504/1‑149 (MQ=255)
                       gAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTg                                                                                                                               <  7:245254/148‑1 (MQ=255)
                         gCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGTGCCGTCGACCTGCATTCAGTTCg‑cggcgcgcag                                                                                                                            <  6:129832/149‑6 (MQ=255)
                         gCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGc                                                                                                                             >  4:51591/1‑148 (MQ=255)
                            gCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCaa                                                                                                                          >  6:309458/1‑148 (MQ=255)
                             cACCTTGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAAc                                                                                                                         <  5:218606/148‑1 (MQ=255)
                             cACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCaa                                                                                                                          <  8:20433/147‑1 (MQ=255)
                               ccTCGAAATCCAGCGCCCGCTCCACCCCGGTCTCTCAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAg                                                                                                                       >  1:59122/1‑149 (MQ=255)
                                cTCGAAATCCAGTGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                     >  1:198601/1‑149 (MQ=255)
                                       tCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGccc                                                                                                              >  2:316052/1‑149 (MQ=255)
                                         cAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCct                                                                                                            >  8:133326/1‑149 (MQ=255)
                                            cgcCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCtctc                                                                                                         <  5:248281/149‑1 (MQ=255)
                                               ccGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                     >  5:396312/1‑134 (MQ=255)
                                               ccGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGgc                                                                                                                     <  6:396312/134‑1 (MQ=255)
                                               ccGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                      >  5:321665/1‑149 (MQ=255)
                                               ccGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGAGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGaa                                                                                                      >  6:158854/1‑149 (MQ=255)
                                                cGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCGTCGACCTGCATTCAGTTCgcgg                                                                                                                                    >  8:21670/1‑116 (MQ=255)
                                                cGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CTAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCCGCTCTCGAAc                                                                                                     <  2:242068/149‑1 (MQ=255)
                                                cGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAAc                                                                                                     <  5:184572/149‑1 (MQ=255)
                                                 gCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAAcc                                                                                                    >  7:304210/1‑149 (MQ=255)
                                                     cACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCa                                                                                                 <  7:74981/148‑1 (MQ=255)
                                                       ccccGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGCACCtcctcc                                                                                              >  7:315233/1‑149 (MQ=255)
                                                       ccccGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                              <  8:343285/149‑1 (MQ=255)
                                                       ccccGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                              >  3:205784/1‑149 (MQ=255)
                                                       ccccGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCACCCCTCTCGAACCTCATcc                                                                                              <  3:51591/149‑1 (MQ=255)
                                                        cccGGTCTC‑CAGCCGCCCATCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATcc                                                                                              >  3:124757/1‑148 (MQ=255)
                                                         ccGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCa                                                                                            >  8:86504/1‑149 (MQ=255)
                                                            gTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTc                                                                                         <  8:304210/149‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCTCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                    <  2:382061/149‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCt                                                                                       <  5:309458/146‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                     <  6:234869/148‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                     <  6:149290/148‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                    <  6:264092/149‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGAt                                                                                    <  6:321665/149‑1 (MQ=255)
                                                                  cAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGAGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGa                                                                                     <  5:158854/148‑1 (MQ=255)
                                                                        cccGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CTAACTGAATGCAGGTCGACGGCGTGGCTAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                              <  1:316052/149‑1 (MQ=255)
                                                                        cccGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACAt                                                                              >  3:208436/1‑149 (MQ=255)
                                                                              gtcATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGccca                                                                         <  8:220704/146‑2 (MQ=255)
                                                                                cATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCAt                                                                      <  3:161445/149‑1 (MQ=255)
                                                                                             ccAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTa                                                         >  5:222180/1‑149 (MQ=255)
                                                                                                 cGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCtg                                                     <  2:59122/149‑1 (MQ=255)
                                                                                                    tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                  <  4:208436/149‑1 (MQ=255)
                                                                                                    tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                  <  4:205784/149‑1 (MQ=255)
                                                                                                    tACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTgg                                                  <  4:124757/149‑1 (MQ=255)
                                                                                                      ccccGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGa                                                 >  8:20264/1‑148 (MQ=255)
                                                                                                        ccGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGatat                                              <  8:118201/149‑1 (MQ=255)
                                                                                                              gCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTc                                         >  1:235246/1‑148 (MQ=255)
                                                                                                               ctctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCaa                                       <  2:198601/149‑1 (MQ=255)
                                                                                                                 ctTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAAc                                      >  1:344931/1‑148 (MQ=255)
                                                                                                                       tCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATc                                >  2:129716/1‑148 (MQ=255)
                                                                                                                       tCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATCCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATc                                >  8:46349/1‑148 (MQ=255)
                                                                                                                       tCACCTGGAAGTCTTCGCTGCCCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGCCGCCGCCCCTCTCGAACCACATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGataac                                >  4:171527/1‑146 (MQ=255)
                                                                                                                         aCCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGaa                             >  3:309326/1‑149 (MQ=255)
                                                                                                                           cTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACa                           <  6:222180/149‑1 (MQ=255)
                                                                                                                            tGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAg                          >  4:159552/1‑149 (MQ=255)
                                                                                                                             ggAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCAACTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCCTAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAgc                         >  7:200640/1‑149 (MQ=255)
                                                                                                                                      cGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTgg                   <  7:20264/146‑1 (MQ=255)
                                                                                                                                        gtGCGCGGCG‑CGTACTGAATGCAGGTCGACGGCGTGGGCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCTAAGCCCGTACCTGTGGATATCGTTCAACGATATCGAACAGCGACTGGGCAtt              <  3:171527/148‑1 (MQ=255)
                                                                                                                                        cTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt              <  7:86504/149‑1 (MQ=255)
                                                                                                                                        cTGCGCGCCG‑CGAACTGAATGCAGGTCGAAGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt              <  7:46349/149‑1 (MQ=255)
                                                                                                                                        atGCGCGCCG‑CTAACTGAATGCAGTGCGACTGCGTGGCCACCAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCAtt              <  8:315233/148‑1 (MQ=255)
                                                                                                                                         tGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTg             >  4:194480/1‑149 (MQ=255)
                                                                                                                                          gcgcgcCG‑CGAACTGAATGCAGGGCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCTCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgc            <  4:101477/149‑1 (MQ=255)
                                                                                                                                           cgcgcCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg           >  2:241172/1‑149 (MQ=255)
                                                                                                                                            acgccgtcgACCTGCATTCAGTTCG‑CGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg           >  3:290730/8‑148 (MQ=255)
                                                                                                                                            acgccgtcgACCTGCATTCAGTTCG‑CGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTgcg           >  4:131139/8‑148 (MQ=255)
                                                                                                                                                cg‑cgAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGc      <  1:129716/149‑1 (MQ=255)
                                                                                                                                                     aaCTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAt  <  6:89947/149‑1 (MQ=255)
                                                                                                                                                     aaCTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAt  >  5:89947/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CAATACCCTGGGCTGTCGAAGTCGAGCTGCACCTCGAAATCCAGCGCCCGCTCCACCCCGGTCTC‑CAGCCGCCCGTCGGCATGCAGGCGCAACCAGCGGTACCCCGGATGCTCTTCGCTCACCTGGAAGTCTTCGCTGCGCGCCG‑CGAACTGAATGCAGGTCGACGGCGTGGCCAACAGGCGCCGCCCCTCTCGAACCTCATCCCATTCCTGATGCACATGCCCCCATAGCAAAGCCCGTACCTGTGGATATCCTTCAACGATATCGAACAGCGCCTGGGCATTGCGCAAGCCGAT  >  NC_002947/5589924‑5590219

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: ‑