Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,996 C→T 6.6% R868H (CGC→CAC)  PP_1983 ← PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,246,9960CT 6.6% 222.0 / 4.2 75R868H (CGC→CAC) PP_1983PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/3);  total (35/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 = 9.92e-01

GGCGCAGGCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTGGGTGGTGTTGCGA  >  NC_002947/2246848‑2247139
                                                                                                                                                    |                                                                                                                                               
ggCGCAGGCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGc                                                                                                                                                 <  8:220629/149‑1 (MQ=255)
 gcgcAGGCGTCTCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCg                                                                                                                                                <  8:13714/149‑1 (MQ=255)
   gcAGGCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGt                                                                                                                                              >  4:197243/1‑149 (MQ=255)
      ggCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGt                                                                                                                                           >  6:189015/1‑149 (MQ=255)
      ggCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGt                                                                                                                                           <  2:123848/149‑1 (MQ=255)
         gTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGCcg                                                                                                                                         >  1:346027/1‑148 (MQ=255)
         gTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGCcg                                                                                                                                         >  1:347209/1‑148 (MQ=255)
              aGGCTTCGTCGAGCCTCCCGCTGCCATCTTCACGGATCAGGCCACTTTCTTCCCGTACCTGGATGTACTGCGAGGGCGGTTGCTGACCCAGTTGCGTGTGGTGCCAGCGCATCAGTACTTCGGCACCGATGATGTGGTTGTCGCGGGCa                                                                                                                                   <  8:233135/149‑1 (MQ=255)
               ggCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGCTCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCAt                                                                                                                                  <  5:275401/149‑1 (MQ=255)
                     gtcgAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCt                                                                                                                            >  3:126596/1‑149 (MQ=255)
                       cgAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGc                                                                                                                          >  6:376600/1‑149 (MQ=255)
                          gaATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGc                                                                                                                       >  1:279360/1‑149 (MQ=255)
                          gaATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGc                                                                                                                       <  4:126596/149‑1 (MQ=255)
                           aaTCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCAc                                                                                                                              >  2:232554/1‑141 (MQ=255)
                           aaTCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCAc                                                                                                                              <  1:232554/141‑1 (MQ=255)
                               cAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGaa                                                                                                                  <  2:279360/149‑1 (MQ=255)
                                      ctaCTTCAAGGAGCAGGCCACTTTCTTCGAGTTCGTGGATGAACTGCGAGGGCGGTTTCTGAGCCAGTTGCGGGGTGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCTGGCATCCACCTGCGGCTGGAagtgcagt                                                                                                           <  2:127184/147‑1 (MQ=255)
                                      cAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAagtgcagt                                                                                                           <  6:112690/149‑1 (MQ=255)
                                       aaCTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAagtgcagtg                                                                                                          >  6:39419/1‑149 (MQ=255)
                                           tCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAg                                                                                                      >  7:265190/1‑149 (MQ=255)
                                             aaGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCt                                                                                                    >  5:19845/1‑149 (MQ=255)
                                             aaGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCt                                                                                                    >  3:229341/1‑149 (MQ=255)
                                             aaGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCt                                                                                                    >  7:49268/1‑149 (MQ=255)
                                             aaGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCt                                                                                                    >  2:99729/1‑149 (MQ=255)
                                              aGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTc                                                                                                   <  3:197243/149‑1 (MQ=255)
                                                  tCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCa                                                                                               <  5:189015/149‑1 (MQ=255)
                                                     ggCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgg                                                                                            <  6:19845/149‑1 (MQ=255)
                                                     ggCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgg                                                                                            <  6:241054/149‑1 (MQ=255)
                                                      gCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGTGGTTGTCGCGGGCATCCACCTGCGGCTGGaa                                                                                                                  >  5:319618/1‑126 (MQ=255)
                                                      gCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGTGGTTGTCGCGGGCATCCACCTGCGGCTGGaa                                                                                                                  <  6:319618/126‑1 (MQ=255)
                                                      gCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACggg                                                                                           >  2:114523/1‑149 (MQ=255)
                                                       ccACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGATCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggc                                                                                          >  1:262099/1‑149 (MQ=255)
                                                       ccACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggc                                                                                          >  3:338608/1‑149 (MQ=255)
                                                       ccACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggc                                                                                          >  7:342846/1‑149 (MQ=255)
                                                          cTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggccag                                                                                       <  1:99729/149‑1 (MQ=255)
                                                           tttCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggc                                                                                          >  1:226734/1‑145 (MQ=255)
                                                           tttCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggc                                                                                          <  2:226734/145‑1 (MQ=255)
                                                             tcttcCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggccaggg                                                                                     <  6:84741/148‑1 (MQ=255)
                                                             tcttcCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggccagg                                                                                      <  5:376600/147‑1 (MQ=255)
                                                               ttcCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGTCACgggccagggc                                                                                    >  1:346106/1‑147 (MQ=255)
                                                                 ccAGTACCTGGATGAACTGCGAGGGCGGTTGCTGAGCGTGTTGCGGGTGGTGGCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACgggccagggcca                                                                                  <  3:180756/147‑1 (MQ=255)
                                                                 ccAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGggccagggccag                                                                                 <  4:145977/148‑1 (MQ=255)
                                                                    gTACCTGGCTGACCGGCGAGGGCGGTTGCTGACCCAGTGGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCGGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGgcgc                                                                             <  8:342846/149‑1 (MQ=255)
                                                                    gTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCGGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGgcgc                                                                             <  4:338608/149‑1 (MQ=255)
                                                                               aaCTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGtttt                                                                   <  8:279546/148‑1 (MQ=255)
                                                                               aaCTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTc                                                                  <  4:229341/149‑1 (MQ=255)
                                                                                aCTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTcc                                                                 <  2:262099/149‑1 (MQ=255)
                                                                                aCTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGGGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTcc                                                                 <  8:342043/149‑1 (MQ=255)
                                                                                  tGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCAt                                                               >  3:297326/1‑149 (MQ=255)
                                                                                    cGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCAt                                                               <  5:364685/147‑1 (MQ=255)
                                                                                          cGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAg                                                         <  6:33967/147‑1 (MQ=255)
                                                                                          cGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAg                                                         >  5:33967/1‑147 (MQ=255)
                                                                                                cTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGTGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTcc                                                                 <  8:174356/133‑1 (MQ=255)
                                                                                                cTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGTGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTcc                                                                 >  7:174356/1‑133 (MQ=255)
                                                                                                   aCCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTg                                               <  8:265190/148‑1 (MQ=255)
                                                                                                     ccAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATgcgc                                                           <  5:345291/134‑1 (MQ=255)
                                                                                                     ccAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATgcgc                                                           >  6:345291/1‑134 (MQ=255)
                                                                                                      cAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTggcgg                                           >  3:61773/1‑149 (MQ=255)
                                                                                                             gggTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTg                                    <  8:49268/149‑1 (MQ=255)
                                                                                                                   tgCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACGTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCAtggt                              <  4:18216/149‑1 (MQ=255)
                                                                                                                   tgCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCAtggt                              <  7:299580/149‑1 (MQ=255)
                                                                                                                    ccacccgcAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCAtggtg                             <  5:343998/144‑1 (MQ=255)
                                                                                                                     ccAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAagtgcagtgc                                                                                                         >  4:193/1‑72 (MQ=255)
                                                                                                                     ccAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAagtgcagtgc                                                                                                         <  3:193/72‑1 (MQ=255)
                                                                                                                        gcgcAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGgtgt                         >  5:140111/1‑149 (MQ=255)
                                                                                                                         cgcAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGTCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGtgtg                        <  2:346106/149‑1 (MQ=255)
                                                                                                                          gcAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGtgtg                        >  5:36351/1‑148 (MQ=255)
                                                                                                                             aCAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTg                                               <  7:187197/122‑1 (MQ=255)
                                                                                                                             aCAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTg                                               >  8:187197/1‑122 (MQ=255)
                                                                                                                                gTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGc                 >  7:215491/1‑149 (MQ=255)
                                                                                                                                   cTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGaa                      <  6:353870/141‑1 (MQ=255)
                                                                                                                                   cTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGaa                      >  5:353870/1‑141 (MQ=255)
                                                                                                                                   cTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTgg              >  2:92417/1‑149 (MQ=255)
                                                                                                                                    ttCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTggg             >  6:168813/1‑149 (MQ=255)
                                                                                                                                            cgacgaTGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTGGGTGGTGTTg     >  2:64274/1‑149 (MQ=255)
                                                                                                                                               cgaTGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTGGGTGGTGTTGCGa  <  1:114523/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGCGCAGGCGTCGCAGGCTTCGTCGAGAATCCAGCTGCCAACTTCAAGGATCAGGCCACTTTCTTCCAGTACCTGGATGAACTGCGAGGGCGGTTGCTGACCCAGTTGCGGGTGGTGCCAGCGCAACAGTACTTCGGCACCGACGATGCGGTTGTCGCGGGCATCCACCTGCGGCTGGAAGTGCAGTGCCAGCTCGCCACGGGCCAGGGCCAGGCGCAGGTCGTTTTCCATGCGCAGGCGCTCGCTGGCGGCCTTCTGCATGGTGGTGTGAAACAGCTGGGTGGTGTTGCGA  >  NC_002947/2246848‑2247139

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