Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 421,165 T→C 5.6% H928R (CAC→CGC)  PP_0348 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947421,1650TC 5.6% 259.4 / 4.1 90H928R (CAC→CGC) PP_0348hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/43);  new base C (2/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGT  >  NC_002947/421017‑421309
                                                                                                                                                    |                                                                                                                                                  
gCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGt                                                                                                                                                    >  6:75899/1‑149 (MQ=255)
  aCGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTgg                                                                                                                                                  <  6:285679/149‑1 (MQ=255)
    gCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGc                                                                                                                                               <  5:327323/149‑1 (MQ=255)
    gCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGc                                                                                                                                               >  1:53027/1‑149 (MQ=255)
         ccAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCCCACTGACAGCGGAGCcggcg                                                                                                                                          <  6:97408/149‑1 (MQ=255)
          cagctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCcggcgg                                                                                                                                         <  8:270621/149‑1 (MQ=255)
            gctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGt                                                                                                                                                    >  2:74054/1‑137 (MQ=255)
            gctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCggtg                                                                                                                                       <  1:284496/149‑1 (MQ=255)
              tcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCggtggt                                                                                                                                     >  2:151270/1‑149 (MQ=255)
                agctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTg                                                                                                                                   <  2:273660/149‑1 (MQ=255)
                agctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGACCGCCGGCCCACTGACAGCGGAGCcggcggcctc                                                                                                                                     <  7:245852/147‑5 (MQ=38)
                       ccccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTGGCCAAAACTGGTCGACTCAGGCCGCCGGCTCCGCTGTCTGGGGGCCGGCGGTGGTTGATCTGGc                                                                                                                            <  3:316158/149‑1 (MQ=255)
                        cccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCATTGGGCCGGCGGTGGTTGATCTgg                                                                                                                             >  4:81932/1‑147 (MQ=255)
                        cccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTgg                                                                                                                             >  5:300082/1‑147 (MQ=255)
                         ccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCTGTGGTTTATCTGGCAg                                                                                                                          <  8:99584/149‑1 (MQ=255)
                              cTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCa                                                                                                                     >  4:94258/1‑149 (MQ=255)
                               tATGCCCGTGACCGGGCACGGCTGCTACGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAg                                                                                                                    >  2:282372/1‑149 (MQ=255)
                                 tGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGc                                                                                                                  <  5:75899/148‑1 (MQ=255)
                                 tGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCACCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGc                                                                                                                  <  8:248944/148‑1 (MQ=255)
                                    ccGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGAGTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAgt                                                                                                                         <  4:232566/139‑1 (MQ=255)
                                    ccGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAgt                                                                                                                         >  3:232566/1‑139 (MQ=255)
                                     cGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGt                                                                                                              >  6:195357/1‑148 (MQ=255)
                                     cGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCAGCTCCGCTGTCAGTGGGCCGGCCGTGGTTGATCTGGCAGTGTCAGCCAGt                                                                                                              >  3:76312/1‑148 (MQ=255)
                                         aCCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCActg                                                                                                         <  2:281679/149‑1 (MQ=255)
                                          ccGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCActgc                                                                                                        >  8:118554/1‑149 (MQ=255)
                                           cGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGGCAGTGGGCCGGCGGTGGTTTATCTGGCAGTGTCAGCCAGTCActgct                                                                                                       >  6:115333/1‑149 (MQ=255)
                                               cACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTgcgc                                                                                                   <  7:118554/149‑1 (MQ=255)
                                                aCGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGACAGTCACTGATGCGcc                                                                                                  >  6:256201/1‑149 (MQ=255)
                                                  ggCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcct                                                                                                >  1:221954/1‑149 (MQ=255)
                                                    ctgctgCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgc                                                                                              <  2:221954/149‑1 (MQ=255)
                                                      gctgcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgcct                                                                                            >  6:74644/1‑149 (MQ=255)
                                                       ctgcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgcctg                                                                                           <  2:133131/149‑1 (MQ=255)
                                                         gcgcAACTCCCCCACTGGCAGTGCCTTCAATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGCCCGGCGGTGGTTGATCTCGCAGTCTCAGCCAGTCACTGCTGCCCcctcaatggc                                                                                         >  8:252887/1‑142 (MQ=255)
                                                          cgcAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTgg                                                                                          <  1:282372/147‑1 (MQ=255)
                                                             aaCTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGTTGTTGCCAATACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAc                                                                                     <  5:323803/149‑1 (MQ=255)
                                                               cTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  5:74644/149‑1 (MQ=255)
                                                               cTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  8:107918/149‑1 (MQ=255)
                                                               cTGGCACACTGGCAGTGCCTTGCAGCAGCGGATGGACCGGCTGTTGCCAAACCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGGCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  5:115333/149‑1 (MQ=255)
                                                                      aCTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAg                                                                            >  7:183987/1‑149 (MQ=255)
                                                                       cTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGc                                                                           >  5:137624/1‑149 (MQ=255)
                                                                          gCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGATCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        >  1:231759/1‑149 (MQ=255)
                                                                            aGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTg                                                                      <  1:151270/149‑1 (MQ=255)
                                                                             gTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTgg                                                                     <  2:82697/149‑1 (MQ=255)
                                                                                     cATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTc                                                             <  3:94258/149‑1 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgcctg                                                                                           >  2:57606/1‑117 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgcctg                                                                                           <  1:57606/117‑1 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        >  4:139269/1‑136 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        <  3:139269/136‑1 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGt                                                           >  3:36715/1‑149 (MQ=255)
                                                                                       tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCATATCAGCCACTGGCCGTTTTCGt                                                           >  5:78933/1‑149 (MQ=255)
                                                                                        cAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCATATCAGCCACTGGCCGTTTTCGt                                                           <  6:78933/148‑1 (MQ=255)
                                                                                         aGCTGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTgg                                                         <  2:53027/149‑1 (MQ=255)
                                                                                         aGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCAGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTgg                                                         <  4:76312/149‑1 (MQ=255)
                                                                                             gCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCGCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgca                                                     <  6:300082/149‑1 (MQ=255)
                                                                                               tGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagc                                                   >  3:245101/1‑149 (MQ=255)
                                                                                                 gACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagca                                                  <  5:195357/148‑1 (MQ=255)
                                                                                                  aCCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCAcc                                                >  2:64902/1‑149 (MQ=255)
                                                                                                    cGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTg                                                          <  6:218423/137‑1 (MQ=255)
                                                                                                    cGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTg                                                          >  5:218423/1‑137 (MQ=255)
                                                                                                          ttGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCAc                                          <  1:323482/147‑1 (MQ=255)
                                                                                                            gCCAAAACTGGTCGACTGAGGCCGCCGCCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCACTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                      >  7:37934/1‑149 (MQ=255)
                                                                                                             ccAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                      <  8:284827/148‑1 (MQ=255)
                                                                                                                aaaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcaccac                                  <  3:77696/149‑1 (MQ=255)
                                                                                                                 aaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGACTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCcccccc                                 >  2:212964/1‑149 (MQ=255)
                                                                                                                     ggTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCgg                             >  2:320860/1‑149 (MQ=255)
                                                                                                                      gTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGACGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACAGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCCTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCCGc                            >  7:187752/1‑149 (MQ=255)
                                                                                                                       tCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCg                           >  8:108691/1‑149 (MQ=255)
                                                                                                                         gACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGtt                         >  2:142900/1‑149 (MQ=255)
                                                                                                                         gACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGtt                         >  5:222801/1‑149 (MQ=255)
                                                                                                                          aCTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGGCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATGAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTg                        <  3:135067/149‑1 (MQ=255)
                                                                                                                          aCTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTg                        <  7:252887/149‑1 (MQ=255)
                                                                                                                            tGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGcc                      >  1:322047/1‑148 (MQ=255)
                                                                                                                             gAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     >  6:255312/1‑149 (MQ=255)
                                                                                                                             accaccgccgGCCCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     >  4:130639/5‑149 (MQ=255)
                                                                                                                             accaccgccgGCCCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     <  3:130639/145‑1 (MQ=255)
                                                                                                                              ccaccgccgGCCCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  2:42797/4‑149 (MQ=255)
                                                                                                                              aGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    <  2:303786/149‑1 (MQ=255)
                                                                                                                              aGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  3:305039/1‑149 (MQ=255)
                                                                                                                              aGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  3:310041/1‑149 (MQ=255)
                                                                                                                               ggCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGAGCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  1:84916/1‑149 (MQ=255)
                                                                                                                                 ccgccgGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCaa                 <  5:88287/149‑1 (MQ=255)
                                                                                                                                 ccgcagGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGAACCGCCCACGGCACCACCCCGGCGTTGCTGGCac                 >  6:106351/1‑148 (MQ=255)
                                                                                                                                   gccgGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACg               >  6:91630/1‑149 (MQ=255)
                                                                                                                                     cgGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGc             <  4:245101/149‑1 (MQ=255)
                                                                                                                                      ggCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACgccg            <  2:322047/149‑4 (MQ=255)
                                                                                                                                       gCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACgccgg           <  5:91630/149‑5 (MQ=255)
                                                                                                                                        cTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAt          <  6:137624/149‑1 (MQ=255)
                                                                                                                                           cGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       >  5:149130/1‑149 (MQ=255)
                                                                                                                                           cGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       <  4:36715/149‑1 (MQ=255)
                                                                                                                                            gCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       >  1:253282/1‑148 (MQ=255)
                                                                                                                                              tGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCaa    <  6:222801/149‑1 (MQ=255)
                                                                                                                                                tCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGt  <  1:142900/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGGGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGT  >  NC_002947/421017‑421309

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