Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,320,489 G→C 5.6% A251G (GCC→GGC)  PP_2041 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,320,4890GC 5.6% 233.8 / 4.0 71A251G (GCC→GGC) PP_2041LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATGCG  >  NC_002947/2320341‑2320637
                                                                                                                                                    |                                                                                                                                                     
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTgg                                                                                                                                                       >  6:83110/1‑149 (MQ=255)
    cAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAg                                                                                                                                                    <  5:366695/148‑1 (MQ=255)
     aGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGgc                                                                                                                                                 >  3:228747/1‑149 (MQ=255)
      gcgcTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGgcg                                                                                                                                                >  6:316525/1‑149 (MQ=255)
       cgctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATAAGCCGCGCGCTGGCAGgcgc                                                                                                                                               >  6:331999/1‑149 (MQ=255)
        gctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGcgcg                                                                                                                                              >  4:215405/1‑149 (MQ=255)
         ctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCCTGCCAGCGCGCgg                                                                                                                                             <  5:137682/149‑1 (MQ=255)
           cgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGcc                                                                                                                                           <  7:160208/149‑1 (MQ=255)
            gctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGcc                                                                                                                                           <  5:316525/148‑1 (MQ=255)
             ctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGc                                                                                                                                            >  7:92665/1‑146 (MQ=255)
                cgcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCCTGCCAGCGCGCGGctgata                                                                                                                                       <  4:123724/148‑6 (MQ=255)
                 gcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATcc                                                                                                                                     >  6:352570/1‑149 (MQ=255)
                 gcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATcc                                                                                                                                     <  8:29521/149‑1 (MQ=255)
                  cACCGCATCCACCAGCCGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCca                                                                                                                                    <  6:258459/149‑1 (MQ=255)
                  cACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCca                                                                                                                                    <  8:344232/149‑1 (MQ=255)
                   aCCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGAGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcag                                                                                                                                   >  4:330128/1‑149 (MQ=255)
                     cGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcagca                                                                                                                                 <  1:115598/149‑1 (MQ=255)
                           caccaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAgcggcg                                                                                                                           <  7:270475/149‑1 (MQ=255)
                           caccaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAgcggc                                                                                                                            <  2:362405/148‑1 (MQ=255)
                            cacacccgTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGcggcgg                                                                                                                          <  5:397723/142‑1 (MQ=255)
                                 cagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGt                                                                                                                     >  3:282519/1‑149 (MQ=255)
                                 cagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGt                                                                                                                     >  3:279103/1‑149 (MQ=255)
                                       gcgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGa                                                                                                               <  3:148158/149‑1 (MQ=255)
                                        cgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCATGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAt                                                                                                              <  2:252217/149‑1 (MQ=255)
                                            ccAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtggtg                                                                                                          >  6:13152/1‑149 (MQ=255)
                                                gTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgccg                                                                                                      >  3:115540/1‑149 (MQ=255)
                                                gTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgccg                                                                                                      >  1:280941/1‑149 (MQ=255)
                                                   gcgcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgcca                                                                                                   <  7:268561/149‑1 (MQ=255)
                                                    cgcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGGTGGTGGCCgccag                                                                                                  <  4:228747/149‑1 (MQ=255)
                                                     gcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccagg                                                                                                 >  2:209507/1‑149 (MQ=255)
                                                      cATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccaggc                                                                                                >  2:338993/1‑149 (MQ=255)
                                                              cggcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACgg                                                                                        <  4:279103/149‑1 (MQ=255)
                                                              cggcggcgAACGACAGCGGCTTGGGCTGCGAGGGGCCCAAGGCGATGGTGCGGGCGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACgg                                                                                        <  4:115540/149‑1 (MQ=255)
                                                                gcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCTCCAGGCCAGCACGGAc                                                                                      <  2:256382/149‑1 (MQ=255)
                                                                gcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGAc                                                                                      <  2:250899/149‑1 (MQ=255)
                                                                 cggcgAACGACAGCTGATTGTGCTGCGTGTGGGCCAAGGGTATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCTGCCAGATCCAGCTGCTGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGAct                                                                                     <  5:69888/149‑2 (MQ=255)
                                                                         gACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGt                                                                             <  8:92665/149‑1 (MQ=255)
                                                                          aCAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGt                                                                             >  6:107807/1‑148 (MQ=255)
                                                                                  ttGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCggg                                                                    >  3:102569/1‑149 (MQ=255)
                                                                                     tgctgaGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGc                                                                 <  1:264668/149‑1 (MQ=255)
                                                                                        tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGa                                                              >  1:241648/1‑149 (MQ=255)
                                                                                          cGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGa                                                              <  4:282519/147‑1 (MQ=255)
                                                                                          cGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGAct                                                            >  4:287852/1‑149 (MQ=255)
                                                                                             gTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgct                                                         <  4:102569/149‑1 (MQ=255)
                                                                                             gTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgct                                                         <  5:83110/149‑1 (MQ=255)
                                                                                                 gCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTc                                                     <  5:13152/149‑1 (MQ=255)
                                                                                                 gCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTc                                                     >  4:7280/1‑149 (MQ=255)
                                                                                                  ccAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCg                                                    >  3:164511/1‑149 (MQ=255)
                                                                                                   cAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGt                                                   <  1:338993/149‑1 (MQ=255)
                                                                                                    aaGGCGATGGTGCGGGGGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGGGATACCGGGGGCCGACTGCTGTTCGTa                                                  <  2:266466/149‑1 (MQ=255)
                                                                                                     aGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAg                                                 >  5:235033/1‑149 (MQ=255)
                                                                                                        cGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTg                                              <  3:7280/149‑1 (MQ=255)
                                                                                                              tGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTc                                         >  3:290662/1‑148 (MQ=255)
                                                                                                               gCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                       <  7:267810/149‑1 (MQ=255)
                                                                                                                 gggTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgt                                     >  4:364734/1‑149 (MQ=255)
                                                                                                                  ggTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgtg                                    >  5:81892/1‑149 (MQ=255)
                                                                                                                      aTGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTcac                                <  5:352570/149‑1 (MQ=255)
                                                                                                                          gcgGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCg                            <  3:287852/149‑1 (MQ=255)
                                                                                                                             gTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                       <  4:390617/135‑1 (MQ=255)
                                                                                                                             gTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                       >  3:390617/1‑135 (MQ=255)
                                                                                                                                ccACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTc                       >  2:113171/1‑148 (MQ=255)
                                                                                                                                 cACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCa                     <  6:318745/149‑1 (MQ=255)
                                                                                                                                  aCTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCAt                    >  1:362374/1‑149 (MQ=255)
                                                                                                                                    tATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATgg                  <  3:215405/149‑1 (MQ=255)
                                                                                                                                       tggCCGCGCCTGCCAGCGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGCCTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAg               >  3:91951/3‑149 (MQ=255)
                                                                                                                                       tggCCGCGCCTGCCAGCGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAg               >  4:10888/3‑149 (MQ=255)
                                                                                                                                           cgcgcgCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTg           >  7:260510/1‑149 (MQ=255)
                                                                                                                                            gcgcgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGc          >  2:316595/1‑149 (MQ=255)
                                                                                                                                            gcgcgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGc          <  3:330128/149‑1 (MQ=255)
                                                                                                                                             cgcgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCt         <  2:280941/149‑1 (MQ=255)
                                                                                                                                             cgcgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCt         <  3:364734/149‑1 (MQ=255)
                                                                                                                                                 cTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGAt     <  1:316595/149‑1 (MQ=255)
                                                                                                                                                    gCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACATCGGTTTCCATGGCAGCCTGCTCGATgcg  >  1:160632/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATGCG  >  NC_002947/2320341‑2320637

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: