Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,716,352 G→T 5.5% S528I (AGT→ATT)  PP_4173 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,716,3520GT 5.5% 239.6 / 3.4 73S528I (AGT→ATT) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base T (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATGAAG  >  NC_002947/4716208‑4716495
                                                                                                                                                 |                                                                                                                                                
gCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTc                                                                                                                                              >  6:41566/1‑149 (MQ=255)
  cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCCACCGGAGggtggt                                                                                                                                            >  4:267787/1‑149 (MQ=255)
  cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                            >  5:239000/1‑149 (MQ=255)
  cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                            <  8:237191/149‑1 (MQ=255)
  cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                            <  5:41566/149‑1 (MQ=255)
     ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTgg                                                                                                                                          <  4:101245/148‑1 (MQ=255)
     ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGa                                                                                                                                         <  1:111971/149‑1 (MQ=255)
          gACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGTCAATCCGGTCAACCGGAGTGTCGTGGAAGcc                                                                                                                                     <  1:192361/148‑1 (MQ=255)
          gACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcct                                                                                                                                            <  8:265988/141‑3 (MQ=255)
          gACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcct                                                                                                                                            >  7:265988/1‑139 (MQ=255)
           aCCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCAt                                                                                                                                  <  6:147408/149‑1 (MQ=255)
            ccTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCAtg                                                                                                                                 <  2:230270/149‑1 (MQ=255)
            ccTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCAtg                                                                                                                                 <  6:239000/149‑1 (MQ=255)
             cTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCAtgc                                                                                                                                >  2:13462/1‑149 (MQ=255)
                 aCCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCAtgctgc                                                                                                                             <  3:267787/148‑1 (MQ=255)
                  cccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGCCGCCGCCAACGACTGTACGCACGACGGTGCCGGACAATGAACAGGTACTGCTGGTGGAGGACAATCCGCTCAACCGGAGTGTCGTGGAAGCCATGCTgcg                                                                                                                            >  6:104827/1‑148 (MQ=255)
                   ccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTgcg                                                                                                                            >  5:7807/1‑147 (MQ=255)
                                   gCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAg                                                                                                           >  5:267052/1‑148 (MQ=255)
                                     ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                                        >  2:273807/1‑149 (MQ=255)
                                     ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                                        >  7:185014/1‑149 (MQ=255)
                                               ccTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                               >  4:233361/1‑148 (MQ=255)
                                                cTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                               >  7:187545/1‑147 (MQ=255)
                                                  gCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCagaa                                                                                             >  2:54493/1‑147 (MQ=255)
                                                  gCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCagaa                                                                                             <  1:54493/147‑1 (MQ=255)
                                                   cccccGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGAt                                                                                          >  7:171855/1‑149 (MQ=255)
                                                     cccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCa                                                                                                <  6:112865/141‑1 (MQ=255)
                                                     cccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCa                                                                                                >  5:112865/1‑141 (MQ=255)
                                                       cgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGcc                                                                                      <  6:30649/149‑1 (MQ=255)
                                                        gtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGcc                                                                                      >  1:216226/1‑149 (MQ=255)
                                                            cgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTg                                                                                   >  1:33326/1‑147 (MQ=255)
                                                                 aCGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGAACGGCAGAAGATGGCCTTGAAGCCGt                                                                            >  8:296671/1‑149 (MQ=255)
                                                                     cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                               <  7:92969/126‑1 (MQ=255)
                                                                     cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                               >  8:92969/1‑126 (MQ=255)
                                                                     cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAg                                                                        <  1:273807/149‑1 (MQ=255)
                                                                        cAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTg                                                                     <  1:13462/149‑1 (MQ=255)
                                                                                 tACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCt                                                                      <  7:206184/139‑1 (MQ=255)
                                                                                 tACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCt                                                                      >  8:206184/1‑139 (MQ=255)
                                                                                 tACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTcggccgg                                                            <  2:216226/149‑1 (MQ=255)
                                                                                    gaacgaCGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAg                                                         >  8:49491/1‑149 (MQ=255)
                                                                                        gacgGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCt                                                     >  5:313194/1‑149 (MQ=255)
                                                                                                ccGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTtgctg                                               <  7:49491/148‑1 (MQ=255)
                                                                                                ccGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTAGGCCGGCAGCGCTTtgctgc                                              >  3:145792/1‑149 (MQ=255)
                                                                                                   gACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGt                                           >  3:62459/1‑149 (MQ=255)
                                                                                                      aaTGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCt                                        >  8:215066/1‑149 (MQ=255)
                                                                                                      aaTGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTAGTGGAAGCCATGCTGCGCAGCCTGGGGCTCGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCt                                        >  3:150616/1‑149 (MQ=255)
                                                                                                        tGATCAGGTACTGCTGGTGGAGGACAATCGGGTCCACCGGAGTGTCGTGGAAGCCATGCTGCGCTGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                                      <  5:104827/149‑1 (MQ=255)
                                                                                                        tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                                      <  4:62459/149‑1 (MQ=255)
                                                                                                        tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGACAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGGAGCGCGTTGCTGCCGTGCGGa                                      >  4:4503/1‑149 (MQ=255)
                                                                                                          aTCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGACAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGAt                                     <  3:4503/148‑1 (MQ=255)
                                                                                                           tCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTgg                                                                                                              >  7:115663/1‑74 (MQ=255)
                                                                                                           tCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTgg                                                                                                              <  8:115663/74‑1 (MQ=255)
                                                                                                             aGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGa                                  <  6:313194/148‑1 (MQ=255)
                                                                                                             aGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGGCGTGGAAGCCATGCTGCTCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGAt                                 <  4:273427/149‑1 (MQ=255)
                                                                                                              ggTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGAtt                                <  2:294828/149‑1 (MQ=255)
                                                                                                              ggTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGGCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGAtt                                <  6:267052/149‑1 (MQ=255)
                                                                                                               gTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTg                               >  6:44595/1‑149 (MQ=255)
                                                                                                                  ctgctgGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTgccg                            <  3:233361/149‑1 (MQ=255)
                                                                                                                     ctgGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTgccgcc                          >  3:61887/1‑148 (MQ=255)
                                                                                                                             ggaCAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCAcc                 >  4:25959/1‑149 (MQ=255)
                                                                                                                                cAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTgg              <  8:171855/149‑1 (MQ=255)
                                                                                                                                cAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTACCTCACCTgg              <  4:61887/149‑1 (MQ=255)
                                                                                                                                 aaTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTggg             <  5:225765/149‑2 (MQ=255)
                                                                                                                                  aTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGAc            <  5:44595/149‑1 (MQ=255)
                                                                                                                                   tCCGGTTGACCGGATTGTCGTGGAAGCCTATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATgg                                   >  1:178842/1‑126 (MQ=255)
                                                                                                                                   tCCGGTTGACCGGATTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATgg                                   <  2:178842/125‑1 (MQ=255)
                                                                                                                                   tCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGGCGCCTGCCTGACCTGGACg           <  4:150616/149‑1 (MQ=255)
                                                                                                                                   tCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACg           >  7:120191/1‑149 (MQ=255)
                                                                                                                                      ggTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGc         >  1:149310/1‑148 (MQ=255)
                                                                                                                                       gTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTa       <  3:25959/149‑1 (MQ=255)
                                                                                                                                            ccGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATg     >  8:277746/1‑146 (MQ=255)
                                                                                                                                            ccGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATGaa   <  7:277746/148‑1 (MQ=255)
                                                                                                                                            ccGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATGAAg  <  8:185014/149‑1 (MQ=255)
                                                                                                                                            ccGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATGAAg  <  2:149310/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTATGAAG  >  NC_002947/4716208‑4716495

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: