Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,196 G→T 5.3% A56D (GCC→GAC)  TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,1960GT 5.3% 312.5 / 5.5 95A56D (GCC→GAC) TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (47/43);  new base T (2/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

TTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAC  >  NC_002947/5297048‑5297338
                                                                                                                                                    |                                                                                                                                              
ttGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCggg                                                                                                                                                <  1:36229/149‑1 (MQ=255)
ttGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCggg                                                                                                                                                <  5:152094/149‑1 (MQ=255)
     aGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                           <  6:1950/149‑1 (MQ=255)
         cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAAc                                                                                                                                        <  1:204787/148‑1 (MQ=255)
         cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAAc                                                                                                                                        <  1:25067/148‑1 (MQ=255)
          gcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTg                                                                                                                                      <  7:220900/149‑1 (MQ=255)
              gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                                                                                  >  6:226638/1‑149 (MQ=255)
              gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                                                                                  <  3:34697/149‑1 (MQ=255)
                   tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATg                                                                                                                             >  1:370178/1‑149 (MQ=255)
                   tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATg                                                                                                                             >  5:254830/1‑149 (MQ=255)
                   tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGCAGTTGATCGCCCGCTCGATCAACTACaccgt                                                                                                                                >  5:155368/1‑142 (MQ=255)
                   tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGCAGTTGATCGCCCGCTCGATCAACTACaccgt                                                                                                                                <  6:155368/146‑5 (MQ=255)
                     gCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                           >  7:27560/1‑149 (MQ=255)
                     gCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                           >  3:51685/1‑149 (MQ=255)
                       ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTTCAGGTTATGGTTg                                                                                                                         >  6:315049/1‑149 (MQ=255)
                       ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         >  1:347550/1‑149 (MQ=255)
                        cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                        >  5:25253/1‑149 (MQ=255)
                         tCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                        >  5:287197/1‑148 (MQ=255)
                           gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGTGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCa                                                                                                                      >  7:9906/1‑148 (MQ=255)
                           gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                                                                                       <  2:370178/147‑1 (MQ=255)
                           gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCa                                                                                                                      <  1:65880/148‑1 (MQ=255)
                               ttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGa                                                                                                                  >  2:189086/1‑148 (MQ=255)
                                 gAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGaca                                                                                                                <  6:287197/148‑1 (MQ=255)
                                  aaGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGTGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACg                                                                                                              <  8:9906/149‑1 (MQ=255)
                                    gACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAAGGCCATCGATTCGTTGTCGGCGGCGCTGTGATCGAACACGACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATGAACTGCAGGTTATGGTTGGCATAGACACGCa                                                                                                            <  4:145010/149‑1 (MQ=255)
                                     aCCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAg                                                                                                           <  1:280749/149‑1 (MQ=255)
                                       cTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                         >  3:49131/1‑149 (MQ=255)
                                       cTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                         >  5:306172/1‑149 (MQ=255)
                                       cTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                         <  7:115351/149‑1 (MQ=255)
                                           gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                      <  5:226638/148‑1 (MQ=255)
                                           gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                      <  1:272719/148‑1 (MQ=255)
                                             cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                                   >  3:367215/1‑149 (MQ=255)
                                              tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                  >  4:46711/1‑149 (MQ=255)
                                                ggCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                      <  6:207924/143‑1 (MQ=255)
                                                ggCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                      >  5:207924/1‑143 (MQ=255)
                                                 gCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                                                                               <  4:367215/149‑1 (MQ=255)
                                                 gCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                                                                               >  4:34366/1‑149 (MQ=255)
                                                  ccACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                                                                               >  7:29363/1‑148 (MQ=255)
                                                    aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                           >  5:262856/1‑118 (MQ=255)
                                                    aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                           <  6:262856/118‑1 (MQ=255)
                                                       tcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGtt                                                                                                                                >  5:150392/1‑110 (MQ=255)
                                                       tcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGtt                                                                                                                                <  6:150392/110‑1 (MQ=255)
                                                                cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAggg                                                                                 <  4:95946/148‑1 (MQ=255)
                                                                cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                                <  7:30673/149‑1 (MQ=255)
                                                                cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                                <  3:46711/149‑1 (MQ=255)
                                                                   gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgag                                                                               <  2:347550/147‑1 (MQ=255)
                                                                   gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                                <  4:49131/146‑1 (MQ=255)
                                                                     aaTCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAc                                                                           >  6:88401/1‑149 (MQ=255)
                                                                      gccgccGCCAACGAATCGCTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAcc                                                                          >  4:348100/6‑149 (MQ=255)
                                                                           cGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCa                                                                     <  8:29363/149‑1 (MQ=255)
                                                                               aGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                  >  3:243351/1‑148 (MQ=255)
                                                                                  gATTCGTTGGCTGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAAc                                                              <  4:51685/149‑1 (MQ=255)
                                                                                   aTTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                             <  1:189086/149‑1 (MQ=255)
                                                                                     tCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACgg                                                            >  4:77757/1‑148 (MQ=255)
                                                                                           gcggcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaagcaa                                                     <  4:65545/149‑1 (MQ=255)
                                                                                            cggcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                    <  8:17994/149‑1 (MQ=255)
                                                                                             ggcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAgg                                                   <  3:77757/149‑1 (MQ=255)
                                                                                              gcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGCGCACGCCCAGCAGGGTGAGGGAGCTACCCGCCAGGCCAACGGCAAGCAAGGc                                                  >  3:44664/1‑149 (MQ=255)
                                                                                                 gcgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                                                                            <  8:181667/104‑1 (MQ=255)
                                                                                                 gcgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                                                                            >  7:181667/1‑104 (MQ=255)
                                                                                                    cTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCt                                            <  1:6/149‑1 (MQ=255)
                                                                                                     tGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTg                                           >  2:71545/1‑149 (MQ=255)
                                                                                                       tCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGTGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                         <  6:306172/149‑1 (MQ=255)
                                                                                                       tCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                         <  5:173813/149‑1 (MQ=255)
                                                                                                         aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                        <  6:80693/100‑1 (MQ=255)
                                                                                                         aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                        >  5:80693/1‑100 (MQ=255)
                                                                                                         aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                       >  1:95718/1‑149 (MQ=255)
                                                                                                         aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                       >  5:145301/1‑149 (MQ=255)
                                                                                                          tCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTg                                      >  3:333157/1‑149 (MQ=255)
                                                                                                           cGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTgg                                     >  8:131579/1‑149 (MQ=255)
                                                                                                               caccacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCAc                                 >  2:120864/1‑149 (MQ=255)
                                                                                                                accacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACg                                >  1:340828/1‑149 (MQ=255)
                                                                                                                 ccacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACgg                               <  7:131579/149‑1 (MQ=255)
                                                                                                                  cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                              >  3:70927/1‑149 (MQ=255)
                                                                                                                  cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                              >  7:170781/1‑149 (MQ=255)
                                                                                                                      gCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGCCGCTGAGGTGGCCACGGCCCAg                          >  2:297206/1‑149 (MQ=255)
                                                                                                                       ccGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                  >  5:89133/1‑108 (MQ=255)
                                                                                                                       ccGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                  <  6:89133/108‑1 (MQ=255)
                                                                                                                       ccGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGTCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAg                          >  3:14943/1‑148 (MQ=255)
                                                                                                                        cGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACCGCAAGCAAGGCGACGCTGAGGTGGCCACgg                               >  5:384526/1‑142 (MQ=255)
                                                                                                                        cGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACCGCAAGCAAGGCGACGCTGAGGTGGCCACgg                               <  6:384526/142‑1 (MQ=255)
                                                                                                                         gCTTATACGGTGTAGTTGATCGAGCGGGCGATCAAATGCAGGTTATGTTTGGCATAGACACGAAGGGCGAGCACGCCCAGCAGGGAGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCATCAc                       >  4:353608/1‑149 (MQ=255)
                                                                                                                                cGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAgc                  <  2:95718/147‑1 (MQ=255)
                                                                                                                                  gtgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAg              >  5:95840/1‑149 (MQ=255)
                                                                                                                                   tgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCGCGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg             <  7:110539/149‑1 (MQ=255)
                                                                                                                                   tgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg             <  3:34366/149‑1 (MQ=255)
                                                                                                                                   tgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg             >  5:173906/1‑149 (MQ=255)
                                                                                                                                       gTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCg         >  7:18719/1‑149 (MQ=255)
                                                                                                                                       gTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCg         <  6:375790/149‑1 (MQ=255)
                                                                                                                                        ttGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg        >  4:92514/1‑149 (MQ=255)
                                                                                                                                        ttGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg        <  6:362529/149‑1 (MQ=255)
                                                                                                                                         tGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCggg       >  7:131444/1‑149 (MQ=255)
                                                                                                                                           aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgcg     >  4:307303/1‑149 (MQ=255)
                                                                                                                                              gAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAc  <  4:333157/149‑1 (MQ=255)
                                                                                                                                              gAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAc  <  6:95840/149‑1 (MQ=255)
                                                                                                                                               aGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg             >  5:300679/1‑137 (MQ=255)
                                                                                                                                               aGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg             <  6:300679/137‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAC  >  NC_002947/5297048‑5297338

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