Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,716,345 A→G 5.6% N526D (AAC→GAC)  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,3450AG 5.6% 215.9 / 2.0 73N526D (AAC→GAC) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (2/2);  total (34/39)
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.59e-01

GCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTA  >  NC_002947/4716199‑4716490
                                                                                                                                                   |                                                                                                                                                  
gCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAAc                                                                                                                                                  <  2:134853/149‑1 (MQ=255)
    tGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGa                                                                                                                                              <  5:377/149‑1 (MQ=255)
    tGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGa                                                                                                                                              <  7:293101/149‑1 (MQ=255)
    tGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGa                                                                                                                                              <  6:147806/149‑1 (MQ=255)
         gCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTc                                                                                                                                         >  6:41566/1‑149 (MQ=255)
           cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCCACCGGAGggtggt                                                                                                                                       >  4:267787/1‑149 (MQ=255)
           cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                       >  5:239000/1‑149 (MQ=255)
           cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                       <  5:41566/149‑1 (MQ=255)
           cTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGt                                                                                                                                       <  8:237191/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                                                                                                   >  7:185014/1‑149 (MQ=255)
                                              ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                                   >  2:273807/1‑149 (MQ=255)
                                                        ccTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                          >  4:233361/1‑148 (MQ=255)
                                                         cTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                          >  7:187545/1‑147 (MQ=255)
                                                           gCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCagaa                                                                                        <  1:54493/147‑1 (MQ=255)
                                                           gCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCagaa                                                                                        >  2:54493/1‑147 (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                                                                                                         <  8:115663/74‑1 (MQ=255)
                                                                                                                    tCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTgg                                                                                                         >  7:115663/1‑74 (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)
                                                                                                                                                   |                                                                                                                                                  
GCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACGGCTA  >  NC_002947/4716199‑4716490

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: