Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,595,528 C→G 7.1% V66L (GTG→CTG)  pde ← 3',5'‑cyclic‑nucleotide phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,595,5280CG 7.1% 986.1 / 48.3 295V66L (GTG→CTG) pde3',5'‑cyclic‑nucleotide phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (124/150);  new base G (12/9);  total (136/159)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 308 total aligned reads displayed.

CACCCGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATCGG  >  NC_002947/5595379‑5595677
                                                                                                                                                     |                                                                                                                                                     
caccCGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCAc                                                                                                                                                       >  1:22371/1‑150 (MQ=255)
      ccACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGc                                                                                                                                                 <  1:6322757/150‑1 (MQ=255)
      ccACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGc                                                                                                                                                 >  1:716732/1‑150 (MQ=255)
            gcCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCAGCGTGGCATCCta                                                                                                                                           <  1:1335249/150‑2 (MQ=255)
             ccGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTgg                                                                                                                                          <  2:3730311/150‑1 (MQ=255)
               gATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCAGCGTGGCATCCtatga                                                                                                                                        <  1:6790875/150‑5 (MQ=255)
                 tatCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACa                                                                                                                                      <  1:2507202/150‑1 (MQ=255)
                   tCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAgg                                                                                                                                    >  1:6769067/1‑150 (MQ=255)
                    cAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGt                                                                                                                                   >  1:4276490/1‑150 (MQ=255)
                     aGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTc                                                                                                                                  <  2:6083246/150‑1 (MQ=255)
                       tCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGc                                                                                                                                >  1:16262/1‑150 (MQ=255)
                       tCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGc                                                                                                                                >  1:2169867/1‑150 (MQ=255)
                       tCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGc                                                                                                                                >  2:5068124/1‑150 (MQ=255)
                         aCCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGccc                                                                                                                              >  1:4413220/1‑150 (MQ=255)
                          ccGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCg                                                                                                                             >  1:3389550/1‑150 (MQ=255)
                           cGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGt                                                                                                                            <  1:4038408/150‑1 (MQ=255)
                            gCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTg                                                                                                                           <  2:5327259/150‑1 (MQ=255)
                            gCCTGAACTGGTGCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCAGGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTg                                                                                                                           >  1:4804150/1‑150 (MQ=255)
                             ccTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGc                                                                                                                          >  1:5938589/1‑150 (MQ=255)
                             ccTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGc                                                                                                                          <  2:1235068/150‑1 (MQ=255)
                              cTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTg                                                                                                                                           <  1:288795/132‑1 (MQ=255)
                              cTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTg                                                                                                                                           >  2:288795/1‑132 (MQ=255)
                                gCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGcaca                                                                                                                       >  2:4552121/1‑150 (MQ=255)
                                 cACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcag                                                                                                                      >  1:2537728/1‑150 (MQ=255)
                                  aCCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagc                                                                                                                     >  1:5705210/1‑150 (MQ=255)
                                  aCCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagc                                                                                                                     <  1:6864019/150‑1 (MQ=255)
                                  aCCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagc                                                                                                                     <  1:5124317/150‑1 (MQ=255)
                                  aCCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagc                                                                                                                     >  2:3622168/1‑150 (MQ=255)
                                   ccAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagca                                                                                                                    <  1:5865068/150‑1 (MQ=255)
                                   ccAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagca                                                                                                                    <  2:2017557/150‑1 (MQ=255)
                                    cAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCAcagcag                                                                                                                   >  2:5670239/1‑150 (MQ=255)
                                     aGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACagcagc                                                                                                                  <  2:554613/150‑1 (MQ=255)
                                      gTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTg                                                                                                                                                   >  1:1772539/1‑116 (MQ=255)
                                      gTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAgcagca                                                                                                                 >  1:2313016/1‑150 (MQ=255)
                                       ttCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCaa                                                                                                                <  1:2843143/150‑1 (MQ=255)
                                        tCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAg                                                                                                               <  2:3071679/150‑1 (MQ=255)
                                          aGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTc                                                                                                             <  2:3389550/150‑1 (MQ=255)
                                            gcgcCACCTGCACCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGa                                                                                                           >  1:2062700/1‑150 (MQ=255)
                                               ccACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACaa                                                                                                        <  1:2670069/150‑1 (MQ=255)
                                                cACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAg                                                                                                       <  2:6516708/150‑1 (MQ=255)
                                                   cTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCgg                                                                                                    >  1:5475540/1‑150 (MQ=255)
                                                    tGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGgc                                                                                                   >  1:2096376/1‑150 (MQ=255)
                                                      cgcCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGgctg                                                                                                 <  1:822468/150‑1 (MQ=255)
                                                         cATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTc                                                                                              >  1:4527703/1‑150 (MQ=255)
                                                         cATCACATGAGCCTCCTCAAGGATCCACCCCAGCAACAGAGTCCCTACGGAAAACGGCGCGGTCAGCTCAAGAAACCCTTCATAGGATGACACGCTGGCATCCTGTGACAGGTCGCCCGTGCACAGCAGCACGTCGACAAGCGGCTGCTc                                                                                              >  2:5425087/1‑150 (MQ=255)
                                                             aCCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCgcgc                                                                                          >  2:860951/1‑150 (MQ=255)
                                                             aCCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCgcgc                                                                                          <  2:6648654/150‑1 (MQ=255)
                                                              ccTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGcgcc                                                                                         <  1:3314372/150‑1 (MQ=255)
                                                                    ccTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCt                                                                                   >  2:743015/1‑150 (MQ=255)
                                                                      tCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTgg                                                                                 >  1:5658344/1‑150 (MQ=255)
                                                                      tCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTgg                                                                                 >  1:3662947/1‑150 (MQ=255)
                                                                      tCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTgg                                                                                 <  1:2783333/150‑1 (MQ=255)
                                                                        gTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGc                                                                               <  1:6035857/150‑1 (MQ=255)
                                                                         tCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGcc                                                                              >  2:4121697/1‑150 (MQ=255)
                                                                           aTGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCaa                                                                            <  1:332098/150‑1 (MQ=255)
                                                                           aTGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCaa                                                                            <  1:5314482/150‑1 (MQ=255)
                                                                                 gCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTg                                                                      <  1:614328/150‑1 (MQ=255)
                                                                                 gCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTg                                                                      <  1:308150/150‑1 (MQ=255)
                                                                                  cccGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTgg                                                                     <  2:6008113/150‑1 (MQ=255)
                                                                                   ccGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCACGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGgc                                                                    >  2:4985930/1‑150 (MQ=255)
                                                                                    cGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGgcg                                                                   >  1:3494194/1‑150 (MQ=255)
                                                                                         gccagcGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGc                                                              <  2:2270941/150‑1 (MQ=255)
                                                                                             gcGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGtc                                                          >  2:4764389/1‑150 (MQ=255)
                                                                                                 gTGGGTACGGCAAACGGTGCGGTCAGCTCCCGACACGCTTCATAGGAGGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCgg                                                      <  1:1267558/150‑1 (MQ=255)
                                                                                                    ggTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTc                                                   >  1:383508/1‑150 (MQ=255)
                                                                                                      tACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCtt                                                 <  2:6833752/150‑1 (MQ=255)
                                                                                                      tACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCtt                                                 <  1:4584442/150‑1 (MQ=255)
                                                                                                        cGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGa                                               <  1:4848690/150‑1 (MQ=255)
                                                                                                         ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAg                                              >  1:1448715/1‑150 (MQ=255)
                                                                                                         ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAg                                              <  2:6255787/150‑1 (MQ=255)
                                                                                                         ggCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAg                                              <  1:981426/150‑1 (MQ=255)
                                                                                                          gCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAgg                                             <  2:76910/150‑1 (MQ=255)
                                                                                                           cAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTgg                                                                                                                                                  >  2:2143255/1‑48 (MQ=255)
                                                                                                           cAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGc                                            <  2:2446711/150‑1 (MQ=255)
                                                                                                            aaaCGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTTCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGcc                                           <  2:4384023/150‑1 (MQ=255)
                                                                                                            aaaCGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGcc                                           >  2:1588197/1‑150 (MQ=255)
                                                                                                              aCGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCa                                         >  1:101680/1‑150 (MQ=255)
                                                                                                               cGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAg                                        >  2:1295669/1‑150 (MQ=255)
                                                                                                                 gTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCa                                      <  1:84569/150‑1 (MQ=255)
                                                                                                                    cGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCAtgc                                   >  2:6717794/1‑150 (MQ=255)
                                                                                                                     ggTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCAtgct                                  <  2:1437280/150‑1 (MQ=255)
                                                                                                                            tCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGcc                           >  1:4111121/1‑150 (MQ=255)
                                                                                                                            tCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGcc                           >  1:1874999/1‑150 (MQ=255)
                                                                                                                             cccGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCg                          >  2:1691681/1‑150 (MQ=255)
                                                                                                                                   aCGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCgg                    >  1:700084/1‑150 (MQ=255)
                                                                                                                                    cGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGc                   <  2:1933823/150‑1 (MQ=255)
                                                                                                                                    cGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGc                   <  2:6136263/150‑1 (MQ=255)
                                                                                                                                     ctgtcccaGGATGCCAGCGTGGCATCCTAGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCg                  >  2:5584034/8‑150 (MQ=255)
                                                                                                                                       gtcccaGGATGCCAGCGTGGCATCCTAGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGaa                >  1:5067255/6‑150 (MQ=255)
                                                                                                                                          ccaGGATGCCAGCGTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAg             >  1:4936285/3‑150 (MQ=255)
                                                                                                                                           caGGATGCCAGCGTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAgg            >  2:487390/2‑150 (MQ=255)
                                                                                                                                            aGGATGCCAGCGTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGt           >  2:6744677/1‑150 (MQ=255)
                                                                                                                                            aGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGt           <  2:3922315/150‑1 (MQ=255)
                                                                                                                                             ggATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTg          <  1:2196647/150‑1 (MQ=255)
                                                                                                                                             ggATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTg          >  1:290730/1‑150 (MQ=255)
                                                                                                                                              gATGCCAGCGTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTgg         <  2:4936285/150‑1 (MQ=255)
                                                                                                                                              gATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTgg         <  1:5248152/150‑1 (MQ=255)
                                                                                                                                              gATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTgg         <  1:6630529/150‑1 (MQ=255)
                                                                                                                                                    aCGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATCg   <  2:5175107/150‑1 (MQ=255)
                                                                                                                                                    aCGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATCg   >  1:518169/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCCGCCACGCGCCGATATCAGTCACCGCCTGCACCAGTTCAGGCGCCACCTGCGCCATCACCTGAGCCTCGTCATGGTTGCCCGGCAGCCAGCGAGTGGGTACGGCAAACGGTGCGGTCAGCTCCCGAAACGCTTCATAGGATGCCACGCTGGCATCCTGGGACAGGTCGCCCGTGCACAGCAGCAAGTCGACAAGCGGCTGCTCGCGCCGCACCTGGGCCAACACGTGGCGCAGGCTGTCTCGGGTCTTGAGGCCCAGCATGCTGCCAGCCGGGTCGGCGAACAGGTGGGCATCGG  >  NC_002947/5595379‑5595677

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