Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,885,049 C→A 8.3% Q172K (CAG→AAG) ‡ PP_2540 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,885,0490CA 8.3% 187.6 / 6.4 60Q172K (CAG→AAG) ‡PP_2540short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base A (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GCGGTGGTAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTT  >  NC_002947/2884902‑2885196
                                                                                                                                                   |                                                                                                                                                   
gCGGTGGTAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGgcagca                                                                                                                                                    <  4:295376/149‑1 (MQ=255)
  ggtggtAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACTGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCagcagc                                                                                                                                                  <  6:195703/149‑1 (MQ=255)
      gtAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTc                                                                                                                                               <  6:192366/148‑1 (MQ=255)
        aGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTgg                                                                                                                                            >  7:55879/1‑149 (MQ=255)
         gATGGCATGGCCGAGAGAGGTATTGGCCGGGTGGTGAATATAACCTCAGGCGCGGTGAAGGCACCGATCGACAGGCTGGGATTATCCAAGGGTGCCCGCAGTGGCGGGCCCGGATTCATAGCCGCGCTGGCCCGGCACCGGCGGCTggg                                                                                                                                           >  2:290350/1‑148 (MQ=255)
             gcatggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                       <  4:101860/149‑1 (MQ=255)
               atggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTg                                                                                                                                     >  5:15343/1‑149 (MQ=255)
               atggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTg                                                                                                                                     >  8:2325/1‑149 (MQ=255)
                      aGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgtg                                                                                                                              <  7:2325/149‑1 (MQ=255)
                       gCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGa                                                                                                                             <  5:60092/149‑1 (MQ=255)
                             ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTca                                                                                                                       <  8:54709/149‑1 (MQ=255)
                             ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                       >  4:303043/1‑133 (MQ=255)
                             ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                       <  3:303043/133‑1 (MQ=255)
                                  gCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCaacaac                                                                                                                  >  5:106045/1‑149 (MQ=255)
                                       gtggtgAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGcc                                                                                                                                      >  1:259977/1‑123 (MQ=255)
                                       gtggtgAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGcc                                                                                                                                      <  2:259977/124‑2 (MQ=255)
                                            gAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACACCCTGTTGCCTg                                                                                                        >  4:66178/1‑149 (MQ=255)
                                              atatCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      <  6:15343/149‑1 (MQ=255)
                                                  cACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                               >  2:55657/1‑120 (MQ=255)
                                                  cACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                               <  1:55657/120‑1 (MQ=255)
                                                   aCCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGtt                                                                                                 <  1:78125/149‑1 (MQ=255)
                                                     cTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTAACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                               >  3:171255/1‑117 (MQ=255)
                                                     cTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTAACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                               <  4:171255/117‑1 (MQ=255)
                                                           cgcgGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACgg                                                                                         <  7:20858/149‑1 (MQ=255)
                                                                tGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                              <  7:200201/139‑1 (MQ=255)
                                                                    ggCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGc                                                                                <  6:106045/149‑1 (MQ=255)
                                                                     gCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                              <  7:58308/134‑1 (MQ=255)
                                                                     gCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                              >  8:58308/1‑134 (MQ=255)
                                                                       aCCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCaa                                                                                                                                  <  6:8661/96‑1 (MQ=255)
                                                                       aCCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCaa                                                                                                                                  >  5:8661/1‑96 (MQ=255)
                                                                         cGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCaa                                                                              >  8:292813/1‑146 (MQ=255)
                                                                         cGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAg                                                                           <  7:292813/149‑1 (MQ=255)
                                                                                   gCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      >  4:234078/1‑112 (MQ=255)
                                                                                   gCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      <  3:234078/112‑1 (MQ=255)
                                                                                    cTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGAGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAgg                                                                 <  8:55879/148‑1 (MQ=255)
                                                                                     tGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGCCACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCt                                                               >  5:227249/1‑149 (MQ=255)
                                                                                      gggTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTg                                                              >  3:115659/1‑149 (MQ=255)
                                                                                         tttATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCgg                                                           >  7:229958/1‑149 (MQ=255)
                                                                                            aTCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      >  5:95599/1‑103 (MQ=255)
                                                                                            aTCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      <  6:95599/103‑1 (MQ=255)
                                                                                             tCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTc                                                                                                >  8:230763/1‑108 (MQ=255)
                                                                                             tCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTc                                                                                                <  7:230763/108‑1 (MQ=255)
                                                                                               cAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGggtggt                                                      <  8:229958/148‑1 (MQ=255)
                                                                                                aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCt                                                                                  <  2:306474/119‑1 (MQ=255)
                                                                                                aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGa                                                    >  4:134660/1‑149 (MQ=255)
                                                                                                aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGa                                                    >  4:207422/1‑149 (MQ=255)
                                                                                                 aCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGAt                                                   <  3:134660/149‑1 (MQ=255)
                                                                                                    gTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGAtgtg                                                <  3:153746/149‑1 (MQ=255)
                                                                                                       cccGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTgg                                               >  1:12034/1‑147 (MQ=255)
                                                                                                        ccGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAg                                            >  7:4165/1‑149 (MQ=255)
                                                                                                          ggaGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTGAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTc                                          <  1:290350/147‑1 (MQ=255)
                                                                                                              tGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      <  6:35111/85‑1 (MQ=255)
                                                                                                              tGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                      >  5:35111/1‑85 (MQ=255)
                                                                                                                  cTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGcc                                    >  4:270937/1‑147 (MQ=255)
                                                                                                                        ggATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcg                            >  3:168980/1‑149 (MQ=255)
                                                                                                                        ggATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcg                            >  6:40050/1‑149 (MQ=255)
                                                                                                                          aTTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcggc                          <  3:270937/149‑1 (MQ=255)
                                                                                                                               ttGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGTCGCGCCTGCAACCGACCTTGAAGGCTGCGGGCGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAg                     >  7:205928/1‑149 (MQ=255)
                                                                                                                                 gCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAgg                    <  8:4165/148‑1 (MQ=255)
                                                                                                                                  cccgCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACACAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCACCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGc                   >  6:27849/4‑148 (MQ=255)
                                                                                                                                        tGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCCTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCg            >  4:133416/1‑149 (MQ=255)
                                                                                                                                            ccGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGt        >  3:97539/1‑149 (MQ=255)
                                                                                                                                                 agcagcGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTCGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTAGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGt   >  8:266601/1‑149 (MQ=255)
                                                                                                                                                  gcagcGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGtt  <  3:207422/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGGTGGTAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTT  >  NC_002947/2884902‑2885196

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